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Abstract 



X-ray observations of active galactic nuclei and black hole binaries offer a unique lab- 
oratory for testing the general relativity in strong gravity regime, for studying accretion 
physics around black holes, and for constraining properties of accreting black holes. In 
this Thesis, we discuss black hole spin measurements employing the relativistic iron line 
profiles in the X-ray domain. 

We investigate the iron line band for two representative sources - MCG -6-30-15 (ac- 
tive galaxy) and GX 339-4 (X-ray binary). We compare two relativistic models of the 
broad iron line, LAOR and kyrline. In contrast to LAOR, the kyrline model has the 
spin value as a variable parameter. However, the LAOR model can still be used for evalua- 
tion of the spin if one identifies the inner edge of the disc with the marginally stable orbit. 
We realise that the discrepancies in the results obtained with the two models are within 
general uncertainties of the spin determination using the skewed line profile when applied 
to the current data. This implies that the spin is currently determined entirely from the 
position of the marginally stable orbit while the effect of the spin on the overall line shape 
would be resolvable with higher resolution X-ray missions like IXO (International X-ray 
Observatory). 

We show that the precision of the spin measurements depends on an unknown angular 
distribution of the disc emission. Often a unique profile is assumed, invariable over the 
entire range of radii in the disc and energy in the spectral band. However, an improper 
prescription for the directionality profile affects the parameters inferred for the relativistic 
broad line model. We study how sensitive the spin determination is to the assumptions 
about the intrinsic angular distribution of the emitted photons. We find that the uncertainty 
of the directional emission distribution translates to 20% uncertainty in the determination 
of the radius of marginally stable orbit. 

By assuming a rotating black hole in the centre of an accretion disc, we perform 
radiation transfer computations of an X-ray irradiated disc atmosphere (NOAR code) to 
determine the directionality of outgoing X-rays in the 2-10 keV energy band. Based 
on these computations, we find that from the simple formulae for the directionality, the 
isotropic case reproduces the simulated data with the best accuracy. The most frequently 
used limb darkening law favours higher values of spin and, in addition, a steeper radial 
emissivity profile. We demonstrate our results on the case of XMM-Newton observation 
of MCG -6-30-15, for which we construct confidence levels of chi-squared statistics, and 
on the simulated data for the future X-ray IXO mission. 

Furthermore, we present a spectral analysis of an XMM-Newton observation of a 
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Seyfert 1.5 galaxy IRAS 05078+1626 being the first X-ray spectroscopic study of this 
source. The lack of the significant relativistic blurring of the reflection model component 
suggests the accretion disc to be truncated at a farther radius (inner disc radius Rin > 
60Rg). 

As a by-product of our reduction of the XMM-Newton data, we find that a careful 
treatment of the raw instrumental data is necessary to obtain the highest quality data. As a 
crucial step, we consider the correct re-binning of the data reflecting the energy resolution 
of the used instrument. Photon pile-up is another problem which may occur in the data of 
very bright sources, and it may significantly affect the spectral shape. 

In summary, we found relativistic iron line models to be a feasible method for mea- 
suring the spin of black holes at all scales - from solar-mass microquasars to giant black 
holes of billions solar masses in distant quasars. Some useful constraints on spin are 
achievable already from X-ray spectra of currently operating instruments. However, our 
simulations with the tentative IXO response show a significant improvement in the accu- 
racy of spin measurements in the future. 

Keywords: black holes - accretion disc - active galaxies - Galactic X-ray binaries 
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Preface 



In questions of science the authority of a thousand is not worth the humble 
reasoning of a single individual. 

Galileo Galilei, 1632 

Black holes are objects which were first created in minds of theoretical astrophysicists, 
and for a long time, they were not supposed to exist in the real Universe, in the endless 
world of stars and planets. Even Albert Einstein who developed the theory, which allowed 
the existence of black holes, did not believe that the nature would be so crazy to give 
permission to such objects to form. However, theoretical astrophysicists calculating the 
details of the stellar collapse, such as Chandrasekhar, Oppenheimer, Snyder etc., predicted 
that the collapse of too massive stars could not be stopped by any means and that it must go 
on to create space-time singularities, which were later named by Wheeler as black holes. 
Meanwhile, with no idea of connection, astronomers classified radio bright galaxies as a 
peculiar group distinguished from the standard galaxies by extremely strong radio power. 
Presently, we generally believe that accretion on a rapidly rotating super-massive black 
hole is the process behind such an enormous power. 

The first suggestions of black holes as real celestial objects came with the development 
of X-ray astronomy in 1960s. Riccardo Giaconni, one of the pioneers of X-ray astronomy, 
won the Nobel prise in physics in 2002 for opening the X-ray window to the Universe. 
The discoveries of pulsars, quasars and X-ray binaries started a fruitful life of a new 
branch of astronomy, astronomy of ultra-compact objects. Since that time, black holes 
have fascinated many people around the world, not only scientists. Studying these objects 
brings light to the physics of stellar collapse, galaxy formation, accretion physics and 
behaviour of matter in the strong gravitational field. 

Astrophysical black holes are actually very simple objects being characterised only by 
their mass and angular momentum (spin). In the current knowledge, it seems that there 
are two populations of black holes according to their mass - stellar black holes of mass of 
several solar masses and super-massive black holes of mass of millions to billions solar 
masses. The distribution of their spin is still unknown. Measuring the spin is difficult 
because the effect of the spin quickly decreases with the growing distance from the black 
hole. Matter feels the black hole spin only within several gravitational radii. However, the 
black hole spin plays an important role in black hole energetics and evolution. The spin is 
assumed to be responsible for generation and up-keeping of the powerful relativistic jets. 
The information about the spin value on a statistically significant sample of black holes is 
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important in the understanding of formation and growth of black holes. It can significantly 
help to answer the question if the observed spin value is natal or if the rotating black holes 
are spun-up via accretion. 

The innermost parts of black hole accretion discs can be uncovered in high energetic 
radiation, such as X-rays or y-rays. This is for two reasons. First, such an energetic radi- 
ation can originate only under the extreme conditions close to a black hole. Second, any 
weaker radiation is often efficiently absorbed by surrounding matter and cannot reach a 
distant observer. Fast development of X-ray detectors in the last two decades has allowed 
astronomers to provide spectra with an unprecedented sensitivity, and so constrain accre- 
tion flows within a few gravitational radii and measure the black hole spin. There are 
currently several methods of spin measurements. Modelling of relativistic iron line repre- 
sents one of them, and it is particularly suitable because it is applicable to black holes at 
all mass scales. 

In this Thesis, we will look into the inner black hole accretion disc with relativis- 
tic models of iron line. In Section [2l we summarise the basic concepts of this method, 
and compare two models which employ diff"erent ways to determine the spin value. In 
Section [3l we study in detail the role of the angular emissivity on the spin value mea- 
surement. Both analyses are provided on the current XMM-Newton data whose reduction 
is described in Section HI and on the simulated data of a next generation X-ray mission 
with a significantly higher energy resolution. In Section |4l we also present the results of 
our X-ray spectroscopic study of the Seyfert 1.5 galaxy IRAS 05078-1-1626, which repre- 
sents the first X-ray spectroscopy measurement of this source. The achieved results are 
discussed at the end of each Chapter. Main conclusions of the Thesis are summarised in 
Section [5] Future perspectives are pronounced in Section [6l 

Most of the results pre sented in the Thesis were published in two papers in Astro nomy 
and Astrophysics ioi irnal (ISvoboda g? a/.ll2009l: ISvoboda. Guainazzi & Karasll2010|) . and 
several proceedings (jSvoboda et a/.ll2008alJbl : ISvobodall2009[) . 



In Prague, April 2010 



Jin Svoboda 
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Chapter 1 
Introduction 



1.1 Black holes 



Black holes are regions of space-time in which the gravitational well is so deep that no 
particle or even light can escape. First ideas of existence of such objects from which a 
distant observer cannot get any s i gnal d ue t o the escape vel ocity greater than speed of 
light were formulated by Mi chelll (|1784h and iLaplacd (117991) (recent English translation 
can be found in Stephani 2003'). However, real f oundations for b lack hole theory were laid 



19161). The first exact so lution 



by discovery of general theory of relativity by Einstein 

of Einstein's equations describing a black hole was found by iSchwarzschildl (| 1 9 1 61) for a 

point mass assuming spherical symmetry. The radius of the event horizon is accordingly 

called Schwarzschild radius: 

2GM 
r, = 2rg = ^, (1.1) 



where r^ is the gravitational radius which is equivalent to the mass if the geometrised units 
with c = G = 1 are used. The Schwarzschild metric in Schwarzschild coordinates is: 



ds^ = -\\ - ^\df + 



1 



('-?) 



dr^ + ?dQ}, 



(1.2) 



where d^ = dO^ + sin^ Odip^, t is the time measured by an observer at infinity at rest, r is 
the radial coordinate, 6 is latitude, and ip is the azimuthal angle. 



The solution for rotating black holes was found almost half a century later by iKen 



(| 19631) . The Kerr's metr ic in the Boyer-Lindquist spheroi dal coordinates {t, r, 6, (f) and 
geometrised units is (e.g. IMisner. Thome & Wheeleilll973L chap. 33): 



, o AZ , , A sin , , , o 2,9 ^ .„9 

A Z 



{dip - codty + -dr^ + I.de^, 



(1.3) 



where the metric functions are: A(r) = r^ - 2r + a^, S(r, 6) = r^ + a^ cos^ 6, A{r, 6) 



ir^ + a^\ - A(r) a^ sin^ 9, and a)(r, 9) = 2arjA{r, 9); a denotes the specific rotational an- 
gular momentum (spin) of the central body. 
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The conversion factor from the angular momentum 7phys (in physical units) to the 
angular momentum J (in geometrical units) reads: J = (G/c^)7phys. The geometrised 
dimension of J is the square of the length [cm^]. It is convenient to make all geometrised 
quantities dimensionless by scaling them with the appropriate power of mass M. The 
dimensionless specific angular momentum, a = J/M^, spans the range -1 < a < 1, where 
the positive/negative value refers to the motion co/counter-rotating with respect to the (f- 
coordinate. We will further assume co-rotational motion only (a > 0). The magnitude of 
a is thought to be less than unity in order to have a regular horizon and avoid the case of 
naked singularity. The value of the outer horizon is: 



1/2 



r+ = 1 + (l - ay . 



(1.4) 



Fo r the Keplerian angular velocity of the orbital motion, we obtain (IBardeen. Press & Teukolsky 
I972I) : 

"kW = ^7^- (1-5) 

r^'^ + a 

For the linear velocity with respect to a locally non-rotating observer, we have: 



^(r) 



Ai/2 (^3/2 + ^) 



(1.6) 



o 
o 

> 

!a 

o 



0.8 



0.6 



0.4 



0.2 




Figure 1.1: Orbital velocity "Vir) of co-rotational motion near a rotating black hole, as 
given by formula (|1.6I) for three values of the black-hole dimensionless angular momen- 
tum parameter a (spin). The thick coloured part of each curve indicates the range of radii 
above the marginally stable orbit, r > r^sici), where the circular motion is stable. The thin 
curve indicates an unstable region at small radii. 



1.2. Marginally stable orbit 



The velocity at the marginally stable orbit reaches a considerable fraction of the speed 
of light c (=1) and has a similar value 'V{r) ^ (0.5 - 0.6) c for any value of the angular 
momentum, see Figure fTTTl For large spin, a small dip develo ps in the velo city profile near 



the horizon. Although it is an interesting feature (see Stuchli k et al.\\2005\) . its magnitude 



is far too small to be recognised with current observational facilities. 

The Schwarzschild and Kerr metric represent solutions of Einstein's equations rele- 
vant for astrophysical black holes which are assumed to be electrically neutral. Never- 
theless, they can be further generalised by taking of the electric charge into the consid- 
eration. The corresponding solutions are Reissner-Nordstrom metric for a non-rotating 
charged black hole, and Kerr-Newman metric for a rotating charged black hole, respec- 
tively. Besides the mass, angular momentum and electric charge, black holes do not 
have any other parameters, which is often called "no-hair theorem" (this statement gets 
its name from a comment by the famous astrophysicist John A. Wheeler 1968, see also 



Misner. Thome & Wheeled (llQVSh ). 



1.2 Marginally stable orbit 

The marginally stable orbit (Vs), sometimes also called innermost stable circular orbit 
(ISCO), is the closest orbit to the centre of a black hole where the orbit of a test particle is 
stable. Below this orbit, only unstable or unbound orbits coming from infinity may occur. 



The position of ISCO depends on the value of the spin (|Bardeen. Press & Teukolsky 



19721) : 

= 3+Z2-[(3-Zi)(3+Zi+2Z2)]^ (1.7) 



' ms 



Zi = 1 -I- (1 - a^)3[(l -I- a) 3 -I- (1 - ap] and Z2 = (3a^ -I- Z^p. Notice that r^nsia) spans the 
range of radii from Vs = 1 for a = I (the case of a maximally co-rotating black hole) to 
''ms = 6 for a = (static black hole). Figure [L2] illustrates the relation (11.71) graphically. 
It is generally supposed that the rotation of the astrophysical black holes is li mited by 



an eq uilibrium value, a = 0.9982, because of capture of photons from the disc (iThome 



19740 . This translates to Vs - 1-23. For a hypothetically higher value of the spin than 



a = 1, the radius of th e marginally stable orbit decreases to Vs = 2/3 and then increases 



again (iStuchlikll 19801) . 

The ISCO is an important quantity in the standard accretion disc theory because the 
inner edge of the accretion disc is assumed to coinc ide with it. However, this rnay no t 



be satisfied precisely under realistic circumstances (Beckwith. Hawley & Krolikll2008|) 



The magnitude of t he resulting error on spin m easurements (see Section [T. 5. 21) was con- 



strained recently by [Reynolds & FabianI (|2008r) who applied physical arguments about the 



emission properties of the inner flow. It is very likely that this discussion will have to con- 
tinue for some time until the emission properties of the general relativistic MHD flows are 
fully understood. 
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Figure 1.2: Relation between spin a and marginally stable orbit r„ 



1.3 Black hole accretion 

1.3.1 Accretion power 

Black holes are interesting objects for the theoreticians since they represent a natural lab- 
oratory for testing the theory of general relativity, but they are also important objects for 
astronomers due to their interaction with the surrounding matter. Nowadays, black holes 
cannot be observed directly as the sensitivity of gravitational wave detectors is still in- 
sufficient. However, there is an increasing amount of observations of the electromagnetic 
radiation of the matter accreting on black holes. These observations can constrain the 
black hole parameters as well as the accretion physics. 

The accretion power is due to conversion of the gravitational potential energy into 
radiation. For a body of mass M and radius R we can estimate the energy released by the 
accretion of a particle with mass m as: 



A£ = 



GMm 
R 



(1.8) 



We can compare the power of the accretion process with the energy output by nuclear 
fusions. In the case of hydrogen burning, we obtain AE ^ 6x lO^^erg g"^ For a neutron 
star with M « Mq and R « 10 km, the energy output by accretion on the stellar surface 
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is AEacc ~ 10^°ergg"^ This means that more energy per a mass unit is released by the 
accretion on a compact object than from the nuclear synthesis. Accretion onto super- 
massive black holes in quasars powers the most luminous sources in the Universe. 

The maximal luminosity produced by the spherical accretion on a star is restricted by 

the Eddington limit: 

AnGMm„c 
Lm^ = -, (1.9) 

where nip is the mass of a proton (i.e. a nucleus of hydrogen) and crj is Thomson cross 
section. The Eddington luminosity represents the highest possible luminosity on the con- 
ditions of spherically symmetric accretion M and totally ionised accreting material. Li 
this case, the opacity is dominated by electron scattering and we suffice only with Thom- 
son cross section, u For accretion powered objects the Eddington limit implies a limit on 
the steady accretion rate. 

If all the kinetic energy of infalling matter is conversed to radiation at the surface of 
the central body, then the accretion luminosity Lace is given by: 

GMM 

^acc = —^. (1.10) 

In the case of black hole accretion, the radius R does not refer to the surface but a natural 
choice is a Schwarzschild radius. The uncertai nty of the value of the luminosity can b e 



parametrised by a dimensionless efficiency rj (|Shakura & Sunyaevlll973l : lSoltanlll982|) . 
which measures how efficiently the rest mass energy of the accreted material is converted 

to radiation: 

GMM . , 

L,,. = ri——^riMc\ (1.11) 

K 

Comparing eq. (11.111 ) with energy released by the b urning of hydrogen we get 77 = 0.007 



for nuclear synthesis. IShakura & SunyaevI (|1973 i) mentioned rj x 0.06 in the case of 



Schwarzschild's metric, and in Kerr's metric 77 can achieve 40%. The estimation of real- 
istic value for rj in black hole a ccretion is an irnportant problem, a reasonable guess for 



it would appear to b e 77 « 0.1 (|Yu & Tremaine 



Marconi et a/.ll2004|) 



2OO2I : lElvis. Risaliti & Zamoranil 12002 



1.3.2 Accretion discs 

The dominant accretion process for compact objects involves disc accretion. Most of 
the accreting matter, which is gas supplied by a donor star in binary systems or the host 
galaxy in quasars, possesses sufficient angular momentum to go into the orbit around the 
black hole, forming an accretion disc. The previous relations calculated for spherically 
symmetric case will serve as convenient approximations and estimations. Accretion disc 
physics includes many processes, including gravity, hydrodynamics, viscosity, radiation 
and magnetic fields. The angular momentum of matter in an accretion disc is gradually 



'The radiation influences mostly electrons which are of the less mass than protons, vice-versa the gravity 
more attracts the massive particles. Electrons and protons hold together due to electromagnetic interaction. 
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transported outwards by stresses (turbulent, magnetic, etc.). The time scale for redistribut- 
ing angular momentum is long compared with either the orbital or radiative time scales, 
which allows matter to gradually spiral inwards. 



Steady thin accretion discs (Shakura-Sunyaev) 



Steady thin accretion discs represent the standard disc solution found by lShakura & Sunyaev 



(|1973[) . The model is applicable when the characteristic values of luminosities are sub- 
Eddington, L/L^dd ^ 1, i-C- when the vertical component of the thermal disc radiation 
alone cannot support the matter against gravity at substantial altitudes above the disc 
plane. This condition may be violated if the thermal motion is predominant. The thermal 
criterion for the Shakura-Sunyaev thin disc is: 

^^«l. (1.12) 

where k^ = 1.38 x 10^^ JK^ is the Boltzmann constant, T is the gas temperature, and 
fi is the mean mass per particle. When the condition (11.121 ) is not fulfilled, luminosities 
may reach Eddington values, L/L^dd * 1, and the radiation-pressure force becomes com- 
parable to that of the gravity. The height of the disc ceases to be small and the thick-disc 
solution needs to be considered. 

When the assumption of geometrical thinness is justified, only radial advection is 
dominant and the equations for the conservation of mass and angular momentum Q can 

be written as: 

5S d 
r— + —{rllv,.) = 0, (1.13) 



and 



where 



dt dr 



^ (Sr^^) + -(Sv.r^Q)) = — — , (1.14) 

ot or 2n or 



-, dQ 
Gir,t) = 2nr\X—- (1.15) 

or 



defines the torque, S is the surface density, v^ is the velocity in the radial direction, and v 
is the kinematic viscosity. Combining eqs. (11.131) - (|1.15l) and using eq. (11.51) we obtain 
the equation governing time evolution of surface density in the Keplerian disc: 

f = -|-(V^|-(v2:V^)). (1.16) 

ot r or or 

The kinematic viscosity v may be a function of local variables in the disc. The 
eq. (11.161) has the form of non-linear diffusion equation governing the behaviour of 2(r, t). 
Given a solution for 2(r, t) the radial velocity is: 

v, = ----(vSV^). (1.17) 



1.3. Black hole accretion 



Finally, we need some prescription for v to close the system of equations and to fully 
determine the radial structure of the accretion disc. All of the qualities currently ignored, 
such a detailed micro-physics, enter into the problem via v. After radial integration of the 
equations (11.131) and (11.141) and assuming time-steady discs (^ = 0), we get: 

rSvr = const. (1-18) 

and 






i-'y' 



(1.19) 



At r = r,„, the viscous torque G{r) vanishes. It is generally assumed that this radius co- 
incides with the marginally stable orbit because below it the matter losses the centrifugal 
support. 

The essential idea of the accretion process is dissipation of energy. The dissipation is 
caused by viscous torques and is given by: 



D(r) 



dr 



Using eq. (11.151) we obtain: 



1 / da\ 



For the Keplerian steady disc, we get the relationship for D{r) independent of v: 

3GMM 



Dir)^ 



%nr^ 



i-ly'^ 



(1.20) 



(1.21) 



(1.22) 



The luminosity produced by the disc between radii ri and r2 is given by: 



L(ri 



,r2)=2 r 



D{r)2nrdr, 



(1.23) 



and using eq. (11.221) : 



L(ri,r2) = 



3GMM n 



I-ItF 



dr 



(1.24) 



Setting ri = r;„ and r2 -^ oo, we obtain the luminosity for the whole disc: 

GMM 1 



■^disc 



2r;, 



= ^L„ 



(1.25) 



The thin disc is characterised by no motions or accelerations in the z-direction and the 
relevant Euler's equation has the form: 



Idp _ d 

p dz dz 



GM 



(r2 + z2)' 



(1.26) 
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where p is the pressure. For a thin disc (z «; r) we can write: 

Ip GMh 



p h r 



3 



(1.27) 



where h is the typical vertical scale height of the disc. From eq. (11.271) we can estimate 
the thickness of the disc as: 

/^ « ^ = — , (1.28) 

where vk is the local Keplerian velocity and c^ = J- is the sound speed. Hence, we can 
write the criterion for the thin disc as: 

c.,«Vk. (1.29) 

It means that the accretion disc is geometrically thin when the Keplerian velocity is highly 

supersonic. 

The total pressure of the disc is the sum of gas and radiation pressure. The equation 

of state has the form: 

pkriT Act a 
P = t-^ + —T\ (1.30) 

/imp 3 c 

The relationship for the central temperature can be derived from evaluating the heat 
loss per unit area by the radiative transport with the thermal energy by viscous dissipation 
given by eq. (11.221) : 



4cr^4 ^. X 3GMM 
— T^ = D(r) 



1-1^ 
r 



3r Snr^ 

where r is the optical depth of the disc, which can be defined as: 



(1.31) 



T = KrpH = /crE, (1.32) 

where a-r is the total Rosseland mean opacity. For hot discs around compact objects, 
free-free transitions and Thomson scattering contribute mostly to the opacity. The above 
estimation of the central temperature value is valid only for the case of an optically thick 
disc with T » 1 . 

Up to now, we have considered equations of the mass conservation, the angular mo- 
mentum conservation, energy conservation, hydrostatic equilibrium, equation of state, 
and equation of radiative transport for the thin steady disc. In order to study the detailed 
physical structure of such a disc, or any aspect of its time-variability, as well as stabil- 
ity, the knowledge of the viscosity v is required . A simple and very useful model of the 



viscosity is a-model (IShakura & Sunyaevlll973|) 



V = acsh, (1.33) 

where cr is a dimensionless parameter, which can take the value from to 1 . In the first 
approximation, a is a constant for a given accretion disc. 



1.3. Black hole accretion 



The viscous stress /^ exerted in the v'-direction can be defined for a Keplerian disc as: 



L = i^riQ.K, 



(1.34) 



where 77 is the coefficient of dynamic viscosity. For turbulent motion, it can be expressed 
as: 

?7 ~ PVturb^turb, (1.35) 

where p is the matter density, Vturb < Cs is the velocity of turbulent cells relative to the 
mean gas motion and lt„rh < h is t he size of the largest turbulent cells (|Landau & Lifshitz 
1959l : IShapiro & Teukolskyi 1983b . Using eq. (11.281) we can estimate the value of viscosity 



stress as: 



In general, we may write: 



L < ipC.h) ^K ~ PC,^ = P. 



U = <^^' 



(1.36) 



(1.37) 



where a is the same parameter as in eq. (11.331) . From eq. (11.171) . the radial velocity can be 
expressed by this model as: 

V, «-«^^«c,. (1.38) 

r r 

Thus, the radial inflow is very subsonic. 

The steady thin disc solution (or also called Shakura-Sunyaev's solution) allows to 
express the central density p{r), the surface density 2(r), the central pressure P(r), the 
disc height hir), the radial drift v,(r), the central temperature T{r), and the optical depth 
T(r) as functions of the parameters M, M, and a. 



Accretion flows in sub- and super-Eddington regime 

The condition for the steady thin disc accretion is not always accomplished. The ac- 
cretion rate is the main factor which constrains the shape of the accretion flow. In the 
low (sub-Eddington) accretion regime, the cooling mechanism via rad iation ceases to 



be sufficiently efficient and advection mecha n ism takes place instead (llchimarul 11977 



iNaraya n &Yil 1 19941 : lAbramowicz et al\ 1 19951 : iNarayan & McClintocM I2OO8I and refer- 
ences therein). In the super-Eddington regime, the gas pressure is strong enough to ex- 
pand sig nificantly the disc vertically and the height of the d isc needs to be taken into 
account (JAbramowicz et a/.|[l988l : lYuanlllOOll : ISadowskill2009|) . 

Generally, we may write the energy conservation equation per unit volume as: 



dS 



(1.39) 



where p is the density, T is the temperature, 5 is the entropy per unit mass, t is the time, 
and q^ and q- are the heating and cooling rates per unit volume. Since all the entropy 
stored in the gas is advected with the flow, the left-hand of eq. (11.391 ) may be replaced 
by a quantity ^adv which represents cooling rate via advection. The heat energy released 
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by viscous dissipation is partially lost by radiative cooling q^ and partially by advective 
cooling ^adv: 

q+ =q-+q^dy. (1-40) 

For the Shakura-Sunyaev thin disc is q » ^adv If, oppositely, q^ «c q^^v, the gas is 
radiatively inefficient and the accretion flow is under-luminous (L «c O.lMc^). Such an 
accretion flow is known as ADAF (advection-dominated accretion flow), or RIAF (radia- 
tively inefficient accretion flow). 

There are two distinct regimes of advection-dominated accretion flow. The first one 
occurs when the cooling time scale is much larger than the accretion time scale, t^ooi » 
tac,c.- It is the standard case of the ADAF, a self-similar solution of which is described 



by iNarayan & Yil (|1994|) . The second regime corresponds to very high scattering op- 



tical depth when the radiation is unable to diff'use out of the system, i.e. the photon 
diff'usion time is much larger than the accretion time sc ale, ?diff ^ ^acc- This radiation- 
trapped regime was briefl y discussed by lBegelmanl (Il979() and then developed in detail by 



Abramowicz et al.\ (|l988|) as the "slim disc" model. 



Accretion state of an accretion disc evolves in time because the accretion rate varies 
and the accretion flow regimes interchange accordingly. We know from the observations 
that the black hole binaries change the s pectral properties rapidly and s everal different 
X-ray states were defined (for review see 



Remillard & McClintock 2006. and references 



th erein). A close connection between th e states and accretion flow regimes was discussed 



by lEsin. McClintock & NarayanI (119971) . The value of the accretion rate is the main cri- 
terion for appearance of a certain X-ray state. Their schematic sketch is shown in Fig- 
ure II. 3[ which illustrates the link between the accretion states and the mass accretion 
rate. However, this simple picture is rather far from the complete description of the spec- 
tral states behaviour, not taking into account dynamic properties of the corona including 
strong flares, hysteresis of the states when interchanging etc. 

The main conclusion of the unification scheme is that below some critical value of 



the accretion rate, such as merit ~ 0.08 niEdd ( Esin. McClintock & NarayanI 119971) . the 



accretion disc does not extend to the innermost stable orbit, but, instead, it is truncated 
at some further radius whose position is inversely proportional to the accretion rate. The 
standard thin (Shakura-Sunyaev) accretion disc occurs only in the high state. In the very 
high state, the radiation pressure inflates the disc in the close neighbourhood to the black 
hole, and consequently, the advection plays again an important role in the inner flows 
(slim disc solution). 

Thermal radiation of accretion discs 

The radiation temperature follows from thermodynamical considerations: 

kT,,a ^ hv, (1.41) 

where h = 6.626 x 10"^"^! s is the Planck constant and v is the mean frequency. The value 
of this temperature can be estimated by comparing it with the black-body temperature T\, 
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Figure 1.3: Schematic sketch of the configuration of the accretion flow in difi'erent spec- 
tral states shown as a function of the mass accretion rate in Eddington units. The ADAF 
is indicated by dots. The thin disc (horizontal bars) extends to the innermost stable orbit 
only in the high state. In lower states, the truncation of the disc occurs, and the transition 
radius depends on the accretion rate. Figure is adopted from lEsin. McClintock & Narayan 
(Il997h . 
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and the thermal temperature Tn^: 



L V 



■^acc 



UnR^ 



(T 



(1.42) 



GMm„ 
T. = -^, (1.43) 

where cr = 5.67x10"* J s~' m"^ K""* is the Stefan-Boltzmann constant. These temperatures 
are two extreme cases for Trad- If the accretion flow is optically thick, the radiation reaches 
thermal equilibrium before escaping and r^ad « T'b- On the other hand, if the conversion 
to radiation is direct (from optically thin material), then r^d ~ T'th- It gives the limits for 
the radiation temperature: 

n < Trad < r,h. (1.44) 

For a solar-mass neutron star we will give the range for photon energy: 

IkeV < hv < 50MeV. (1.45) 

Thus, we can expect accreting neutron stars and black holes to appear as X-ray emitters, 
or possibly y-ray sources. 

If we suppose that the viscously dissipated energy (eq. 11.221) is radiated as a black- 
body spectrum (Shakura-Sunyaev disc), using eq. (11.421) we get the relationship for the 
surface temperature of the disc: 



T{r) 



_ (3GMM 



l-(^ 
r 



(1.46) 



For a fixed ratio between the source luminosity and Eddington limit and a scaled radius 
r/M the temperature of the disc depends on the mass as: 

r(r)ocM-^. (1.47) 

The temperature of the innermost region of an accretion disc surrounding a solar-mass 
black hole isT ^ 10^ K. Using eq. (11.411) . the corresponding spectral energy is E x I keV. 
Thus, the thermal component of the black hole accretion disc in an X-ray binary occurs 
in the soft X-rays. For a super-massive black hole, the accretion disc temperature is 
T X 10^ - 10^ K for the mass range M « 10^ - 10^ Mq. The thermal disc component has 
its maximum in the ultraviolet energy band and the spectral energy is £ < 0.1 keV for 
M « 10^ Mq, and E < 0.01 keV for M « 10'^ Mq. 

Time dependent discs 

Accretion discs are fueled by a variable amount of the accreting matter, and thus, the 
time variability of accretion discs is a natural consequence. Let us estimate diff"erent time 
scales: dynamical, thermal and viscous. The dynamical (or orbital) time scale can be 
defined as: 

t,yn = 0.-'. (1.48) 
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The thermal time scale, on which the dissipated energy is radiated from the accretion disc, 
is given by: 

^th « ^ = 4^. (1.49) 

a ail 

The viscous time scale which is linked with the redistribution of the angular momentum 

is defined as: 



^2 ,^.2 



Comparing the time scales we get for a thin disc: 



» tth > td- 



yn- 



(1.50) 



(1.51) 



On dynamical time scales we can well consider the temperature and a parameter in- 
dependent of time. Although the a-prescription is very successful in many applications, 
there are some limitations of this approach. Especially, the presence of the magnetic field 
plays an important role in the accretion disc stability. The present-day theoretical mod- 
els assume that the turbulence, which is respon sible for the angular mom entum transport, 
may be due to a magneto-rotational instability (IBalbus & Hawleylll99l|) . 



1.3.3 Accretion disc atmosphere and disc reflection 

X-ray emission of both, black hole binaries and active galaxies, is characterised by a 
power-law component with an exponential cutoff at high energy (E ^ 300 keV). Its origin 
is suggested to be due to multiple inverse Compton scatterings of the "seed photons" 
from the accretion disc (UV photons in the case of AGN, soft X-ray photons in the case 
of black hole binary) i n the optically thin accret ion disc atmosphere, so called "corona" 
(|Thorne & Pricdl 19751 : iHaardt & Maraschilll99lb . The corona is believed to consist of hot 
relativistic electrons which are possibly heated up by the magnetic dissipation processes. 
These processes may be caused by amplification of the magnetic field due to convective 
motions and differential rotation within a hot inner region of th e accretion disc, resulting 
into flaring events in the places of magnetic re-connections (IGaleev. Rosner & Vaiana 
19791 : iHaardt. Maraschi & Ghisellinill 19941 : ICzemv & Goosmannll20ol 

Geometrical properties of the corona are still uncertain, as well as the distribution of 
the Comptonising electrons (t hermal, non-thermak or rn ixed). The corona might be "sand- 
wiching" the accretion disc (IHaardt & Maraschilll99 1'). or rather be locally centralised 

with a typical size of a few of r^. In this geometry ("sphere -i- disc" ), the corona is irradi- 

ated by soft photons from the cooler outer parts of the accretion disc (IShapiro. Lightman & Eardley 
19761 : IHaardt. Maraschi & Ghisellinill 19941 : Istem et a/.lll995l : lDoveg?a/.lll997h . 

Some X-ray photons produced in the corona may escape directly to the observer and 
then be detected as the primary power-law radiation, but some of the m may illuminate the 
disc a nd be reflected from its surface before reaching the observer (Basko. Sunyaev & Titarchuk 
1974|) . The illuminating radiation is partly absorbed in the disc medium and partly re- 
radiated from the accretion disc. The reprocessed spectrum is characterised mainly by the 
Compton hump at E > lOkeV and the fluorescent iron line at E > 6 - 7keV, as seen in 
Figure [T31 or Figure [L9l 
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incident hard 
1 X-ray photon 




ejected K electron 



K alpha photon ]* 
E = 6.4keV 




iron atom before photo- ionization 



iron atom after piioto-ionization 



Figure 1.4: A schematic sketch of the photoionisation of a neutral iron atom followed 
by production of a Ka photon. The three atomic shells are shown, from inside: K, L, 
M. Two electrons of the N shell are grouped with the M shell electrons in the picture for 
simplicity. 



The Compton hump is the result of the increased importance of Compton scattering 
compared to the bound-free absorption. The cr oss-section of the bound-free ab sorption 
CTbf decreases namely with increasing energy (|Morrison & McCammonI Il983h . except 
for absorption edges, whi le Compton cross-section ctc is significant up to £ > 50keV 
(ILightman & Whitelll988h. This effect m akes th e overall spectral harden ing and forming 
of the Compton hump JUghtman & White.19881 : lOuilbert &Reeslll988h . 

In addition to the reflection continuum, intrinsically narrow features, absorption edges 
and fluorescent emission lines, are other significant imprints of the reflection from an ac- 
cretion disc. These spectral features have a great potential in investigation of the prop- 
erties and localisation of the original source of radiation. They are defined with specific 
energies and so, detecting of energy shifts and broadening of the features is an important 
observational tool which allows us to investigate the innermost region of a black hole 
accretion disc. 

A schematic picture of photoionisation followed by production of an X-ray photon is 
shown for the case of a neutral iron atom in Figure [L4l The incident radiation must have 
high energy (hard X-ray) to be able to kick ofi" an inner electron. The minimal energy 
of the incident photon is equal to the binding energy of the electron in the atom, and it 
corresponds to the absorption edge in the reprocessed spectrum. The originally neutral 
atom is photoionised and the electron vacancy is immediately filled up by an electron 
dropping down from a higher atomic level, L shell or M shell. This is accompanied by a 
release of the energy equivalent to the difference between the energy levels either in the 
form of radiation (fluorescence), or by ejecting of an outer electron (Auger effect), q 

The fluorescent emission line is named according to the shell from which the electron 
was ejected by the photoionisation, K line if the ejected electron is from the K shell, L 



^The Auger effect may be interpreted as a further, inner, photoionisation. A photon released by dropping 
down of an electron from higher to lower level (from L to K shell) is absorbed by an outer electron. The 
photon has enough energy to eject the electron from the atom. 
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Figure 1.5: Left: Disc reflection spectrum originated by irradiation of a neutral slab of 
gas by an incident power law radiation (dashed line). Individua l fluorescent lines are ap- 
parent the iron line being most prominent. Figure adopted from iRey noldsl (| 1 997h . Right: 
Reflection spectra for diff"er ent ionisation st ates (from bottom to top ^ = 30 - 10000). 
Figure adopted from Ross. Fabian & YoungI (|l999,) . 



line if it is from the L shell etc. Further division follows the origin of an electron filling 
the vacancy in the inner shell. If it comes from a neighbouring shell a Greek letter a 
is added to the name, and gradually with higher levels (J3, y,...). The Ka line, shown 
in Figure ll.4[ means that a K electron was ejected by the photoionisation, and an L- 
shell electron jumped into the vacancy. This is a more probable transition than filling the 
vacancy by an M-shell electron. The probability ratio of Ka to K/3 appearance is well 
defined quantity and it is about 7.4 : 1. Similarly, if the vacancy by photoionisation is 
produced in the L shell, ha photon is radiated if an M shell electron fills the vacancy. 

The probability of the fluorescence occurrence, i.e. probability that an X-ray photon 
is radiated from the atom after photoionisation, is characterised by the fluorescence yield 
CO. This quantity depen ds strongly on the atomic number, approximately as a» « Z"* (e.g. 
Bambynek et Q/.lll972r) . Taking into account cosmic abundances as well, the iron fluo- 
rescent lines are expected to be particularly strong. Moreover, the iron Ka and ¥ifi lines 
occur at the energies (6.4 keV and 7.06 keV, respectively) where both, the thermal radia- 
tion of an accretion disc and the reflection continuum (represented mainly by the Compton 
hump) are minimal, and where only a few absorption lines of lighter elements occur. This 
all makes the iron K lines to be relatively easily observable spectral features (see left panel 
of Figure 1 1.51 ). Properties of the iro n Ka line in the X - ray ill u minated cold accretion disc 
pict ure were extensiv ely studied by lGeorge & FabianI (|l99l|) : lMatt. Perola & Pirol (|l991|) 
and lMattg^oZl (Il992h . 

The energy and the intensity of the fluorescent line depends significantly on the ioni- 
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sation state (IRoss & Fabianll 19931 : Matt. Fabian & Rossll 19931 : iRoss & Fabianll2005l) . The 



ionisation parameter may be defined as: 

m = —^, (1-52) 

where F^ is the flux received per unit area of the disc at a radius r, and n{r) is the co- 
moving electron number density. 

The intrinsic energy of the fluorescent line monotonically increases with higher ioni- 
sation state because the binding energies of the inner shells increase with ionisation. In the 
case of iron, the energy of the Ka is very close to 6.4 keV up to Fe XVII, then the energy 
value increases up to 6.7 keV for FeXXV (helium-like iron atom), and finally reaches 
the value E = 6.97 keV for Fe XXVI (hydrogen-like iron atom) when the last electron 
is ejected by photoionisation and another electron is caught up by the fully ionised atom 
and drops down to the innermost K shell. The energy of the K photon depends also on the 
sub-level of the dropping down electron, which makes the K emission line to be a dou- 
blet with energies ^Kai = 6.404 keV and E^ai = 6.391 keV This difference is, however, 
very small and beyond the resolution abilities of the detectors on-board the current X-ray 
satellites. 

The intensity of the fluorescent line is maximal for ^ ^ 1000 - 3500. This is due to 
two eff'ects. First, the fluorescent yield is higher for ionised iron atoms than for neutral 
ones. Second, the incident radiation which is strong enough to ionise iron atoms is not so 
much photoabsorbed by dissociated lighter elements. The fluorescence does not occur if 
all the matter is too highly ionised, so that all electrons are unbound. The iron fluorescent 
line cannot be produced if ^ > 5000 erg cm s"^ Reflection spectra for different ionisation 
are shown in the right panel of Figure [T75l 



1.4 Observational evidence of accreting black holes 
1.4.1 Stellar-mass black hole binaries 

The first black hole candidates to be identified were celestial bodies of small size and mass 
only a few solar masses, in close orbits to ordinary companion stars emitting intense and 
rapidly flickering X-rays. This emission is attributed to the radiation of inward- spiralling 
matter in the form of a ccretion disc (see Section[l.3.2l) . The fir st established black hole bi- 



nary was Cygnus X-1 (|Webster & Murdinlll972 : Bolton 1 19721) . This object is persistently 



bright in X-ray s because its companion is a blue super-giant of spectral type 09.7Iab 



(|Walbornlll973l) which fuels the accretion onto black hole by large amount. These types 
of objects are classified as High-Mass X-ray Binaries (HMXBs) and next to the Cygnus 
X-1, two black hole binaries in the Large Magellanic Cloud, LMC X-1 and LMC X-3, 
belong to this category. 

More frequently observed black hole binaries are, however, Low-Mass X-ray Binaries 
(LMXBs), or also called X-ray novae, which are transient and change rapidly the spectral 
state according to the accretion rate of the infalling matter. Two of them are especially 
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remarkable, GRS 1915+105 belongs to the brightest X-ray objects on the sky (excluding 
the Sun), and since its eruption in August 1992 it remained very bright for more than one 
decade; and GX 339-4, which undergoes frequent outbursts followed by very faint states. 
In total, about two dozens of black hole binaries are confirmed and two other dozens are 
the candidates. For review abou t X-ray properties of the observed black hole binaries see 



Remillard & McClintocM (l2006l) . 



1.4.2 Active galactic nuclei 

Another class of the astrophysical black holes are super-massive black holes (SMBH) 
which settle in the dynamical centres of galaxies. The immense nuclear activity of some 
galaxies was detected already at 1930s by radiotelescopes, but the spatial resolution was 
poor to measure the position of the radiation source. In the early 1940s, Carl Seyfert 
discovered intense nuclear activity together with the presence of highly ionised and ex 



treme ly broad (up to 8500 km s ^) optical emission lines in a sample of galaxies (|Seyfert 



19431) . Consequently, a new subclass of galaxies, Seyfert galaxies, was established. 



Another windows to the Universe, especially X-ray and infrared, revealed that some 
galactic nuclei are enormously bright compared to the rest of galaxies in the whole spec- 
tral energy range. These objects are commonly denoted as Active Galactic Nuclei (AGN). 
Some AGNs are faint radio sources like M31, but other nuclei, like the one in a quasar 
3C 273, belong to radio-loud AGNs, which are characterised by coUimated jets of ener- 
getic particles spanning millions of light years into the space. 

Such an activity of galactic nuclei is attributed to the accretion on the super-massive 
black hole with the mass of several millions to billions of solar masses, Mbh « 10^ - 



10 Mq, (|Reeslll984|) . The precise value of the central black hole mass is measured from 
the velocity dispersion of the stars orbiting near to the centre (see Section [T.5. II for more 
details). 

Extensive mapping of properties of the individual AGNs across the whole electromag- 
netic spectrum lead to the origin of several different empirical AGN subclasses, which 
gradually evolved into a realisation that a unification into a single family of intrinsi 



cally similar active gala xies may be possible (lAntonucci & Milleiill985HAntonuccilll993 



Urry & Padovanilll995|) . The different appearance is mainly due to the orientation effect, 
see Figure 11.61 The activity in the radio spectral energy range, so called "radio loud- 
ness", is the only large distinction, probably connected with the presence or absence of 
the relativistic jet. 

Particularly interesting subclass of the AGNs are Seyfert galaxies which belong to the 
low-luminosity radio-quiet AGNs. Its significant role in the unification scheme is mainly 
due to the relatively frequent occurrence at low redshift (close to our Galaxy), which en- 
ables a spatial resolution unachievable for the distant quasars. The Seyfert galaxies are 
divided into two main groups, Seyfert 1 and Seyfert 2, according to the presence or ab- 
sence of broad lines, especially Balmer lines of hydrogen. According to the unification 
scenario, as seen in the Figure II. 6[ the broad line region of Seyfert 2s is obscured by a 
torus surrounding the central region as the nucleus is seen under a high value of the incli- 
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Figure 1.6; Standa rd picture of the AGN unified scenario (jAntonuccil 119931 : 
Urry & Padovanilll995h . All the observed active galaxies have the same intrinsic structure 
but are viewed under a different inclination. (QSO = quasar, Sey = Seyfert galaxy, BLRG 
= Broad Line Region Galaxy, NLRG = Narrow Line Region Galaxy, OVV = Optically 
Violent Variables, BL Lac = BL Lac Objects) 



nationHI The torus does not enter the line-of-sight between us and the central broad line 
region when we see the nucleus under a low value of the inclination. The typical value 
of the inclination angle of Seyfert Is is around 30 degrees. More sensitive de tectors con- 
struct ed in 1980s enabled fainter division of Seyferts into spectral subclasses (lOsterbrock 
19891) . The spectral type of a Seyfert galaxy may be expressed as (IBlandford et a/.ll 19901) : 



Spectral type = 1 + 



narrow-line flux 
total flux in lines 



nO.4 



(1.53) 



In spite of the unquestionable success of the standard unification scheme by lAntonucci 



(| 19931) . some observational facts are not explained within this standard AGN picture. Es- 
pecially, broad absorption lines {^ 0.1c) in about 10% of quasars and highly ionised out- 
flows (v « 1000km s"') in narrow absorption lines in a h alf of S eyfert galaxies belong to 
the issues which are beyond the standard picture. Hence, l ElvisI (2000) proposed a differ- 
ent structure for quasars which can be also applied to other AGNs, see Figure 11.71 The 
warm highly ionised medium, WHIM, is in the form of an outflow from the disc and has a 



''Standard convention is that the incHnation angle is zero when we see the disc/torus along the symmetric 
axis. We say that we see the disc "face-on". Oppositely, if the observer is in the disc plane, the incUnation 
is 90 degrees and the disc is "edge-on". 
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Figure 1.7: Proposed structure of quasar by lElvisl(|2000l) . The figure is divided into four 
quadrants which illustrate the following (clockwise from top left): the opening angles of 
the structure, the spectroscopic appearance to a distant observer at various angles (BAL 
= broad absorption lines, NAL = narrow absorption lines), the outflow velocities along 
diff"erent lines-of-sight, and the typical column densities of the absorber. 



conical configuration. The opening angle of the cone is about 70-90 degrees and its width 
is at maximum about 10-20 degrees. If we translate it to the probability of observation of 
a quasar through the conical sheath then broad absorption lines of the WHIM should be 
detected in about 10-20% quasars. This corresponds to the proportional representation 
of the so called BAL quasars in the observations. The spectra are characterised with high 
values of column density of the ionised absorber Nu « 10^"^ cm"^. 



1.4.3 Black hole binaries versus active galaxies: similarities and dif- 
ferences 

Both, black hole binaries and active galaxies, are powered by black hole accretion. The 
masses of black holes in binaries and galactic nuclei are very diff"erent (see previous Sec- 
tions). The typical length and time scale differ accordingly (eq. ll.ll) . However, the phys- 
ical mechanism of the energy balance seems to be the same, and so, the AGNs can be 
interpreted as scaled-up Galac tic black holes, or vice- vers a. Galactic black holes are of- 
ten called as "microquasars" (JMirabel & Rodrigueall998 b. The unified picture of both 
types of objects is illustrated in Figure [L8] where another object, coUapsar, is added to the 
same family of accretion powered objects, as well. The coUapsar occurs when a massive 
star undergoes the final gravitational collapse, and it appears only in a short flash. If the jet 
is aligned with the line-of-sight the objects appear as blazars, microblazars, and gamma 
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Figure 1.8: The same physical mechanism can be responsible for three different types of 
objects: microq uasar (left), quasar (middle), and coUapsar (right). Figure is adopted from 
MirabelldlOOeh . 
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Figure 1.9: Model of intrinsic X-ray spectrum of a Seyfert 1 AGN (left) and a stellar-mass 
black hole binary (right). Relativistic smearing is not included in the figure. The figure 
illustrates similarity of the AGN and BHB X-ray spectrum. Hard X-rays are dominated 
by the power-law component (red curve) which partially reflects on the disc (green curve). 
The reflection component is dominated by the Compton hump at £ a: 20 - 40 keV and the 
iron fluorescent line. In the figure, the accretion disc of the AGN is assumed to be cold 
while the disc of BHB is ionised, and so, the shape of the reflection component differs 
accordingly. The soft excess (blue curve) observed in several Seyfert galaxies is shown 
here as a black-body component but it also may arise through disc reflection (for more 
details, see Section [T.4. 31 ). Figure is adopted from .Miller. (2007^) . 
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ray bursts, respectively. 

Besides the similar appearance of A GNs and BHBs the temporal beha viour seems to 
be the same, only proportionally scaled (IMcHardy et a/.ll2006l : IUttleyll2008 V The minimal 
time scale for the spectral variability is limited by the black hole horizon (eq. ll.ll) : 

T >'' -^^ n 54^ 

-t variability ^ — - T- ■ U-J'+J 



The typical time scale for stellar-mass black hole binaries is of order of milliseconds, for 
super-massive black holes in AGNs it is of order of minutes to hours. 

The very rapid flux variability is a characteristic feature of both kinds of objects 
and was constrained in several observations. The RXTE (Rossi X-ray Timing Explorer) 
satellite is suitable for studying millisecond variability of BHBs (see e.g. the review by 
van der Klisll2006L and references therein). The variability of A GNs was also reported in 
several works (see e.g. Ilwasawa et a/.lll999l: Turner g? a/.lll999l: iMiniutti & Fabian 2004 : 



Vaughan & FabianI 120041; luttlev & M cHardvl bOOSi : iTumer et alhood IPonti et a/.lbood : 
Gierlihski et a/.l2008l:lde Marco et a/.l2009i) although more sensitive instruments or longer 
exposure times would be eligible to increase the accuracy. 

The time-averaged spectrum of a black hole binary evolves as the spectral states 
change, but in general, it may be interpreted as a compositio n of power-la\y radia tion, 
thermal multi-colour black-body radiation of the accretion disc (|Mitsuda et q/.II1984|) . and 
reflection radiation from the disc, all absorbed by a local (possibly ionised) matter and 
the interstellar Galactic gas. The model spectrum of BHB is shown in the right panel of 
Figure [L9l The model spectrum of an AGN is very similar (left panel of Figure [L9l ), but 
the thermal component is dominated in the UV energy range, and, thus, it is obscured by 
the interstellar gas, which is opaque to UV photons. Nevertheless, the high energy tail of 
a black-body spectrum can reach soft X-rays if the temperature of the disc is high enough. 
Indeed, a "soft excess" is observed in several Seyfert galaxies which may be interpreted 
by the therm al radiation of the disc (as shown in the left panel of Figure 1 1.91 ). 

However. ICrummy et q/.I(|2006|) and lGierlihski & Done! (|2006l) argued that the thermal 
radiation interpretation is doubtful since the temperature would be the same for a big sam- 
ple of AGN with distinct masses, and therefore, they suggested alternative explanations. 
The "soft exces s" may be due to reflection by the ionised surface of the accretion disc 
(|Crummy et a/.l l2006). For a cold disc, the reflection is dominated by a plenty of emission 
lines in soft X-rays. In this case, the "soft excess" may also occur when the reflection ra- 
diation comes from the innermost regions of an accretion disc and the individual lines are 
blended due to relativistic smearing. Another interpretation of the "s oft excess" is that it 
appea rs due to the partially ionised and Doppler smeared absorption (|Gierlihski & Done 
2006h . 



1.5 Measuring black hole parameters 

Astrophysical black holes are characterised only by their mass and angular momentum. 
Although electrically charged black holes also represent the solution of Einstein's equa- 
tions, their existence in the Universe is very unlikely since the electromagnetic repulsion is 
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about 42 orders of magnitude stronger than the gravitational attraction. Thus, any charged 
black hole will be immediately neutralised by an accreting matter. 



1.5.1 Mass 

Masses of astrophysical black holes are known with relatively good precision and cus- 
tomarily categorised in three groups of stellar-mass black holes (< 30Mq), intermediate 
mass b lack holes (10^-10^ Mp), and super-massive black holes (> 10^ Mp) in galactic 
nucle i (lCasaresll2007l : IZi61kowskill2008l : l^tergaard & Osmei]l2009l : [Czemv & Nikolajuk 
20091). The intermediate mass black holes may occur in centres of globular star clusters 

and represent a plausible explanation of so called ultra-luminous X-ray sources (ULXs) 

(lKaaretg?Q/.l200ll : iMiller & Hamiltonl2002l : iMiller & Colbertll2004l : Istrohmaver & Mushotzkv 
20091). but their exi stence has not yet been unambiguously accepted (e.g. iKingI l2004l : 
Soria & Ghoshll2009 ). 

The mass of the black hole Mbh in a binary system can be determined from the knowl- 
edge of the inclination i of the system and the mass ratio between two components of the 
binary (the mass of the companion star can be estimated from the spectral type of the 
star). The mass function is given by: 



/(M) = 



P^ 
2nG 



Mbh sin i 

(1 + m' 



(1.55) 



where P is the orbital period, Kq = v sin / is the maximal line-of-sight Doppler velocity 
of the companion star, and Mq is the mass of the companion star, which is usually lower 
than the black hole mass (if it is not a super-giant as in the case of the Cyg X-1 system). 
The mass of black holes in galactic nuclei can be derived from the dynamical prop- 
erties of the stellar neighbourhood. The first velocity moment of the coUisionless Boltz- 
mann equation, the Jeans equation, gives t he mass as the function of radius (see e.g. 
Binnev & Tremainell2008l : ICamenzindll2007 ): 



2 2 

rv rat 



Jlnp, Jlncr^ / cr- 
d\nr d\nr \ al 






(1.56) 



where v is the rotational velocity, a] are velocity dispersions and p* is the density of 
stars. All these quantities are measurable, but several problems exist. First, we observe 
these quantities projected on the sky, and de-projection is rather complicated process. 
Second, the most galaxies are not spherically symmetric. The high quality of the data 
are required. Observations performed by the Hubble Space Telescope can serve to extract 
the line-of-sight velocity dispersions (LOSVD). Then, an orbit based approach, known as 
Schwarzschild's method is applied to de-project the quantities. Finally, maximum entropy 
models are used to account for axisymmetry instead of the spherical symmetry, and the 
potential is derived from the profile of the surface brightness. The velocity maps can 
be achieved with higher precision from rn icrowave maser emission of water molecules 
(iMivoshi et a/.lll995l : lOreenhill et q/.II2002|) . 
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In order to estimate the black hole mass one can use the fact that the stars around more 
massive nucleus orbit with higher velocities. The empirical relation between the black 
hole mass and velocity dispersion of the nearby stars, "A/ - cr" relati on, was found on the 
samp l e of galaxies with rel atively known values of the central mass CFerrarese & Merritt 
2000l : lGebhardt et a/.ll200ol) : 



log(MBH/Mo) = or +^6 log ((T, /200 km s"^) 



(1.57) 



E arlier, the correlation between t he central mass and the bulge luminosity was discov- 
ered (IKormendy & Richstonelll995l : iMagorr ian et al.\\l99S\) . This relation is often called 
"Magorrian" relation. It has larger scatter than the "M - cr" relation, and hence it is used 
less often at present. However, both correlations have important implications for theo- 
ries of the galaxy and bulge formation and the interactions of the stars with the central 
super-massive black hole. 

In fact, the value of the mass determined using such relations is rather an order-of- 
magnitude estimation than the precise measurement due to the relatively large scatter of 
the correlations. Our own Galaxy is an exception in the mass determination of black 
holes in galactic nuclei because it is so close (^ 25 000 ly) that the individual stars or 



biting around the central black h ole are observed (Eckart & Genzell 1 19961 : iGhez et al. 



1998L |2005|; iGillessen et al.\ 12009) . The latter one presents the black hole mass to be 



^BH = (4.3 + 0.2stat ± 0.3sys) X 10^ Mo, where the uncertainty is given by the statistical 
and systematical error, in which the uncertainty of the distance is included. The suitable 
view to measure dynamical properties of the closest stars to the Galactic centre is due to 
the most sensitive infrared detectors at the largest telescopes (VLT, Keek's telescopes). 

Th ough the infrared o bservations reveal some accretion activity in the form of the 
flares (IGenzel et al.\ l2003h . there is no observational evidence for a standard Shakura- 
Sunyaev accretion disc in the Galactic centre. The advection dominated accretion flow 
is likely to explain the X-ray sp ectrum and also the luminosity of the Galactic centre 
(|Narayan. Yi & Mahadevanlll995|) . 



1.5.2 Angular momentum (spin) 

The mass of a black hole is relatively easy to measure because the attractive gravitational 
force reaches to the large distance and aff"ects therefore orbital movement of the com- 
panion star in a binary system or the surrounding stars and gas in the central region of 
a galaxy, respectively. More challenging is to measure the value of the spin of a black 
hole because the spin causes curvature of the space-time in a detectable level only within 
a few gravitational radii around the black hole. Despite this small outreach, the black 
hole spin plays an important role in the black hole energetics, especially, it is assumed 
to be responsible in g enerating and up-keeping of the powerful relativistic jets (" Penrose 
19691 : iBlandford & Znajek.1977.) . The information about the spin value on a statistically 
significant sample of black holes is important in the understanding of the formation and 
the growth of black holes. It can significantly help to answer the question if the observed 
spin value is natal or if the black hole ro tation is accelerated via the accretion (see e.g. 
King & Kolbll 19991: IVolonteri et Q/.ll2005h . 
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There are several observational methods which provide a good opportunity to explore 
the innermost region of an accretion disc, and thus, to constrain the spin value. The 
summary of them is listed and briefly discussed: 



Continuum fitting 

In the thin Shakura-Sunyaev disc, the temperature decreases with the distance as given 
by eq. (I1.46I ). When the innermost edge of the disc corresponds to the last stable circular 
orbit, which should be the case in at least one accretion state (high/soft), its position can be 
determined from the spectral fitting of the thermal component. This is the basic concept 
of the method which is called X-r ay continuum fitting method and which was first carried 
out by IZhang. Cui & ChenI (|l997|) to measure black hole spin. 

The temperature does not depend only on the inner edge of the disc, but also on the 
accretion rate. The thermal spectrum is further distorted by the rotational and gravitational 
frequency shif t, and therefore, the fully relativistic model is acquired for this purpose 
(JLi et a/.ll2005L kerrbb model). The spectr al hardening facto r f^oi = Tcoi/Tss as a function 
of the Eddington scaled disc luminosity (|Davis et al.\\2005\) and a switch parameter for 
zero/nonzero torque condition at the inner edge are included in the kerrbb model. The 
spectral hardening factor (or also call ed colour correction) fr ni plays an important role 
especially in higher accretion rates (Shimura & Takaharalll995|) . 



Iron Ka line profile 

The iron Ka line profile originates by the reflection of hard X-rays on the accretion disc 
surface. The combination of three efl'ects makes the iron Ka line easily detectable in X- 
ray spectra of black hole binaries as well as active galaxies. First, the fluorescent yield 
is higher for heavier elements (^ Z^ , see Section 11.3.31) . Second, the relative cosmic 
abundances of iron are high compared to other heavy elements. And last, the energy of 
iron Ka line occurs in the spectral range of a simple continuum (see Figure [L9l) . 

The profile of the intrinsically narrow emission line is distorted by thermal motion, 
Compton scattering, Doppler broadening, and relativistic efl"ects on the radiation includ- 
ing relativistic Doppler shift, lensing, gravitational redshift, and time delay (if the radia- 
tion source is not steady in time). The Doppler broadening due to the rapid orbital motion 
and the relativistic effects are much more significant than the thermal or Compton broad- 
ening if the emitted radiation comes from the inner region of an accretion disc (within a 
few hundreds of r^), maki ng the line profile extremely smeared even for a non-rotating 
Schwarzschild black hole (IFabian et a/.lll989|) . The effects are amplified in the case of a 
rotating black hole since the marginally stable orbit shifts closer to the bl ack hole (se e 
Sec. II. H and Fig. II. 21) . and are maximal for a maximally rotating black hole (lLaodll99l|) . 

The advantage of this method is that the line profile is completely independent of 
black hole mass, and is \yidelv applicable to BHBs as w ell as to AGNs (for reviews see 



Revnolds & Nowakll2003l : lMiUed 



20071 : lGuainazzill2009l) . The line profile depends on the 



source properties - its geometrical position, which is influenced by the black hole spin 
(the innermost edge of a disc), and orientation. This makes this technique a suitable tool 



1 .5. Measuring black hole parameters 



25 



for investigation of the nature of the innermost region of the accretion disc, and also for 
measurement of the inclination angle of the disc. This method is more widely investigated 
and discussed in the following Sections of the Thesis. 



Quasi-periodic oscillations 

X-ray quasi-periodic oscillations (QPOs) are transient phenomena associated with the 
non-thermal states and state transitions. There are two kinds of QPOs - low frequency 
quasi-periodic oscillations (LFQPOs) at roughly 0.1-30 Hz in the power density spectra 
(PDS), and high frequency quasi-periodic oscillations (HFQPOs) at roughly 40-450 Hz 
in PDS. The typical frequency of HFQPO corresponds approximately to the orbital fre- 
quency at the ISCO, and ar e thus relevant for the spin measurement. Although several 
models were suggested (e.g.lStella & Vietrill 19981 : iTitar chuk. Lapidus & M uslimovl 1 9981: 



Abramowicz & Kluzniakll2001l : iRezzoUa et a/.l 12003 : Torok et al.\ l2005l : iHorak & Karas 



20060 . the satisfactory description of QPOs is not well established. 



Variability and reverberation 

Small size of a "hole" in the inner accretion disc implies the variability on the relatively 
short time-scales. Temporal changes in the primary radiation are translated to the re- 
flection radiation with a certain time lag r. The characteristic time lag is approximately 
"T ~ rms/c. A clear evidence for such a lag, a 30 s reverberation lag between direct X-ray 
continuum and FeL emission accompanying the relativistic reflection, was reported by 
Fabian et al.\ (120091) . This measurement was possible thanks to a high overabundance of 
iron in this particular galaxy and a very long exposure time by XMM-Newton satellite. 



Polarimetry 

The last method of the spin measurement is via X-ray polarimet ry. The polarisation 



of X-r ays from accretion discs around black holes was studied by iLightman & Shapiro 
(| 19751 ). The thermal emission is polarised due to Thomson scattering in a disc atmo- 
sphere. I Connors. Stark & PiranI (|l980.) showed that the polarisation features are strongly 

afi'ected by general relativistic effects. Other authors have considered efi'ects of magnetic 

fields on the resulting polarisation (e.g. lAgol. Blaes & Ionescu-Zanettilll998l : ISilant'ev & Gnedin 
2008h. Especiallv. ro t ation of the polarisation ang le is a sensitive quantity. Recent studies 
(|Dovciak et a/.ll2008l : ISchnittman & KrolikI 120091) illustrate how the polarisation features 
depend on the spin value in the thermal state of an accreting black hole. Although the 
models capab le to compute Stokes parameters of a polarised accretion disc spectrum are 
ready for use (IDovciak. Karas & Yaqoobll2004|) . any X-ray polarimeter useful for this type 
of measurement has not yet been launched. This method is hence promising in the future. 



Chapter 2 



Relativistic lines 



2.1 Line profiles in strong gravity regime 

In classical astronomy, broad emission lines are useful diagnostic tool for measuring 
of the temperature (thermal broadening) or velocity dispersion (Doppler broadening) of 
an observed system. However, X-ray spectroscopy of some X-ray binaries and active 
galaxies has revealed iron fluorescent lines so broad and asymmetric that their profiles 
cannot be explained in terms of classical physics and instead, a complex fully general 



relativistic approach needs to be taken into accoun t (Fabian etal., ,19891: iTanaka et ah 
19951 : [Reynolds & Nowakll2003l : iNandra et a/.ll2007l: iMiller 2007k The broad iron lines 



are supposed to originate in close neighbourhood of a black hole where the strong gravita- 
tional redshift occurs and the orbital velocities reach a considerable fraction of the speed 
of light (see Fig. O)- 



2.1.1 Frequency shift 

The frequency shift, g, is defined as the ratio of the observed frequency Vobs to the intrinsic 
emitted frequency Vem: 

Vobs=^Vem- (2.1) 

In high-energy astronomy, energies are usually used instead of frequencies: 

h c 
E = hv=—, (2.2) 

A 

where h = 6.626 x 10"^"* J s is the Planck constant, and A is the wavelength. Typically 
used units are kiloelectronvolts (1 keV= 1.602x 10"^^ J) for energies, or Angstroms (lA = 
10"^" m) for wavelengths. The conversion relation is: 

12.4 
^[keV] = -— . (2.3) 

4a] 
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In the case of a rotating accretion disc, the g-factor is a function of the position on the 
disc and the emission angle g = g{R, ip, Qg). In Newtonian approach: 



5 class 



1- 



Vi 



los 



(2.4) 



where Vios is the line-of-sight velocity. For Keplerian orbital velocity, Vios 
In special relativistic approach: 



1 



,?STR - 



(i-^r 



f^cos v^jsin^. 



(2.5) 



r 



where y = 



V^ 



is the Lorentz factor. The observed frequency is shifted even when the 



line-of-sight velocity is zero (Vios = 0, "face-on" disc). Since ^str = 7~' < 1 in this case, 
the frequency is shifted to lower energy. This effect is called as the transverse Doppler 
shift. 

In general relativity, g-factor may be expressed within the approximatio n of geometri 



cal op tics, in terms of four-momentum of photons p^ and four-velocities u^ (|Cunningham 



19751 



f 



8 






(2.6) 



Using u'^obs = [-1, 0, 0, 0] (four- velocity of the observer): 

-(PtXbs 



(p,um\ 



(2.7) 



The frequency shift g and the emission angle 9e may be expressed using constants of 

motion as: 

C „ ^V^ 



S 



-3/2^' 



6e = arccos 



(2.8) 



where !B = I + ar ^^^ , C = I - 3r ^ + la r ^^^; $ and ri are constan t s of rnotion connected 
with symmetries of the Kerr space-time fNovikov & Thomelll973 : Karas 2006)^ 



T he computation of g-factor c an be provided numerically ( e.g. iDov ciak. K aras & Yaqoob 



20041 : iBeckwith & Done.2004; Brenneman & Reynoldsll2006h . or some feasible a pproxi 



Beloborodov 



mations were developed for Schwarzschild metric to perform fast calculations (e.g 
2002l : IPechacek et a/.lllooi) . 

The contours of constant frequency shift g around a non-rotating, a = 0, and a max- 
imally rotating, a = 1, black hole are shown for two different values of the inclination 



angle in Figure 12 .1 Kan atlas of ^-factor for different values of spin and inclination may be 
found in lDovciakll2004r) . The value of the frequency shift g results from two effects, grav- 
itational redshift and Doppler shift. The effect of gravitational redshift quickly increases 
when approaching the black hole, and it is infinite at the black hole horizon {g = 0). Near 
the black hole, the gravitational redshift clearly dominates over the Doppler shift. 
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Figure 2.1: Contours of redshift factor g{r,(p) near a non-rotating Schwarzschild black 
hole, a = 0, (top) and a maximally rotating black hole, a = 1, (bottom), depicted in the 
equatorial plane ix,y). The black hole and the accretion disc rotate counter clock-wise. 
A distant observer is located towards the top of the figure. The inner region is shown up 
to r = 20 gravitational radii from the black hole. The black hole is denoted by a dark 
filled circle around the centre, for Schwarzschild black hole the circle around represents 
the marginally stable orbit. Two cases of different observer inclinations are shown. Left: 
6o = 30 deg. Right: do = 70 deg. The colour bar encodes the range acquired by g{r, (f), 
where g > I corresponds to blueshift (approaching side of the disc), while g < 1 is for 
redshift. 
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The magnitude of the Doppler shift depends on the line-of-sight velocities. Its effect 
is also maximal at the closest orbit to the black hole where the orbital speed is maxi- 
mal (see Figure 11.21) . Therefore, the maximal redshift always occurs at the innermost 
orbit. However, the position of the maximal blueshift is farther from the black hole due 
to the gravitational redshift, and its exact position depends on the inclination angle. For 
lower inclinations, the Doppler effect is weaker and hence, it overcomes the gravitational 
redshift at a farther radius. Sufficiently far from the black hole, two distinct regions of 
blueshift and redshift are separated according to the rotation and the direction of the line- 
of-sight velocity to the observer. The level of the frequency shift decreases with the radius 
as the orbital velocity decreases as well. 

2.1.2 Intensity 

The specific intensity / is defined as energy emitted at some given frequency v into the 
element of solid angle Q in unit of time t. 

■ ""^ (2.9) 



dMvdQ' 



where A^ is the total number of photons. 

The value of the intensity detected by an observer depends on the frequency shift of the 
emitted radiation. This relation may be derived from the Liouville's theorem which states 
that the phase-space volume Y is invariant to the canonical transformations representing 
the time evolution of the system: 



X 



d"pd''q = C, (2.10) 



where n is the dimension, q expresses coordinates and p the conjugated momenta. In our 
application, the Liouville's theorem states that the phase-space density: 

N 
n = -, (2.11) 

is a constant of the Lorentz transformation. The element of the phase-space volume is 
given by: 

Y = d^pd^x = 4np^dpdQcSdt. (2.12) 

Thus, 

"~ 47rp2dpdQcSdt' ^ ■ -* 

Substituting from (eq. l2.13|) into (eq. 12.91) and using p = '-f'. 

^^4nh^^ (2.14) 



and 



C3 



— = const. (2.15) 
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Thus, using eq. (12.11 ) 

-'obs ^obs 



/em V, 



3 



(2.16) 



em 



This fact implies that the radiation from a matter approaching to the observer (blueshift) 
is boosted while the radiation from a matter receding from the observer (redshift) is dimin- 
ished. The level of the intensity amplification (resp. diminution) depends on the values of 
the line-of-sight velocity. 

2.1.3 Line profiles from accretion discs 

For line emission from an accretion disc, we can assume an axial symmetry and separa- 
bility of the radial and angular emissivities in some cases. Then, the line emission from 
the disc can be written in the form of product: 

I^m{re,fie, ^e) = '^(''e) )"e M(jU^, E^) 6{E^ - Eq), (2.17) 

where r^. is the disc radius (distance from the centre), fi^ = cos ^e is the cosine of the 
emission angle measured from the disc normal direction to the equatorial plane, in the 
disc co-moving frame, i.e. in the local Keplerian frame orbiting with the angular velocity 
Q.^(r). Likewise, the intrinsic energy E^ is measured with respect to the local frame. 
Eq = Eq{^) is the intrinsic energy of the fluorescent line depending on the ionisation state. 
The radial part is being approximated by a power law, 

nre) = r:' (<? = const), (2.18) 

or by a broken power law: 

( n{r,) = r-J' (re < r^) 

^(re) = ^, (2.19) 

[ nr,) = r,""- (re > r^) 

The q parameter is typically q>2 which means that the intrinsic intensity of the radiation 
decreases with the distance. The standard value is ^ = 3, larger v alues may occur only 



under certain conditions in the innermost parts of the disc (see e.g. iNiedzwiecki & Zycki 



2008) 



The angular emissivity law, M{jj.s,Es), defines the distribution of the intrinsic in- 
tensity outgoing from each radius r^ of the disc surface with respect to the perpendic- 
ular direction. The limb darkeni ng law in the form A1(/ie,£'e) = A1(yUe) = 1 + 2.06 /ig 



(|Chandrasekhaiill960l : lLaorlll991|) is most frequently used. However, the choice is some- 
what arbitrary in the sense that the physical assumptions behind this law are not satisfied 
at every radius over the entire surface of the accretion disc. This aspect is studied in detail 
in Section |3] of the Thesis. 

The observed radiation flux from an accretion disc is obtained by integrating the in- 
trinsic emission over the entire disc surface, from the inner edge (r = r^^ to the outer 
edge (r = rout), weighted by the transfer function T{r^,(ps,6o,a) determining the impact 
of relativistic energy change (Doppler and gravitational) as well as the lensing effect for 
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Figure 2.2: Schematic sketch of the transformation from the Boyer-Lindquist coordinates 
of the disc r,(p to the detector coordinates a, /3. Credit by M. Dovciak. 



a distant observer directed along the inclinati on angle 6„ (see ICunninghamll 19751 : lAsaoka 
19891 : iKaras. Vokrouhlickv & Polnarevll 19921 : Earas 2006') : 



-f'obs (9, 



o, CI, iIiq) — I 



r(re, (fis, 9o, a) /(re,//e, £e) d^ r^ dre, 



(2.20) 



where the index 'e' denotes quantities related to the disc and 'o' observed quantities. The 
integration is carried out over all possible values of the frequency shift g and the whole 
surface of an accretion disc. 



Another approach is employed in iDovciak. Karas & YaqoobI (|2004l) where the inte- 
gration is carried out over the coordinates and ^-factor is computed at each place. There 
are two possibilities how the observed flux is expressed: 



(2.21) 



1 . in detector coordinates a, /3: 

Fobs (^o, a, Eo) = I Gi(re, (fe, 6o, a) I{rs,n<„ E^) da djS, 

2. in Boyer-Lindquist coordinates of the disc r, ip: 

I G2(re, i^e, ^o, a) I{rs,iUe, £e) r^ dre dv?e, (2.22) 

Jr,„ 



where Gi{r^, (p^, 6o, a) and G2(re, (fe, 6o, a) are another transfer functions. If / is the specific 
intensity given by eq. (12.91 ) then from eq. (|2.16l) : Gi = g^. The relation for G2 may 
be found from the transformation relation between the a, /3 coordinates and the Boyer- 
Lindquist coordinates of the disc r,(p (see Figure [Z2l) : 
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Figure 2.3: Comparison of iron line models produced with the kyrline model for dif- 
ferent set of parameters. If not further specified the model parameters are a = 0.9982, 
q = 3, Tin = Vs(a), Vit = 400, Oo = 30 deg, Eq = 6.4 keV Left: Three curves correspond 
to the spin a = (black), a = 0.7 (red), and a = 0.998 (blue), respectively. Right: Three 
curves correspond to the inclination 6 = 30 deg (black), 6 = 60 deg (red), and 6 = S5 deg 
(blue), respectively. 
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(2.23) 



where So is the observed element surface (i.e. in the detector frame), Sioc is the local 



element surface of the disc, and £ = ^r- is the lensing factor defined as the ratio of the 

loc 

cross-section of the flux tube at the detector to the cross-section of the same flux tube at 
the disc. Using eqs. (IXIB - dXISl) : G2 = g^nJ- 

As a summary, the resulting line profile from an axisymmetric accretion disc is shaped 
by these parameters (see eqs. l2.17l and l2.22l) : 

1 . inclination angle 9o 

2. inner disc radius rin 

3. outer disc radius rout 

4. radial dependence parameter q 

5. angular emissivity law A1 

6. spin a 

Figure 123) shows t he theoretical line profiles for difl 'erent values of the parameters. 
The KYRLINE model (JDovciak. Karas & YaqoobI |2004|) in the XSPEC fitting package 
(|Amaudlll996l) was used to produce the lines. The inner disc was fixed to the radius 



of the marginally stable orbit r^^ = Vs(6[) which is linked with the spin value by eq. (11.21) . 



2.2. Observational evidence of relativistic iron lines 
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The line profile has a more extended red wing for higher values of the spin because the in- 
ner edge of the disc shifts to lower radii where more extreme gravitational redshift occurs. 
The line profile is broader for higher values of the inclination angle because the line-of- 
sight velocities are higher accordingly implying a stronger Doppler effect. The stronger 
Doppler effect overreaches the gravitational redshift in the region where the matter moves 
towards the observer. The position of the maximal blueshift is closer to the black hole 
proportionally to the strength of the Doppler shift and thus, to the inclination angle (see 
Figure [2?T]) . The peak of the line is shifted to higher energies for higher inclinations. All 
the lines shown in the figure are normalised to the same flux. 



2.2 Observational evidence of relativistic iron lines 



Broad iron lines have been detected in several X-ray spectra of black hole binaries as 
well as active galaxies. The first detection of a moderately broad iron li ne is attributed 
to the EXOSAT observations of Cygnus X -1 (JBarr. White & Pagel ll985b. The redshift 
and broadening of the line was explained by Fabian et all (|1989|) as a result of relativistic 
smearing. They constructed diskline modeUfor the relativistic line around a non-rotating 
black hole. A model of the relativistic line f or a maxirnally rotating Kerr black hole 
was developed independently by lKojimal(| 19911) and lLaon (1199 ih . A strongly asymmetric 
and redshifted line profile was predicted for radiation coming from an inner accretion 
disc around a black hole. This indicated that the interesting insights on the geometrical 
properties of accretion discs in the closest vicinity of a black hole may be derived from 
the shape of the line. 

The first X-ray spectrometer capable for resolving the shape of the line was SIS detec- 
tor (Solid-state Imaging Spectrometer) on-board the ASCA satellite (Advanced Satellite 
for Cosmology and Astrophysics). The first resolv ed broad iron line w as detected in 
the spectrum of a Seyfert galaxy MCG -6-30-15 by iTanaka et al.\ (|1995[) . This source 
exhibited an extremely broad iron line in X-ray spectra of all succeeding satellites (see 
Section 1431) . The one of the most suitable and still operating missions for investigation 
of broad iron lines is the XMM-Newton satellite thanks to its large effective collecting 
area in the iron line band and also a few keV above (see Section 14.1.11) . Recently, the 
Suzaku satellite (launched on 10 July 2005) exceeds XMM-Newton with its broad-band 
coverage (0.4-600 keV) which allows better constraining of the continuum (both, primary 
and reflection). 

Broad iron lines shown in Figure [23] for the case of a black hole binary GX 339-4, and 
in Figure [231 for the case of a Seyfert 1 galaxy MCG -6-30-15 belong to the best exam- 
ples of the extremely relativistically broadened lines. We re-analysed the XMM-Newton 
observations of these sources. For a brief description of both sources and details of our 
re-analysis, see Sections [4.21 and [^31 Other examples of active galax ies with broad iron 



lines may be found in e.g. in lNandra et al.\ (|1997[) : iGuainazzil (|2003[) : Istreblyanska et al. 



'This model is embedded in the XSPEC fitting package and may still be used for spectral analysis of 
non-rotating or slowly rotating black holes or neutron stars. However, be aware that the eff'ect of lensing is 
not included in this model. 
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Energy (keV) 



Figure 2.4: Left: X-ray spectrum of GX 339-4 observed in the high/soft state by XMM- 
Newton/EPN (black), RXTE/PCA (red) and HEXTE (blue) (iMiller gra/.ll2004l) . The 
spectrum is modelled by a multi-colour temperature disk black-body and a power law 
revealing a broad asymmetric iron line. Fe Ka line is sho wn in detail in the b ottom panel. 
Right: X-ray spectrum of GX 339-4 in the low/hard state (JMiller et a/.ll2006h . The XMM- 
Newton data observed during the revolution No. 782 are shown in blue, No. 783 in red 
and the data obtained by RXTE are black. Fe Ka line is shown in detail in the bottom 
panel. 



( 2005 ): 
(2007); 
(2009); 



Reeves et al.\ (120061); iTumer et al.\ (120061) ; Piconcelli et al.\ (120061); iNandra et al. 
Longinotti e t al. ('2008*); 'Markowit z eraU (l2008l . |2009|) ; iBrenneman & Reynolds 



Fabian et ah (^2009) ; S\metaL (i2010l). A statistical d istribution of AGN with 



broad iron lines was studied by iGuainazzi. Bianchi & DovciakI (|2006.) who estimated a 
fraction of 42 + 12 % of well exposed AGN s that exhibit a relativistically broadened 
iron line. Works by iLonginotti et al\ (|2008|) and de la Calle et al. (2009, submitted 
to A&A) represent the continuation of this effort. Different s amples of Seyfert g alax- 
ies observed by the XM M-Newton satellite were studied by Nandra et al\ (120071 ) and 
Brenneman & Reynolds! (|2009) . Both groups concluded that the most X-ray spectra of 
their samples possess a rel ativistically br oadened iro n line. Other ex a mples of black hole 
binaries may be found in loiaz Trigo et al\ (|2007b ; iHiemstra et al\ (120091) ; iMiller et al. 



(|2009l) . Relativistic ally broadened iron line s were also detected in several X-ray spectra 
of neutron stars (see lCackett & Milleiil2010l and references therein). 



2.3. Comparison between the LAOR and KY models 
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Figure 2.5: Relativistic iron line in MCG -6-30-15 observ ed by XMM -Newton (red) and 
Suzaku (black). Figure adopted from lMiniutti et al.\ (|2007h . 



2.3 Comparison between the laor and ky models 

There are several numerical codes for the relativistic disc spectral line, some of them 
were already me ntioned. Th e most widely used model, over almost two decades, has 



been the one by iLaoii (1199 1|) . which includes the effects of a maximally rotating Ken- 
black hole. In other words, the LAOR model sets the dimensionless angular momentum 
a to the ca nonical value of q = 0.9982 - so that it cannot be subject of the data fitting 
procedure. lOovciak. Karas & YaqooQ (|2004[) have relaxed this limitation and allowed a 
to be fitted in the suite of KY models. 



O ther numerical codes have been developed independently by several groups ([Martocchia. Karas & Mi 

2000l:IViergutzlll993l:IZakharov & Repinll2004l:lFuerst & Wul2004l :l Beckwith & Donell2004l: 
Cadez & Calvanil 120051: iBrenneman & Reynoldsl |2006|) using different techniques. The 
last three equi pped their codes with sirn ilar functionality as KY to be used in X-ray spec- 
tra modelling. IBrenneman & Reynoldsl (12006 ) performed useful tests demonstrating that 
KY and their KERRDISK model give compatible results when they are set to equivalent 
parameter values. 

Although the LAOR model does not have the spin value as a variable parameter, it can 
still be used for evaluation of the spin if one identifies the inner edge of the disc with 
the marginally stable orbit (eq. I1.7[ see also Figure [L2l) . The evaluation of the spin from 
extension of the red wing only is, however, not precise because the spin also modulates 
the value of the g factor and thus, it affects the overall shape of the line by itself. This 
systematic error of spin measurements obtained with the LAOR model has not yet been 
constrained. Further, we compare the spin estimation by the LAOR and KYRLINE model 
to constrain this error. 
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Figure 2.6: Comparison of the LAOR (black, solid) and kyrline (red, dashed) model 
for two values of the spin a = 0.9982 (left) and a = 0.7 (right). The other parameters of 
the line are £ = 6.4 keV, q = 3,i = 30°. 



The theoretical line profiles of the LAOR and kyrline models are compared in Fig- 
ure [2i6] for two values of the spin, a = 0.9982 (left panel) and a = 0.7 (right panel). The 
inner radius of the disc is set to coincide with the marginally stable orbit for a given value 
of the spin, i.e. Tin - 1.235 and rin - 3.39, respectively. The angular dependence of 
the emissivity in this example is given by the limb darkening law (/(//e) oc 1 + 2.06 /ie) 
which is the only prescription for the directionality in the LAOR model. The kyrline 
model enables to switch between different emission laws. Further in this Section, we 
used two distinct cases and this name convention: kyrline is equipped with the same 
limb-darkening law as in the LA OR model, an d KYRLINE* uses the limb-brightening law 
in the form IQu^) oc ln(l + ^) (lHaardtlll993l ). 

In the first comparison (left panel of Fig. 12.61) . the difference between the two models 



is noticeable only in the high energy edge of the line profile. This is caused by poorer 
energy resolution of the LAOR model. In the second comparison (right panel of Fig. 12.61) . 
an additional discrepancy appears in the slope of the line profile (look at energies 4.5- 
6 keV). This is the effect of the different frame-dragging taken into account and it is caused 
by the fact that the LAOR model is calculated only for a = 0.9982. 

Relatively to the total line profile, the difference between the LAOR and KYRLINE 
model seems to be rather small. Nevertheless, we further compare both models applied to 
the real data with these goals: 



1 . to investigate the importance of the precise value of spin in the line profile formation 
besides the relation between spin and the marginally stable orbit. 



2. to constrain the uncertainty in the spin determination by the less precise LAOR 
model. 



2.3. Comparison between the LAOR and KY models 



37 



2.3.1 Analysis with current data provided by XMM-Newton satellite 

For the comparison of the LAOR and kyrline model with current data, we chose two 
sources, MCG -6-30-15 and GX 339-4, which exhibited extremely broad iron lines (see 
Figures [Z4l and [231) . We re-analysed the observations performed by the XMM-Newton 
satellite (see Section l4.1.1l) . 



MCG-6-30-15 

The galaxy MCG-6-30-15 is a nearby Seyfert 1 galaxy (z = 0.0 08). The skewed iron 



line has been revealed in the X-ray spectra by all recent satellites (JGuainazzi et al.\ 11999 



Wilms g?a/.ll200 lb iFabian et a/.ll2002l : iLee et a/.ll2002l : iMiniutti et Q/.ll2007h . The XMM 



Newton observed MCG-6-30-15 for a long 350 ks exposure time during summer 2001 
(revolutions 301, 302, 303). We joined the three spectra into one for further analysis. 
Details of the data reduction and spectral re-analysis are described in Section |43l 

We used t he sam e continuum model for the MCG -6-30-15 spectrum as presented in 



Fabian et al.\ (|2002|) : the simple power law component absorbed by Galactic gas matter 
along the line of sight with column density nu = 0.41 x 10^' cm~^. The employed model 
is sufficient to fit the data above « 2.5 keV, which is satisfactory for our goal of the com- 
parison of the relativistic line models. The value of the photon index is F = 1 .90(1). The 
residuals are formed by a complex of a broad iron line and two narrow iron lines - one 
emission line at E = 6.4 keV likely originating in a distant matter (torus) and one absorp- 
tion line at £■ = 6.77 keV which can be explained by a blueshifted absorption originating 
in an outflow. The rest energy of the broad line is E = 6.7 keV, which corresponds to the 
helium- like ionised iron atoms. The spectral complexity in the line band allows an alterna- 
tive explanation - the model with two narrow emission lines at energies E = 6.4 keV and 
E = 6.97 keV. This alternative model leads to the presence of the broad line component 
atE = 6.4 keV. See Section |43] for details. 

A good fit of the broad line was found with a broken power law line emissivity with 
a steeper dependence on the radius in the innermost region. The goodness of the fit is 
constrained by the least squared method. The fit results in 2.5-9.5 keV are presented in 
Table fTJl The;^'^ values give comparable results for all employed models. T^i^^ X^sd ~ 
X^/v ~ 1.2, where v is the number of degrees of freedom which is related to the total 
number of energy bins and model parameters. The six independent parameters of the 
LAOR and KYRLINE models make the global minimum of ;t'^ rather wide with several 
local minima. Each model has a difl^erent tendency to converge to a different minimum. 
Hence, we did not compare only best fits of both models, but also the evaluated spin 
values by the kyrline and LAOR models when the other model parameters correspond 
to each other. The equivalent width of th e line is EW « 700 + 50 eV which is a rather 



high value but consistent with the result of iBrenneman & Reynoldsl (|2006|) . 

The errors presented in the table are evaluated while the other parameters of the model 
are fixed. However, the realistic errors are higher because the model parameters further 
depend on other parameters of the line and continuum models. To catch up these relations 
we produce various contour graphs focusing on the determination of the spin value, taking 



38 



CHAPTER 2. RELATIVISTIC LINES 



Table 2.1: Results of iron line models for MCG -6-30-15 in 2.5-9.5 keV. 



parameter 


KYRLINE 


KYRLINE* 


LAOR best 


LAOR loc.min. 


a/M 


94+0.02 

"•^^-0.03 


95+002 


98+"-°2 


0.96.0.01 


i [deg] 


26.7(7) 


31.5(7) 


35.7(5) 


26.8(5) 


E [keV] 


6.67(1) 


6.60(1) 


6.48(1) 


6.66(1) 


q\ 


4.9(1) 


3.7(1) 


4.8(1) 


4.7(1) 


qi 


2.84(4) 


2.11(4) 


2.50(3) 


2.87(3) 


n 


5.5(2) 


18.3(5) 


6.6(2) 


5.1(2) 


x'lv 


175/148 


174/148 


170/148 


174/148 



Notes: The errors in brackets are related to the last significant digit of the number, and they 
correspond to 90% confidence level calculated while the other model parameters are fixed. 



into account the other parameters of the used model, see Figure 12.71 The dependence of 
thcx^ value on the value of the spin (kyrline) or the inner disc radius (LAOR) is shown 
in the left column of the figure. The contour graphs for the spin and the inclination angle 
are shown in the middle column. The underlying model was fixed in both cases, and it 
was relaxed for the third evaluation in the right column where the contours for the spin 
and the power-law index are shown. The contour plots reveal quite complicated structure 
of the parameter space with several minima, making the spin estimation rather smeared. 
Taking all of these into account, we obtain for the spin value: 



Uky = OMllf, 



and aiaor = 0.96:°:°t 



Another issue is the smoothness of the contour plots. Some of the sharp or local 
features in the presented contour plots are due to computational problems of the fitting 
method but some of them can be smoothed by an improved choice of the internal com- 
putational parameters. We tried to change some of these parameters. The smoothest 
results were generally obtained when we increased number of fitting iterations together 
with wider steps of the value of the examined parameters. However, one may easily miss 
a true minimum when stepping a parameter value roughly which leads to bad results. 

To avoid this problem we always tried first to find a true minimum and then stepping 
the parameters in such a way that the values corresponding to the true minimum are not 
skipped. Finding the true minimum is rather difficult problem in such a rich space of 
free parameters. Therefore, we usually repeated the fitting procedure several times with 
different initial values of the model parameters and then we compared x^ values of the 
local minima. 



2.3. Comparison between the LAOR and KY models 
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Figure 2.7: The contour graphs show (from left to right) the dependence of the;^'^ value, 
the inclination angle and the power law index on the value of the spin (kyrline) or the 
inner disc radius (LAOR) for the MCG -6-30-15 spectrum in 2.5-9.5 keV. The x-axis is 
oppositely directed in the case of the inner disc radius of the LAOR model as x-variable 
for an easier comparison with the KYRLINE results. The black, red and green contours 
correspond to Icr, 2cr and 3cr, respectively. Top: The results of the LAOR model. Middle: 
The results of the kyrline model with limb darkening. Bottom: The results of the 
KYRLINE* model with limb brightening. 



GX 339-4 



The black hole binary GX 339-4 exhibited a strong broadened i ron line in the 76 k s ob- 
servation in 2002 when the source was in the very high state (IMiller et al\ 120041) . and 
also in two 138 ks observations in spring 2004 when the source was in the low-hard state 
(JMiller et a/.ll2006r) . See also Figure [Z4l adopted from the two papers. However, by our 
re-analysis of the data we found that the longer low/hard state observation is affected by 
a pile-up - see the details in Section I4.1.2[ Hence, we use further only the very high 
state observation from 2002. Details of the data reduction and spectral re-analysis are 
described in Section [4^2] 
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Table 2.2: Results of iron line models for GX 339-4 in 3-9 keV. 



parameter 


KYRLINE 


KYRLINE* 


LAOR 


a/M 


0.69!°:!^ 


0.62_oi4 


f\ 77+0.08 
"•''-0.12 


i [deg] 


19 + 3 


19 + 4 


17 + 4 


E [keV] 


6.97_o.2 


6.97_o.2 


6.97_o.2 


Q 


3.45 + 0.08 


3.35 + 0.08 


3.3 + 0.1 


X'lv 


147/125 


148/125 


148/125 



The fitting results of the relativistic line models in 3-9 keV are summarised in Ta- 
ble 12.21 and Figure 12. 8[ There are two minima found during the fitting procedure. We 
preferred the one which better corresponds to the results obtained by the independent 
radio and infrared measurements which constrained the inclination angle to be / < 26° 
(|Gallo et a/.ll2004|) . The dependence of the goodness of the fit on the spin value is shown 
in the left column of the figure. The contour graphs for the spin and the inclination angle 
are depicted in the middle column, and for the spin and the power law photon index in the 
right column of the figure. The derived spin value is then: 

a^y = 0.69!°;}^ and aiaor = 0.77!°:J°. 

In both cases, MCG -6-30-15 and GX 339-4, the fitting with the LAOR model resulted 
in slightly higher values for the spin, but consistent with the values of the kyrline model 
within the general uncertainties of the spin estimation. The spin value of MCG -6-30-15 
is better constrained thanks to its high value and more counts in the line than in the case 
of GX 339-4 (see Table Q. 

Table 2.3: Count rates of the XMM-Newton observations in 2 - lOkeV. 





MCG -6-30-15 


GX 339-4 


number of counts [xlO''] 


0.97 


3.56 


net counts/s 


4.98(1) 


1547(1) 


model counts/s 


5.02 


1546 


line counts/s 


0.20(1) 


5.1(1) 


line counts 


4.37 X 10^ 


1.15 X 10^ 



2.3. Comparison between the LAOR and KY models 
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Figure 2.8: The contour graphs show (from left to right) the dependence of the;t'^ value, 
the inclination angle and the power law index on the value of the spin (kyrline) or the 
inner disc radius (LAOR) for the GX 339-4 spectrum in 3-9 keV. The x-axis is oppositely 
directed in the case of the inner disc radius of the LAOR model as x-variable for an easier 
comparison with the KYRLINE results. The black, red and green contours correspond to 
Icr, 2cr and 3cr, respectively. Up: The results of the LAOR model. Middle: The results 
of the KYRLINE model with limb darkening. Down: The results of the kyrline* model 
with limb brightening. 



2.3.2 Simulated data of next generation X-ray satellites 



Further, we applied the LAOR and KYRLINE models on the data with significantly higher 
quality supposed to be achieved by on-coming X-ray missions. The presently planned 
International X-ray Observatory (IXO, see e.g. ' Sookbinden l20 1 (J ) arose from the merg- 
ing of the former XEUS and Constellation-X missions. Because the details of the IXO 
mission have not been fixed at the moment when this work originated, w e used a prelim- 
inary response matrix of the former XEUS mission (|Amaud et q/.I 120091) . We generated 
the data for a kyrline model with a rest energy of the line E = 6.4 keV and radial disc 
emissivity that follows a power law with the index q = 3. We re-binned the data in order 



42 



CHAPTER 2. RELATIVISTIC LINES 



Table 2.4: Results of the LAOR fit in 3-9 keV in the simulated spectra. 





MCG -6-30-15 (simul.) 


GX 339-4 (simul.) 


parameter 


KY value 


fitted LAOR value 


KY value 


fitted LAOR value 


^in 


2.00 


l-86t°^ 


3.39 


3 20+°°^ 


i [deg] 


26.7 


25.4(3) 


19 


18.6(2) 


£iine [keV] 


6.70 


6.71(1) 


6.97 


6.94(5) 


^1 


4.90 


4.51(3) 


3.45 


3.2(1) 


^2 


2.80 


2.76(2) 


- 


- 


rb 


5.5 


6.2(1) 


- 


- 


Kline 


8.7 X 10"^ 


9.0 X 10-5 


6.5 X 10-3 


6.6 X 10-3 


x'/y 




1364/873 




2298/873 



Note: Fiducial values of the KY model parameters which were used to the data simulation are 
shown in columns "KY value". The inner disc radius is determined from the relation rin = rms(a) 
(eq.O). 



to have a resolution of 30 eV per bin. We then fit the data in the 1-9 keV energy range 
with the LAOR model using the same initial values of the fitting parameters as for the data 
simulation. Due to insufficient resolution of the LAOR model, a significant problem ap- 
pears at the high-energy edge of the broad line. This occurs because the next generation 
instruments achieve much higher sensitivity in the corresponding energy range. 

We examined the artificial data for a grid of values of the angular momentum and the 
inclination angle. The results for different values are shown in Figure [Z9l We excluded 
the higher-energy drop of the lines in order to reveal the differences in the overall shape 
of the line. It is clearly seen from the figures that the effect of the spin on the shape of the 
line is sufficiently resolved with the higher quality data. The contour graphs of the inner 
disc radius and the inclination angle are shown in Figure [2?T0l for three different values of 
the spin, a = 0, a = 0.7, a = 0.9982. For lower values of the spin, the results obtained 
with the LAOR model differ more from the fiducial values. The LAOR model tends to 
underestimate the inner disc radius, or, by other words, to overestimate spin value related 
to the inner disc radius by the relation a = a (r^s = ''in)- 

Further, we produced simulated data for the Seyfert galaxy MCG-6-30-15 and the 
black hole binary GX 339-4 using rather simplified models which were suitable to fit 
the current XMM-Newton data. For MCG-6-30-15 we used a power law model plus a 
KYRLINE model for the broad iron line, absorbed by Galactic gas matter along the line 
of sight: PHABS * (POWERLAW -I- KYRLINE). The parameters of the continuum are the 
column density n^ = 0.4 x 10^' cm"^, the photon index F = 1.9 of the power law, and its 
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normalisation ^r = 5 x 10 "*. The values of the line parameters are summarised in the KY 
value column of Table [Z4l The exposure time was chosen as 220 ks and the flux of the 
source as 1.5 x 10^'^erg cm~^ s~' in the 2-lOkeV energy range (i.e. 1.4 x 10' cts). 

For GX 339-4 we used phabs * (diskbb + powerlaw + kyrline) with n^ = 0.6 x 
10^2 cm-2, kTin = 0.87 keV (Ky^j = 1.4 x 10^), and T = 3 (Kr = 5.6). The exposure time 
was chosen as 75 ks, the flux of the source as 9.3x 10"^erg cm~^ s~' in the 2-lOkeV energy 
range, i.e. 1 .5x 10^ cts. The number of counts is more than two orders of magnitude higher 
than for the observation of the XMM-Newton satellite (see Table 12.31) . The reason is due 
to the loss of 97% of the photons during the burst mode of XMM-Newton observation 
which eliminates the pile-up problem (see Section I4.1.1|) . The next generation X-ray 
missions are supposed to have a calorimeter instead of the CCD camera on-board which 
will get rid of such problems. 

The results of the LAOR fit are shown in Table 12.4] and in Figure IXTTl The broad iron 
line component of the model is plotted in the top panel of the figure. The model continuum 
components are not displayed there to clearly see the deflections of the LAOR model. The 
most prominent discrepancy appears at the higher-energy drop, which is clearly seen in 
the data/model ratio plot in the middle panel of Figure 12.111 The model parameters are 
constrained with small error bars (see contours a vs. i in the bottom panel of Figure IXTTI ), 
which reveals a difference between the kyrline and the LAOR models. 

The spin value derived from the analysis using the LAOR model is: 

n - n QS 8+0003 

«laor,MCG - U-yjO-0.004 

for MCG -6-30-15, while the fiducial value of the spin was afid,MCG = 0.940, and 

- n 7/L+002 
«laor,GX - U. /4_o.02 

for GX 339-4, while the fiducial value of the spin was flfid.ox = 0.70. 



2.4 Discussion of the results 

We investigated the iron line band for two representative sources - MCG -6-30-15 (active 
galaxy) and GX 339-4 (X-ray binary). The iron line is statistically better constrained 
for the active galaxy MCG -6-30- 15 due to a significantly longer exposure time of the 
available observations - see Table [23] for comparison of count rates of the observations. 
The spectra of both sources are well described by a continuum model plus a broad iron 
line model. 

We compared modelling of the broad iron line by the two relativistic models, LAOR 
and KYRLINE. The LAOR model always assumes an extreme Kerr metric (a/M = 0.9982), 
and therefore, it is used for the fitting of different spin values only by identifying the inner 
edge of the disc with the marginally stable orbit. In the kyrline model, on the other 
hand, the spin itself is a fitting parameter and the metric is adjusted to the actual value of 
a/M. 
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The discrepancies between the KYRLINE and laor results are within the general un- 
certainties of the spin determination using the skewed line profile when applied to the 
current data. This means that the spin is currently determined from the position of the 
marginally stable orbit rather than from the overall shape of the line. 

Results with both models are apparently distinguishable for higher quality data, as 
those simulated for the next generation X-ray missions. We find that the LAOR model 
tends to overestimate the spin value and moreover, it has insufficient energy resolution 
which aff"ects the correct determination of the high-energy edge of the broad line. This 
leads to large ;kf^ values by fitting the high resolution data (see Table 12.41) . 

Technically, the kyrline model leads to a better defined minimum of ;^f^ for the best 
fit value. The confidence contour plots for a/M versus other model parameters are more 
regularly shaped. This indicates that the kyrline model has a smoother adjustment 
between the different points in the parameter space allowing for more reliable constraints 
on a/M. The LAOR model has a less accurate grid. Smoother and more precise results 
by the kyrline model are at the expense of the computational speed. Empirically, we 
found that the kyrline model is about lOx slower than the laor model. 
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Figure 2.9: Simulated data for three different values of spin (from top to bottom: a = 0, 
0.7, 0.9982) and two different inclination angles 6 = 30 deg (left), and 6 = 60 deg (right). 
The artificial data are shown with the expected errors (black crosses). Black curve is the 
seed KYRLINE model from which the data were generated. The red curve is the LAOR 
model which resulted from the fitting of the simulated data. Each panel also contains a 
plot with ratios of the data to the LAOR model. 
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Simulated XEUS data for kyrline model with a/M = 0, theta = 30 deg 
fit with the laor model 



Simulated XEUS data for kyrline modei with a/M - 0, theta - 60 deg 
fit with the laor model 
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Simulated XEUS data for kyriine model with a/M - 0.7, theta - 30 
fit with the laor model 
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Simulated XEUS data for kyriine model with a/M - 0.9982, theta = 30 
fit with the laor model 
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fit with the laor model 





nner disk radius [gravit. radii] 



inner disk radius [gravit. radii] 



Figure 2.10: Contour graphs of the inclination angle and the inner disc radius of the 
LAOR model applied to the same artificial data as in Figure [Z9l The fiducial values of the 
studied parameters are i = 30 deg (left), i = 60 deg (right), and rin = r^sia = 0) = 6 rg 
(top), 3.39 Tg (middle), and 1.235 r^ (bottom). The individual curves correspond to 1 cr, 
2 cr, and 3 cr, respectively. 
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Figure 2.11: The simulated spectra for MCG -6-30-15 {left) and GX 339-4 {right) using 
the preliminary response matrix of the XEUS mission. Top: Broad iron line generated 
by the kyrline model (black data) and fitted by the LAOR model (red curve). Middle: 
Data/model ratio when fitting by the LAOR model. Bottom: Contours for the inclination 
angle and the inner disc radius of the LAOR model. The fiducial values are rin = 2 and 
e = 26.7 for the case of MCG -6-30-15 {left), and n^ = 3.39 and 0=19 for the case of 
GX339-A {right). 



Chapter 3 

Role of the emission directionality in the 
spin determination 

3.1 Emission directionality 

Light rays coming from a black hole accretion disc are highly curved in strong gravita- 
tional field and the emission angle given by eq. (12.81) may possess any value, depending 
on the position on the disc from which the radiation is reflected. Figure 13.11 shows the 
contours of the constant local emission angle 6^ from the accretion disc around a non- 
rotating, a = 0, and a maximally rotating, a = 1, black hole taking into account their 
distortion by the central black hole and assuming that the emitted photons reach a distant 
observer at a given view angle Oo- It clearly illustrates that the local emission angle spans 
the entire range, from to 90 degrees for any value of Og. This is due to the combined 
eff'ects of aberration and light bending which grow greatly near the inner rim of the disc. 
In the inner region, the photons are boosted in the direction of rotation and they emerge 
along grazing angles. Notice for the case of a higher disc inclination, 6o = VOdeg, that 
there is a small region behind the black hole (below, in the picture) where the local emis- 
sion angle is very low, i.e. the radiation is emitted in the direction almost perpendicular 
to the disc. This is the clear eff'ect of the light bending. Although the contours are most 
dramatically distorted near the horizon the emission angle is visibly diff"erent from the 
observer inclination even quite far from the horizon, at a distance of several tens rg. This 
is mainly due to special-relativistic aberration which decays slowly with the distance as 
the disc obeys Keplerian rotation at all radii. Asymptotically, 9e(r, <p) -^ B^ for r ^ oo. 

As emission directionality, we call the dependence of the intensity on the emission 
angle. Because of the variety of the emission angle values, the directional distribution of 
the outgoing radiation is among the important aspects that must be addressed. This func- 
tionality depends on how the radiation originates and also on the surrounding conditions. 
In the case of a black hole accretion disc, the directionality of thermal and reflection ra- 
diation may be significantly different. In general, the directionality is a function of radius 
from the black hole and energy of the emitted radiation. In practical application, however, 
a unique profile is standardly assumed, invariably over the entire range of radii in the disc 
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Figure 3.1: Contours of the local emission angle, 6e(r,(p) near a non-rotating 
Schwarzschild black hole, a = 0, (top), and a maximally rotating black hole, a = 1, 
(bottom), depicted in the equatorial plane ix,y). The black hole and the accretion disc 
rotate counter clock- wise. A distant observer is located towards the top of the figure. The 
inner region is shown up to r = 20 gravitational radii from the black hole. The black hole 
is denoted by a dark filled circle around the centre, for Schwarzschild black hole the circle 
around represents the marginally stable orbit. Two cases of different observer inclinations 
are shown, ^o = 30 deg (left) and 6o = 70 deg (right). The colour bar encodes the range 
of 6f.{r, (f) from to 90 degrees. The emission angle is measured from the disc normal 
direction to the equatorial plane, in the disc co-moving frame. 
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and the energy in the spectral band. Particularly, the limb darkening law is frequently 
used since it is embedded in the currently most widely used LAOR model (Sect |2.3l) . 

Limb darkening traditionally refers to the gradual diminution of intensity in the image 
of the surface of a star as one moves from the centre of the image to the edge. It is a 
consequ ence of uneven angular distribution o f the radiation flux emerging from the stellar 



surface (IChandrasekhad 1 1 960l : iMihalasI 1 1 97 8[) . Limb darkening results as a combination 



of two eff"ects: (i) the density of the surface layers decreases in the outward direction, 
and (ii) temperature also drops as the distance from the centre of the star increases. The 
outgoing radiance is therefore distributed in a non-uniform manner. 

The limb darkening law is widely applied also to describe radiation coming from an 
accretion disc around a black hole. Nevertheless, the actual form of the angular distribu- 
tion depends on the physical mechanism responsible for the emission and the geometrical 
proportions of the source. The limb darkening law is relevant for thermal radiation of 
the disc with infinite optical depth. However, already when the optical depth is assumed 
to be finite and also internal heating is taken into a ccount the emission directionality 



significantly changes (see e.g. iFukue & Akizukill2006|) . The case of reflection is substan 



tially diff"erent from the thermal radiation because the initial radiation comes from above 
the disc and passes through the disc atmosphere. The higher layers are characterised by 
higher ionisation and temperature with the electron scattering as the dominant process 
there. The observer looks deeper to the colder layers of the atmosphere when the emis- 
sion angle is low. However, the radiation is there more likely to be absorbed resulting 
in the detection o f weaker reflection rad iation than when the emission angle is higher 
(see discussion in Goos mann et q/.II2007L Sect. 3). This eff'ect is just opposite to the limb 
darkening, and therefore, it is called as limb brightening. 

Given the high velocity of the orbital motion (see Fig. II. II) and the strong-gravity light 
bending near the black hole, the eff'ect of directional anisotropy of the local emission is 
enhanced. As a result of this interplay between the local physics of light emission and 
the global eff"ects of the gravitational field, different processes contribute to the final (ob- 
served) directional anisotropy of the emission. For this reason it is important to describe 
the angular distribution in a correct manner; ad hoc choices of the limb-darkening law 
may lead to errors in the determination of the model best fit parameters, including the 
inaccuracy in the spin parameter which are difficult to control, or they may prevent us 
from estimating the statistical confidence of the model. 

3.1.1 Effects of the emission directionality on the iron line profiles 

In optics, Lambert's cosine law describes an emitter producing a radiation intensity that 
is directly proportional to /ig = cosSg- Lambertian surfaces exhibit the same apparent 
radiance when viewed from any angle 6o. Likewise, Lambertian scattering refers to the 
situation when the surface radiates as a result of external irradiation by a primary source 
and the scattered light is distributed according to the same cosine law. This is, however, 
a very special circumstance; directionality of the emergent light is sensitive to the details 
of the radiation mechanism. For example, the classical result of the Eddington approxi- 
mation for stellar atmospheres states that the effective optical depth of the continuum is 
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T = |, and so the emergent intensity is described by the limb-darkening law, I (pie) oc //g-i-|. 

In the case of a fluorescence iron line produced by an il luminated plane-parallel slab, 
the angular distribution was investigated by va rious authors iBasko 1978; George & Fabian| 
ri991. : .Haardt,1993, : .Ghisellini. Haardt & Matt.l994 : .Revnolds. Nowak & Malonev.200Q) . 
In that case a complicated interplay arises among the angular distribution of the primary 
irradiation, reflection and scattering in the disc atmosphere. Several authors pointed out 
that it is essential for the reliable determination of the mod el parameters to determine 
the angular directionality of the broad line emission correctly. iMartocchia. Karas & Matt 
(1200 0) noted, by employing the lamp-post model, that "...the broadening of the observed 
spectral features is particularly evident when strongly anisotropic emissivity laws, result- 
ing from small h [i.e., the lamp-post elevation above the equatorial plane], are consid- 
ered." 

The important role o f the emission angular directionality was clearly spelled out by 



Beckwith & Pond (|2004) : "...the angular emissivity law (limb darkening or brightening) 



can make significant changes to the derived line profiles where light bending is impor- 
tant" (see their Fig. 9-13). Similarlv, IPovciak et al.\ (l2004ah : IPovciak. Karas & Yaqoob 



(|2004l) and lBeckwithI (|2005b compared the relativistic broad lines produced under differ- 
ent assumptions about the emission angular directionality. However, to verify the real 
sensitivity of the models to the mentioned effect of directionality, it is necessary to con- 
nect the radiative transfer computations with the spectral fitting procedure, and to carry 
out a systematic analysis of the resulting spectra, taking into account both the line and the 
continuum in the full relativistic regime. We report on our results from such computations 



in Section 3.2.4 



Reynolds et a/.l(|2004l sec. 4.3) argue that the combined effect of photoelectric absorp- 



tion in the disc and Compton scattering in the corona more affect the iron line photons 
emerging along grazing light rays than continuum photons. They conclude that the line 
equivalent width should be diminished for observers viewing the accretion d isc at high 



inclination angles. Such a trend can b e seen also in the lamp-post model of iMatt g?a/. 



(|l992 h and^ Martocchia. Karas & Matt! (12000. see their Figure 11). However, in the latter 
work this diminution is less pronounced when we compare it with the case of intrinsically 
isotropic emissivitV; 



More recently, iNiedzwiecki & Zyckil (|2008|) studied the effect of different limb dark- 
ening laws on the iron line profiles. They pointed out that the role of emission direction- 
ality can be quite significant once the radial emissivity of the line is fixed with sufficient 
confidence. However, this is a serious assumption. In reality, the radial emissivity is not 
well constrained by current models. 

The angular dependence of the outgoing radiation is determined by the whole inter- 
connection of various effects. We describe them in more detail below (Sec. 13.2.41) . Briefly, 
the conclusion is such that realistic models require numerical computations of the full ra- 
diative transfer. We have developed a complete and consistent approach to such radiative 
transfer computations in the context of the broad iron-line modelling together with the 
underlying continuum computations. As described in considerable detail below, we per- 
formed the extensive computations which are necessary in order to reliably determine the 
impact of the emission angular anisotropy on spectral fitting results (namely, on the deter- 



52 CHAPTER 3. ROLE OF THE EMISSION DIRECTIONALITY . . . 



mination of the black hole angular momentum). In particular, we describe detailed results 
from the investigation of the model goodness (by employing an adequate statistical anal- 
ysis of the complicated ;^^^ parameter space). Such the analysis has not been performed so 
far in previous papers because detailed stepping through the parameter space and proper 
re-fitting of the model parameters was not possible due to the enormous complexity of the 
models and extensive computational costs. 

In regular stars and their accretion discs, the relativistic effects hardly affect the emerg- 
ing radiation. The situation is very different in the inner regions of a black hole accretion 
disc, where the energy shift and gravitational lensing are significant. The observed signal 
can be boosted or diminished by the Doppler effect combined with gravitational redshift: 
lijio) I lip-s) = g^ (eg. 12. 161 ). where the g values span more than a decade (s ee Fig. 12. II) . 



As mentioned above, many authors have adopted the defining choice (lLaoilll991r) of 
the cosine profile for the line angular emissivity: A\{p.s,Es) = 1-1- 2.06 /ig- This relation 
describes the energy-independent limb-darkening type of profile. However, the choice 
is somewhat arbitrary in the sense that the physical assumptions behind this law are not 
satisfied at every radius over the entire surface of the accretion disc. It has been argued 
that the limb-darkening characteristics need to be modified, or even replaced by some 
kind of li mb br ightenin g in the case of X-ra y irradiated disc atmospheres with Compton 



reflection ( lHaardt,1993; : lGhisellini. Haardt & Matt., 1 994 : IZycki & Czemyiri994h . The lat 



ter should include the energy dependence, as the Compton reprocessing of the reflected 
component plays a significant role. The angle-de pendent computatio ns of the Compton 
reflection demonstrate these effects convincingly (ICzemy et a/.ll2004t) . Indeed, the same 



effect is seen also in our computations, as shown in the right panel of Figure 13.21 The 
increase of emissivity with the emission angle strongly depends on th e ionisation state 



of the reflecting material, so the actual situation can be quite complex (|Goosmann et al. 



2007 ) 



A question arises of whether the current determinations of the black hole angular mo- 
mentum might be affected by the uncertainty in the actual emissivity angular distribution, 
and to what degree. In fact, this may be of critical concern when future high-resolution 
data become available from the new generation of detectors. We have therefore carried 
out a systematic investigation using the KY code to reveal how sensitive the constraints 
on the dimensionless a parameter are with respect to the possible variations in the angular 
part MQus) of the emissivity. 

Figures [2n and [3TT] demonstrate the main attributes of the photon propagation in black 
hole space-time relevant to our problem: the energy shift and the direction of emission 
depend on the view angle of the observer as well as on the angular momentum of the 
black hole. Notice that, near the inner rim, the local emission angle is indeed highly 
inclined towards the equatorial plane where at the same time the outgoing radiation is 
boosted. These effects are further enhanced by gravitational focusing, which we also take 
into account in our calculations. 

It should be noted that the energy shifts and emission angles near a black hole have 
been studied by a number of authors in mutually complementary ways. The figures shown 
here have been produced by plotting directly the content of FITS format files that are 



encoded in the corresponding KY routines (see lDovciakI (|2004|) . where an atlas of contour 



3.1. Emission directionality 



53 



4.5 

4 

3.5 

3 I- \ 
2.5 

2 
1.5 - 

1 

0.5 




limb brightening 



limb darkening 



isotropic 



0.2 



0.4 0.6 




Figure 3.2: Left: directional distribution of the intrinsic emissivity (see Sec. 13.2. ll for de- 
tails). Right, an example of the directional distribution from the numerical computations, 
showing results from the reprocessing model with the continuum photon index Y = 1.9, 
integrated over 2-100 keV energy range (solid line), and in several different energy sub- 
ranges (line styles are indicated in the inset). The latter graph demonstrates the presence 
of limb-brightening effect in comparison with respect to the isotropic emission. This ef- 
fect is clearly visible as a function of the emission angle /ig, although its magnitude is 
smaller than the limb-brightening approximation used in the left panel. See sec. 13.2. 4] for 
details. 



pl ots is presented for differe nt inclinations and spins). Analogous figures were also shown 



m lBeckwith & Pond (|2004|) . who depicted the dependence of the cosine of the emission 
angle on the energy shift of the received photon and the emission radius. 

The above quoted papers concentrated mainly on the discussion of the energy shifts 
and the emission directions of the individual photons, or they isolated the role of relativis- 
tic effects on the predicted shape of the spectral line profile. They clearly demonstrated 
that the impact of relativistic effects can be very significant. However, what is still lacking 
is a more systematic analysis which would reveal how these effects, when integrated over 
the entire source, influence the results of spectra fitting. To this end one needs to perform 
an extensive analysis of the model spectra including the continuum component, match the 
predicted spectra to the data by appropriate spectra fitting procedures, and to investigate 
the robustness of the fit by varying the model parameters and exploring the confidence 
contours. 

One might anticipate the directional effects of the local emission to be quite unimpor- 
tant. The argument for such an expectation suggests that the role of directionality should 
grow with the source inclination, whereas the unobscured Seyfert 1 type AGNs (where the 
relativistically broadened and skewed iron line is usually expected) are thought to have 
only small or moderate inclinations. However, this qualitative trend cannot be used to 
quantitatively constrain the model parameters and perform any kind of precise analysis, 
needed to determine the black hole spins from current and future high-quality data. Such 
an analysis has not been performed so far, and we embark on it here for the first time. 
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3.2 Iron Ka line band examined with different direction- 
alities 

3.2.1 Approximations to the angular emission profile 

We describe the methodology which we adopted in order to explore the effects of the 
spectral line emission directionality. To this end we first employ simple approximations, 
neglecting any dependence on the photon energy and the emission radius. We set the line 
intrinsic emissivity from the planar disc to be described by one of the following angular 
profiles, 

Casel: f ln(l +ju-') (Haardt 1993) 

Case 2: M(jUs) = • 



Case 3: 



1 (locally isotropic emission) (3.1) 

l+2.06;Ue (Laor 1991) 

The three cases correspond, respectively, to the limb-brightened, isotropic, and limb- 
darkened angular profiles of the li ne emis s ion. 



The limb-brightening law by iHaardti (|1993|) describes the angular distribution of a 



fluorescent iron line emerging from an accretion disc that is irradiated by an extended 
X-ray source. The relation was obtained from geometrical considerations and agrees wel l 
with more detailed Monte-Carlo computations (IGeorge & Fabianlll99ll : lMatt et a/.il992|) . 



The physical circumstances relevant for the limb-darkening law are different, and we 
include this case mainly because it is implemented in the laor model and frequently used 
in the data analysis. The isotropic case, dividing all limb-brightening and limb-darkening 
emissivity laws, is included in our analysis for comparison. 

The radial profile of the emission is set to a unique power law, eq. (12.181 ). over the 
entire range of radii across the disc. The directionality formula (|3.1I) of the intrinsic emis- 
sivity and the resulting spectral profiles are illustrated in Fig. 13.21 (left panel). Naturally, 
more elaborate and accurate approximations h a ve bee n discussed in the literature for some 



time. For example, iGhisellini. Haardt & MattI (|1994|) in their eq. (2) include higher-order 



terms in /ig to describe the X-ray reprocessing in the single- scattering Rayleigh approx- 
imation. However, at this stage the first-order terms are sufficient for us to demonstrate 
the differences between the three cases. Later on we will proceed towards numerical radi- 
ation transfer computations that are necessary to derive realistic profiles of the emission 
angular distribution and to keep their energy dependence. 

The dependence of the spectral profiles on the angular distribution M(pe) of the in- 
trinsic emission is shown in figures [331 - 13. 4[ It is apparent from Figure [33] that the limb 
brightening case makes the line profile broader (in the left panel) and the height of the blue 
peak lower (in the right panel) than the limb darkening case for the same set of parameters, 
which can consequently lead to discrepant evaluation of the spin. While Figure [33] shows 
the line profile for an extended disc. Figure \TM deals with a narrow ring. In this case, a 
typical double-horn profile develops. Although the energy of the peaks is almost entirely 
insensitive to the emission angular directionality, the peak heights are influenced by the 
adopted limb darkening law. In the right panel of Figure 13 .41 the ratio of the two peaks 
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Figure 3.3: Left: theoretical profiles of the relativistic line (the KYRLINE model, without 
continuum), corresponding to the three cases in the left panel of Fig. 13.21 The lines are 
normalised with respect to the height of the blue peak. Model parameters are a = 0.9982, 
q = 3, Tin = Vs(a) = 1.23, rout = 400, 6^ = 30 deg, Eq = 6.4 keV. Right: the same as in the 
left panel, but with the normalisation set in such a way that the radiation flux is identical 
in all three profiles. 
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Figure 3.4: Left: The same as in Fig. 13. 3[ but for a narrow ring with r^^ = 4.7, rout = 4.8 
(other parameters are a = 0.9982, q = 3, 9o = 30 deg, Eq = 6.4 keV). The lines are 
normalised with respect to the height of the blue peak. Right: The ratio of the two peak's 
heights for the three cases of the emission angular directionality according to eq. (13.11) 
and for the two values of the spin and the inclination angle (crosses and circles are for 
^o = 30 deg; empty and filled triangles are for do = 60 deg). 



is shown for the three cases of the emission angular directionality according to eq. (13.11) . 
The influence of the directionality is apparent and it is comparable to the effect of the 
inclination angle. The spin value has only little impact for large inclinations. 

In general, we notice that the extension and prominence of the red wing of the line 
are indeed related to t he int rinsic emission directionality. This was examined in detail 
by iBeckwith & Pond (120041) . who plotted the expected broad-line profiles for different 
angular emissivity and explored how the red wing of the line changes as a result of this 



56 



CHAPTER 3. ROLE OF THE EMISSION DIRECTIONALITY . . . 



0.01 




-a 
o 

E 



1.01 



0.99 



.ityttiyrimTi lid "11' 





4 6 

energy [keV] 



4 6 

energy [keV] 



Figure 3.5: Left: Unfolded spectrum generated by the powerlaw + kyrline model 
with parameters of the line identical as in the Figure 13.31 We assumed the limb- 
brightening profile, i.e. Case 1 of eq. (13.11) . Right. Ratio of the simulated data on the 
left to the same model but with the limb-darkening profile (Case 3). The normalisation of 
the line was allowed to vary during the fit. The;^^^^^ changed from 1.04 to 1.30 between the 
two cases. The details about the simulation of the data points are the same as described 
in the section [J.2. 31 



undetermined angular distribution. Similarly, the inferred spin of the black hole must 
depend on the assumed profile to a certain degree. This is especially so because the spin 
is determined by the extremal redshift of the red wing of the line (see Section O. 

However, simple arguments are insufficient to assess how inaccurate the spin deter- 
mination might be in realistic situations, as the spectral fitting procedure involves several 
components extending over a range of energy above and below the iron-line band. We 



illustrate this in Figure 13.51 noticing that the diff"erent prescriptions produce very similar 
results outside the line energy, but they do diff"er at the broad line energy range. The the- 
oretical (background-subtracted) profiles of the relativistic line cannot alone be used to 
make any firm conclusion about the error of the best fit parameters that could result from 
the poorly known angular emissivity. To this end one has to study a consistent model of 
the full spectrum. With a real observation, the sensitivity to the problem of directionality 
(as well as any other uncertainty inherent in the theoretical model) will depend also on 
the achieved resolution, energy binning and the error bars of the data. 



3.2.2 Example: directionality effects in MCG -6-30-15 

We re-analysed a long XMM-Newton observation of a nearby Seyfert 1 galaxy MCG -6- 
30-15 to test whether the diff"erent directionality approximations can be distinguished in 
the current data. The observation took place in summer 2001 and the acquired exposure 
time was about 350 ks. We reduced the EPN data from three sequential revolutions (301, 
302, 303). The details of the reduction of the data are des cribed in Section 14.31 

We applied the same continuum model as presented in lFabian et all (120021) : the power 
law component (photon index F = 1.9) absorbed by Galactic gas (column density n^ = 
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Figure 3.6: The best-fit values of ;^f^ statistics for the spin parameter, which we obtained 
by grad ually stepping a fro m 0.8 to 1. XMM-Newton data for MCG -6-30-15 were em- 
ployed (IFabian et q/.II2002|) . Left: the limb-brightening profile (Case 1). Right, the limb 
darkening profile (Case 3). The dashed line is the 90% confidence level (2 sigma). See 
the text for a detailed description of the model. 



0.41 X 10^^ cm"^). This simple model is sufficient to fit the data above « 2.5 keV, which 
is also satisfactory for our goal of reproducing the overall shape of the broad iron line. 
Other components need to be added to the model in order to fully understand the spectrum 
formation and to decide between viable alternatives, in cluding the presence of outflows 
and a combination of absorption and reflection effects (|Miller. Turner & Reevesll2008al : 
Turner & Milleiil2009|) . However, our goal here is not to give precedence to any of the 
particular schemes. Instead, we rely on the model of a broad line and we test and compare 
different angular laws for the emission. 

The residuals from the simple power law model can be explained by a complex of 
a broad line and two narrow lines. Due to the complexity of the model, it is hard to 
distinguish betwe en the narrow abso rption line at E x 6.75 keV and the emission line 
at £■ « 6.97 keV (JFabian et a/.ll2002b . Although the level of ;^^ values stays almost the 
same, the parameters of the broad iron line do depend on the continuum model and the 
presence of narrow lines (see Section |43] ). By adding an absorption line at E ^ 6.75 keV, 
the best fit rest energy of the broad component comes out to be £ = 6.7 keV, and so 
the broad line originates in a moder ately ionised disc. This result is consistent with 
BaUantvne. Vaughan & FabianI (l2003h . 

Figure 13.61 shows the results of the one-dimensional steppar command performed in 
XSPEC, which demonstrates the expected confidence of a-parameter best-fit values for 
the two extreme cases of directionality. Case 1 and Case 3. The model used in the XSPEC 
syntax is: phabs*(powerlaw -i-zgauss -i- zgauss -i-kyrline). The fixed parameters 
of the model are the column density riu = 0.41 x 10^^ cm"^, the photon index of power 
law r = 1.9, the redshift factor z = 0.008, the energy of the narrow emission line ^em = 
6.4 keV, and the energy of the narrow absorption line ^abs = 6.77 keV. The parameters 
of the broad iron line were allowed to vary during the fitting procedure. Their default 
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Table 3.1: The best- fit spin values inferred for the three cases of the limb darken- 
ing/brightening law, eq. (|3.1I) . for the kyrline model. 



Case 


a{ = 


0.7 


a{ = 


0.9982 


no. 


Oi = 30° 


Of = 60° 


Oi = 30° 


0f = 60° 


1 


56+°*"* 

^•-^"-0.03 


69+"°^ 


92+"°3 

^•^^-0.03 


981+0013 

^•'^° ^-0.031 


2 


66+°"^ 


70+"°2 


> 0.966 


> 0.986 


3 


74+0.05 

^•'^-0.03 


70+"°^ 

^•'^-0.03 


> 0.991 


> 0.993 



The artificial data were generated using the KYRLINE model with isotropic directionality 
and the fiducial values of parameters (denoted by the subscript "f"). See the main text for 
details. 



values were £'broad = 6.7 keV for the energy of the broad iron line, ^o = 30deg for the 
emission angle, qi = 4.5, ^2 = 2 and r^ = 10 rg for the radial dependence of the emissivity 
(the radial part of the intensity needs to be rather complicated to fit the data and can be 
expressed as a broken power law: "Rirg) = r''" for r^ < r^, and 'R(rs) = r~''^ for r^. > r^). 

The determined best-fit values for the spin are virtually the same for both cases, 
independent of the details of the limb-brightening/darkening profile. However, this re- 
sult arises on account of the growing complexity of the model. The differences be- 
tween the two cases become hidden in different values of the other parameters - espe- 
cially in qi, q2 and r^, i.e. the parameters characterising the radial dependence of the 
line emissivity in KYRLINE as a broken power law with a break radius r^,. We find: (i) 
Abroad = 6.60(1), 00 = 31.5(7)deg, qi = 3.7(1), q,, = 2.1(1), r^ = 18(l)rg for Case 1; 
and (ii) Abroad = 6.67(l)keV, O^ = 26.7(7)deg, qi = 5.3(1), qj = 2.8(1), r^, = 4.9(2) r^ 
for Case 3. The errors in brackets are evaluated as the 90% confidence region for a single 
interesting parameter when the values of the other parameters are fixed. The combina- 
tion of three parameters q^, qi, r^ thus adjusts the best fit in XSPEC. Nonetheless, the 
clear differences between the models occur consistently with theoretical expectations: for 
Case 3 the lower values of the spin, a < 0.87, produce larger ;^^ and the best-fit spin can 
reach the extreme value within the 90% confidence threshold. 



3.2.3 Analysis of simulated data for next generation X-ray missions 

In order to evaluate the feasibility of determining the spin of a rotating black hole and 
to assess the expected constraints from future X-ray data, we produced a set of arti- 
ficial spectra. We used a simple model prescription and preliminary response matri- 
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Figure 3.7: Test results of the theoretical fits with cif = 0.7, 6{ = 30 deg and three different 



profiles of the emission directionality - left: limb brightening (Case 1 in eq. 13.11) . middle: 
isotropic (Case 2), right: limb darkening (Case 3). The simulated data were generated 
using the powerlaw + KYRLINE model with isotropic directionality. Top: dependence 
of the best fit;^^ values on the spin value. The horizontal (dashed) line represents the 90% 
confidence level. Bottom: contour graphs of a versus Oo- The contour lines correspond to 
1, 2, and 3 sigma. The position of the minimal value of ;^'^ is marked with a small cross. 
The values of ;^^ corresponding to the minimum and to the contour levels are shown at 
the top of each contour graph. The large cross indicates the position of the fiducial values 
of the angular momentum and the emission angle. Besides the values of spin, inclination 
and normalisation constants of the model, the other parameters were kept fixed at their 
default values: T = 1.9, ^o = 6.4 keV, rjn = Vs, /"out = 400, q = 3, (default values of 
normalisation constants: Ky = 10"^, ^une = lO"'*). 



ces for the International X-ray Observatory (IXO) missionlll Here we limit the energy 
band in the range 2.5-10 keV. The adopted model consists of powerlaw for contin- 
uum (photon index F and the co rresponding normalisation Ky), plus kyrline model 
(IDovciak. Karas & YaqoobI I2004D for the broad line. The normalisation factors of the 



model were chosen in such a way that the model flux matches the flux of MCG -6-30-15. 
In this section, the simulated flux is around 3.1 x 10"'^ ergcm"^ s"^ with about 3% of the 
flux linked to the broad iron line component. The simulated exposure time was 100 ks. 

We generated a set of "fake" spectra (i.e., artificial spectra in the XSPEC terminol- 
ogy). These spectra were produced in a grid of angular momentum values while assuming 
isotropic directionality. Case 2 in eq. (13.11 ). We call the assumed angular momentum of 
the black hole the fiducial spin and we denote it as af. We performed the fitting loop to 



'We used the current version of provisional response matrices available at 
http://ixo.gsfc.nasa.gov/science/responseMatrices.html for the glass core calorimeter, dated Octo- 
ber 30, 2008. We used the energy resolution of 5 eV per bin. 
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Figure 3.8: The same as in Figure [3771 but for a^ = 0.7 and 6{ = 60 deg. 
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Figure 3.9: The same as in Figure [3771 but for Uf = 0.998 and Of = 30 deg. 



these data points using each of the three angular emissivity profiles. Once the fit reached 
convergence, we recorded the inferred spin a. Figures [3771 - 13 . 1 01 show the results in terms 
of best-fit ;\f^ profiles and the confidence contours for two (MSQXQni fiducial values of the 
spin (we assumed at = 0.7, and 0.998). We summarise the values for the inferred spin for 
two inclination angles i = 30 deg and i = 60 deg in Table 13. ll 

The fitting procedure was performed in two diff"erent ways - having the rest of the 
parameters free or keeping them frozen. Obviously the former approach results in an 
extremely complicated ;^f^ space. Therefore, for simplicity of the graphical representation, 
we plot only the results of the second approach which, however, gives broadly consistent 
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Figure 3.10: The same as in Figure [3771 but for a^ = 0.998 and d{ = 60 deg. 



results (though it misses some local minima of ^^). In other words, the plots have the 
parameters of the power law continuum, the energy of the line and the radial dependence 
parameter fixed at F = 1.9, Eq = 6.4 keV, and q = 3. 

The conclusion from this analysis is that the determination of a indeed seems to be 
sensitive within certain limits to the assumed directionality of the intrinsic emission. The 
suppression of the flux of the reflection component at high values of Of. may lead to over- 
estimating the spin, and vice versa. The middle panels of Figs. 13.71 - 13. lOl show the fit 
results for isotropic directionality, which was also the seed model used to generate the 
test data, and so these contours illustrate the magnitude of combined dispersion due to the 
simulated noise and the degeneracy between the spin and the inclination.. The fiducial 
values are well inside the Icr confidence contour in all the graphs in the middle panels. 

However, systematically lower values of the angular momentum are obtained for the 
limb brightening profile and, vice versa, higher values are found for the limb darkening 
profile. The magnitude of the difference is larger for higher values of angular momentum. 



3.2.4 Angular emission profile of the detailed reprocessing model 

The results presented in the previous sections show that using different emission direc- 
tionality approaches leads to a different location of the x^ minimum in the parameter 
space. Doubt about the correct prescription for the emission directionality thus brings 
some non-negligible inaccuracy into the evaluation of the model parameters. The magni- 
tude of this error cannot be easily assessed as a unique number because other parameters 
are also involved. 

A possible way to tackle the problem is to derive the intrinsic spectrum from self- 
consistent numerical computations. This has the potential of removing the uncertainty 
about the emission directionality (although, to a certain degree this uncertainty is only 
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Table 3.2: The best-fit spin and inclination angle values inferred for the three cases of the 
limb darkening/brightening law in the kyl3cr model. 





Case 1 


Case 2 


Case 3 




Ui -- 


= 0.7, Of = 


30° 


a 


60+°°^ 

"•""-0.01 


69-'°*'i 

"•"^-0.01 


0.76:«:«i 


do [deg] 


29 8 +"-2 

^'^•°-0.3 


on 7+0.3 
^^•'-0.3 


29.6!°i 


x^/y 


1.33 


1.27 


1.39 




Uf -- 


= 0.7, Of = 


50° 


a 


65+003 

"•"■^-0.05 


0.73!«:«3 


0-82!«:«^ 


do [deg] 


60.0!«:^ 


60.0!- 


60.3!°;J 


X^/v 


1.87 


1.00 


1.48 




a( = 


0.998, Oi = 


30° 
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956+°*'°^ 

"•'^ -^"-0.005 


1+0 
^-lE-3 


1+" 


0o [deg] 


29.9!|Ji 


29 5+"-^ 


28.7![]:^ 


x^/y 
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5.08 




ai = 


0.998, ^f = 


60° 


a 


982+"°°2 


j+O 


^-6E-5 


Oo [deg] 


59 9+"-i 


6o.r«:i 


60 2+°-^ 

""•^-0.2 


x^/y 


1.24 


1.02 


2.51 



Data were generated using the kyl2cr model. See the main text for details. The quoted 
errors correspond to the 90% confidence level. 
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Figure 3.11: Plots of data/model ratio, where the data are simulated as powerlaw + 
kyl2cr and the model applied to the data is POWERLAW + kyl3cr with a particular 
analytical approach of the directionality. The default angular momentum value is af = 0.7, 
and the emission angle 9j = 30° (left), and 6^ = 60° (right). These parameters and the 
normalisation of the reflection component were allowed to vary during the fitting proce- 
dure. The plotted results correspond to the terminal values of the parameters obtained 
during the;^^^ minimisation process. Other parameters of the model were kept frozen at 
their default values: Y = 1.9, ri„ = Vs, rout = 400, q = 3 and normalisation of the power 
law Kr = 10-2. 
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Figure 3.12: The same plots of data/model ratio, but for a = 0.998 (left: 6{ = 30°, right: 
Of = 60°). 
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Figure 3.13: Results from the test fits with cif = 0.7 and 6{ = 30deg, using the POWER- 
LAW + kyl3cr model applied to the data simulated with powerlaw + kyl2cr. Three 
different profiles of the emission directionality are shown in columns - left, limb bright- 
ening, middle: isotropic, right, limb darkening. Top: dependence of the best fit^^ values 
on the fiducial spin value. The horizontal (dashed) line represents the 90% confidence 
level. Bottom: contour graphs of a versus 6o. The contour lines refer to 1,2, and 3 sigma 
levels. The position of the minimal value of ;^f^ is marked with a small cross. The values 
of ;(f^ corresponding to the minimum and to the contour levels are shown at the top of each 
contour graph. The large cross indicates the position of the fiducial values of the angular 
momentum and the emission angle. Other parameters of the model were kept fixed at 
default values: Y = 1.9, rjn = r^s, ^out = 400, q = 3 and normalisation of the power law 
Kr = 10-2. 



moved to a different level of the underlying model assumptions). In this section, we 
present such results from modelling the artificial data generated by numerical simulations, 
i.e. independently of an analytical approximation of the emission directionality presented 
in the previous sections. 

Let us remind the reader that the orbital speed within the inner < lOrg reaches a 
considerable fraction of the speed of light (Fig. Il.ll ). Beaming, aberration, and the light- 
bending all affect the emitted photons very significantly in this region. Less energetic pho- 
tons come from the outer parts where the motion slows down and the relativistic effects 
are of diminished importance. This reasoning suggests that the analysis of the previous 
section may be inaccurate because the adopted analytical approximations (|3.1I) neglect 
any dependence on energy and distance. 



W e applied the Monte-Carlo radiative transfer code NOAR (see Section 5 in lDumont. Abrassart & Collin 
2000|) for the case of "cold" reflection, i.e. for neutral or weakly-ionised matter. The 



NOAR code computes absorption cross sections in each layer. Free-free absorption and 
the recombination continua of hydrogen- and helium- like ions are taken into account. 
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Figure 3.14: The same as in Figure [3T3l but for a^ = 0.7 and Of = 60deg. 
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Figure 3.15: The same as in Figure [3T3l but for Uf = 0.998 and 6f = 30deg. 



as well as the direct and inverse Compton scattering. The NOAR code enables us to 
obtain the angle-dependent intensity for the reprocessed emission. The cold reflection 
case serves as a reference point that can be later compared against the models involving 
stronger irradiation and higher ionisation of the disc medium. 

The directional distribution of the intrinsic emissivity of the reprocessing model is 
shown in Fig. 13. 21 (right panel). The continuum photon index F = 1.9 is considered and 
the energies are integrated over the 2-100 keV range. Although the results of the radi- 
ation transfer computations do show the limb-brightening eff"ect, it is a rather mild one, 
and not as strong as the Case 1. In the same plot we also show the angular profile of 
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Figure 3.16: The same as in Figure [3T3l but for Uf = 0.998 and 6f = 60deg. 



the emissivity distribution in different energy ranges: (i) 2-6 keV (i.e. below the iron Ka 
line rest energy); (ii) 6-15 keV (including the iron Ka line); (iii) 15-100 keV (including 
the Compton hump). We notice that the energy integrated profile is dominated by the 
contribution from the Compton hump, where much of the emerging flux originates. How- 
ever, all of the energy sub-ranges indicate the limb-brightening effect, albeit with slightly 
different prominence. 

We implemented the numerical results of NOAR modelling of a reflected radiation 
from a cold disc as kyl2cr in the KY collection of models. Furthermore, we produced 
an averaged model, kyl3cr, by integrating kyl2cr tables over all angles. Therefore, 
KYL3cr lacks information about the detailed angular distribution of the intrinsic local 
emission from the disc surface. On the other hand, it has the advantage of increased 
computational speed and the results are adequate if the emission is locally isotropic. Fur- 
thermore, kyl3cr can be a posteriori equipped with an analytical prescription for the 
angular dependence (Cases 1-3 in eq. |3.1| ). which brings the angular resolution back into 
consideration. This approach allows us to switch between the three prescriptions for com- 
parison and rapid evaluation. 

In order to constrain the feasibility of the aforementioned approaches, we generated 
the artificial data using the powerlaw -i- kyl2cr model. The parameters of the model 
are: photon index of the power law F = 1.9 and its normalisation Kr = 0.01, spin of 
the black hole a, inclination angle 6o, the inner and outer radii of the disc rin = r^s 
and Tout = 400, the index of the radial dependence of the emissivity q = 3, and the 
normalisation of the reflection component ^kyi2cr = 0.1. We simulated the data for two 
different values of the spin, a = 0.7 and a = 0.998, and for inclination angles ^o = 30° and 
^o = 60°. The simulated flux of the primary power law component is the same as in the 
previous section. However, now an important fraction of the primary radiation is reflected 
from the disc. The total flux depends on the extension of the disc and its inclination. Its 



3.2. Iron Kg line band examined with different directionalities 67 

value is 4.9 - 5.6 x 10"^' erg cm~^ s"^ for our choice of the parameters. 

As a next step, we replaced kyl2cr by kyl3cr and searched back for the best- 
fit results using the latter model. In this way, using kyl3cr we obtained the values 
of the spin and the inclination angle for different directionalities. The fitting results are 
summarised in Table |3.2l Besides the spin and the inclination angle, only normalisation of 
the reflection component was allowed to vary during the fitting procedure. The remaining 
parameters of the model were kept frozen at their default values. 

The resulting data/model ratios are shown in Figures 13.111 and 13.121 For a = 0.7 
and ^o = 30° the graphs look very similar in all three cases. However, the inferred spin 
value differs from the fiducial value with which the test data were originally created. The 
dependence of the best-fit x^ statistic on the spin and the corresponding graphs of the 
confidence contours for spin versus inclination angle are shown in Figures 13. 13143.161 
again for the three cases of angular directionality. These figures confirm that for the limb- 
brightening profile the inferred spin value comes out somewhat lower than the correct 
value, whereas it is higher if the limb-darkening profile is assumed. 

In each of the three cases the error of the resulting a determination depends on the 
inclination angle and the spin itself. However, we find that the isotropic directionality 
reproduces our data to the best precision. The limb darkening profile is not accurate at 
higher values of the spin, such as a = 0.998, when the resulting ;^^^/v value even exceeds 2. 
The limb darkening profile is characterised by an enhanced blue peak of the line while the 
height of the red peak is reduced (see Figs. 13.31 - 13.41) . Consequently, the model profile is 
too steep to fit the data. This is clearly visible in the data/model ratio plots for a = 0.998 
shown in Figures 13.121 The flux is underestimated by the model below a mean energy 
value ^mean of the line (for a = 0.998 and / = 30 deg ^mean ~ 5 keV) and overestimated 
above ^mean- This fact leads to a noticeable jump in the data/model ratio plot. 

A noteworthy result appears in comparing of the contours produced by the model with 
limb brightening and limb darkening for a = 0.7 and / = 60 deg (Fig. 13.141) . Although 
the former model (limb brightening) gives a statistically worse fit with;^^/y = 1.87 than 
the limb darkening case (x^/v = 1.48), the inferred values of the spin and the inclination 
angle are consistent with the fiducial values within the 3cr level. On the other hand, the 
spin value inferred from the limb darkening model is far from the fiducial (i.e., the correct) 
value. 

It was shown for the case of XMM-Newton data of MCG -6-30-15 that the radial 
emissivity can mimic different angular emissivity (Section 13.2.21 ). Therefore, as a next 
step, we allowed q parameter to vary during the fitting procedure of the artificial data. 
The contour graphs of the spin and the q parameter are shown in Figure 13.171 It is clear 
from the figures that the correct value of the spin can be reached by adjusting of the q 
parameter - for limb brightening by decreasing steepness of radial emissivity, and vice 
versa for limb darkening. However, this is not valid for each set of parameters. For 
a = 0.7, and i = 60 deg, the model with limb darkening overestimates the spin value for 
any value of ^ within 3 cr confidence level. Oppositely, for a = 0.998, and i = 30 deg, the 
model with limb brightening underestimates the spin value for any value of q within 3 cr 
confidence level. 
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Figure 3.17: The contour graphs of the spin a and the parameter q of the radial emissivity 
for different fiducial values of the spin and inclination. From top to bottom: {a = 0.7, 
/ = 30deg}, {a = 0.7, i = 60deg}, {a = 0.998, i = 30deg}, {a = 0.998, / = 60deg}. 
The fiducial value of the radial-emissivity parameter was always q = 3. Three diff'erent 
profiles of the emission directionality are shown in columns - left: limb brightening, 
middle: isotropic, right: limb darkening. The contour lines refer to 1 cr, 2cr, and 3o" 
levels. The position of the minimal value of ;i;^ is marked with a small cross. The values 
of x^ corresponding to the minimum and to the contour levels are shown at the top of 
each contour graph. The large cross indicates the position of the fiducial values of the 
angular momentum and the emission angle. Other parameters of the model were kept 
fixed at default values: Y = 1.9, n^ = rms, ^out = 400 and normalisation of the power law 
Kr = 10-^. 
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3.3 Discussion of the results 



We investigated whether the spin measurements of accreting black holes are affected by 
the uncertainty of the angular emissivity law, MQ^e, E^.), in the relativistically broadened 
iron Ka line models. We employed three different approximations of the angular emissiv- 
ity profiles, representing limb brightening, isotropic and limb darkening emission profiles. 
For the radius-integrated line profile of the disc emission, and especially for higher val- 
ues of the spin, the broadened line has a triangle-like profile. The differences among the 
considered profiles concern mainly the width of the line's red wing. However, the height 
of the individual peaks, also affected by the emission directionality, is important for the 
case of an orbiting spot (or a narrow ring), which produces a characteristic double-horn 
profile. 

We re-analysed an XMM-Newton observation of MCG -6-30-15 to study the emission 
directionality effect on the broad iron line, as measured by current X-ray instruments. We 
showed the graphs of ^^ values as a function of spin for different cases of directionality. 
We can conclude that the limb darkening law favours higher values of spin and/or steeper 
radial dependence of the line emissivity; vice versa for the limb brightening profile. Both 
effects are comprehensible after examining the left panel of Fig. 13.31 The limb darkening 
profile exhibits a deficit of flux in the red wing compared with the limb brightening profile. 
Both higher spin value and steeper radial profile of the intensity can compensate for this 
deficit. 

The higher spin value has the effect of shifting the inferred position of the marginally 
stable orbit (ISCO) closer to the black hole, in accordance with eq. (11.21) . Consequently, 
the accretion disc is extended closer to the black hole. The radiation comes from shorter 
radii and it is affected by the extreme gravitational redshift. Hence, the contribution to 
the red wing of the total disc line profile is enhanced. Naturally, these considerations are 
based on the assumption that the inner edge of the line emitting region coincides with the 
ISCO. Likewise, the steeper radial dependence of the emissivity means that more radia- 
tion comes from the inner parts of the accretion disc than from the outer parts, and this 
produces a similar effect to decreasing the inner edge radius. With the limb brightening 
law the above-given considerations work the other way around. 

We further simulated the data with the povv^erlaw -i- kyrline model. The simple 
model allows us to keep better control over the parameters and to evaluate the differences 
in the spin determination. We used one of the preliminary response matrices for the IXO 
mission and we chose th e flux at a level sim ilar to bright Seyfert 1 galaxies observable by 



current X-ray satellites ( Nandra et a/.ll2007|) . The simulations with the kyrline model 



confirm that the measurements would overestimate the spin for the limb darkening profile 
and, vice versa, they tend to the lower spin values for the limb brightening profile. 

Although the interdependence of the model parameters is essential and it is not pos- 
sible to give the result of our analysis in terms of a single number, one can very roughly 
estimate that the uncertainties in the angular distribution of the disc emission will produce 
an uncertainty of the inferred inner disc radius of about 20% for the high quality data by 
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IXOo We consider this value as realistic. 

In the next step, we applied the NOAR radiation transfer code to achieve a self- 
consistent simulation of the outgoing spectrum without imposing an ad hoc formula for 
the emission angular distribution. We assumed a cold isotropically illuminated disc with 
a constant density atmosphere. We created new models for the KY suite. The results of 
NOAR computations of the cold disc are implemented in the kyl2cr model, while the 
KYL3cr model uses the angle-integrated tables over the entire range of emission direc- 
tions. This enables us to include, a posteriori, the analytical formulae for directionality 
and check how precisely they reproduce the original angle-resolved calculations. We sim- 
ulated the data using the preliminary IXO response matrix with the seed model POWER- 
LAW -I- kyl2cr. Then we replaced the kyl2cr model by kyl3cr to test the analytical 
directionality approaches on these artificial data. 

We found that, for a = 0.998, none of the three assumed cases of the directionality 
profile covered the fiducial values for the spin and the inclination angle within the Scr con- 
tour line (see Figs. 13. 151 - 13.1 71) . The suitability of the particular directionality prescription 
depends on the fiducial values of the spin and the inclination angle. The limb brightening 
profile successfully minimises the;if^ values for a = 0.998 and i = 30deg (though with 
slightly different parameter values), but for a = 0.7 and / = 60 deg it gives the worst fit of 
all studied cases of the directionality laws. 

On the whole, we found that the isotropic angular dependence of the emission inten- 
sity fits best. The model with the limb darkening profile was not able to reproduce the 
data, especially for higher values of the black hole spin. The best fit;^'^/v value exceeds 5 
for a = 0.998 and i = 30 deg, which means a more than three times worse fit than using 
isotropic or limb brightening directionality. The inclination angle was underestimated by 
more than 1 deg, the ^-parameter was overestimated by almost 10%. 

This is an important result because much of the recent work on the iron lines, both 
in AGN and black hole binaries, has revealed a significant relativistic broadening near 
rapidly rotating central black holes (see e.g. Millen 12007 '). In some of these works. 



the limb darkening law was employed and different options were not tested. The mod- 
elled broad lines are typically characterised by a steep power law in the radial part of 
the intensity across the inner region of the accretion disc, as in the mentioned MCG -6- 
30-15 observation. This behaviour has been interpreted as a case of a highly spinning 
compact source wh ere the black hole rotational energy is electromagnetically extracted 



(|Wilms et al.\\200l\) . We conclude that the significant steepness of the radial part of the 
intensity also persists in our analysis, however, the exact values depend partly on the 
assumed angular distribution of the emissivity of the reflected radiation. 

It should be noted that, in reality, the angular distribution of the disc emission is sig- 
nificantly influenced by the vertical structure of the accretion disc. However, our compre- 
hension of accretion disc physics is still evolving. In recent years, several detailed models 
have been developed for irradiated black hole accretion discs in hydrostatic equilibrium 



^Here, we refer to the inner disc radius instead of the spin because the latter is not as uniform a quantity 
as the corresponding ISCO radius to express the uncertainty simply as a percentage value. However, we 
still suppose that the inner edge of the line emitting region of the disc coincides with the marginally stable 
orbit. 
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(see e.g. lNavakshinll200d : lBallantvne. Ross & Fabianll200ll : lR6zanska et Q/.ll2002h . These 
models combine radiative transfer simulations with calculations of the hydrostatic balance 
in the stratified disc medium. Aside from reprocessing spectra, the models provide solu- 
tions for the vert i cal di sc profile of the density, temperature, and ionisation fractions. In 
Goosmann et al.\ (|2007|) the effects of general relativity and advection on the disc medium 
were added. 



In iNayakshinI (l200nh . iGoosmann et aU (l2007h . and iRozahska & Madejl (l2008h the re- 
processed spectra are evaluated at different local emission angles. The shape and normal- 
isation of these spectra depend on various model assumptions. Until we know in more 
detail how accretion discs work, it is hard to choose which is the "correct" reprocess- 
ing model. One could always argue that more physical processes and properties should 
be included in the radiative transfer simulations, such as the impact of magnetic fields, 
macroscopic turbulence, or different chemical compositions of the medium. Using the 
above-mentioned models for an accretion disc in hydrostatic equilibrium may then easily 
become computationally intense. 

For the practical purpose of data analysis, however, the computation of large model 
grids is necessary, which requires sufficiently fast methods. This is why simple con- 
stant density models are most o ften used to analyse the observational data. In fact, 
Ballantyne. Ross & FabianI (1200 1|) have shown that their reprocessing spectra for a strat- 
ified disc medium in hydrostatic equilibrium can be satisfactorily represented by spectra 
that are computed for irradiated constant density slabs. Therefore, we included in our 
present analysis the angular emissivity obtained from the modelling of neutral reprocess- 
ing in a constant density slab. We have not discussed in any further detail the dependence 
on the ionisation parameter, which we expect to be rather important, and it will be possibly 
addressed in a future work. 

We emphasise that the main strategy of our presented research is not intended to find 
the "correct" angular emissivity, as this is still beyond our computational abilities and un- 
derstanding of all the physical processes shaping the accretion flow. Instead, we examined 
the three different prototypical dependencies which are mutually disparate (i.e., the limb- 
darkening, isotropic, and limb-brightening cases), applied in current data analysis, and 
which presumably reflect the range of possibilities. By including these different cases we 
mimic various uncertainties, such as those in the vertical stratification, and we estimate 
the expected error that these uncertainties can produce in the spin determinations. Fur- 
ther detailed computations of reprocessing models and the angular emissivity are needed 
in the future in order to understand its role in different spectral states of accreting black 
holes. 



Chapter 4 

Data reduction and spectral analysis 



4.1 Preliminaries 

4.1.1 XMM-Newton satellite 

X-ray astronomy is essential to study neutron stars and accreting black holes. These 
objects would not be revealed on the sky without X-ray detectors which were first used 
in the early 1960s. It was quite a surprise at that time that the X-ray sky was so different 
from the optical sky. The first detected X-ray source besides our Sun, Scorpius X-1 
(discovered by Aerobee 150 rocket launched on 12th June 1962), was found to be 10,000 
times brighter in X-rays than in its optical emission. X-ray astronomy is not possible from 
the Earth surface because our atmosphere blocks out all X-rays. Only telescopes above 
the atmosphere can detect X-ray radiation. Following the first "rockoons" experiments 
(hybrid of a rocket and a balloon) many X-ray satellites have been launched. 

The XMM-Newton satellite (X-ray Multi Mirror) belongs to the largest scientific 
satellites ever launched and it represe nts a cornerstone m ission of the European Space 



Agency's Horizon 2000 programme ( J ansen et al.\ |2001|) . It was launched by the Ari 



ane 504 rocket on 10th December 1999, and it is planned to operate still several more 
years. Its large effective area (4700 cm^) is its major advantage compared to other recent 
satellites (like Chandra). The satellite is in a highly eccentric orbit (between 7000 and 
114,000km) allowing long uninterrupted exposure times (^ 40 hours). Observing time 
on XMM-Newton is being made available to the scientific community, applying for ob- 
servational periods on a competitive basis. The data are stored in the publicly available 
archivqll. 

The XMM-Newton satellite carries three CCD cameras for X-ray spectroscopy and 
imaging - Europea n Photon Imaging Ca niera, EPIC, which includes two MOS detectors 



(|Tumer et q/.II2001|) and one PN detector (IStriider et a/.ll2001|) . two spect rometers for high 



resolu tion X-ray spectroscopy - Reflection Grating Spectrometer, RGS Aden Herder et al. 



20011). and one o ptical/UV imaging and grism spectroscopy instrument - Optical Monitor, 



OM (IMason et al. 2001). All these instrument are able to operate simultaneously. The 



'http://xmm.esac.esa.int/xsa 
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Table 4.1: Capabilities of detectors on-board XMM-Newton satellite. 



Instrument 


EPIC PN 


EPIC MOS 


RGS 


OM 


Energy bandpass 


0.15- 12keV 


0.15-12keV 


0.35- 2.5 keV 


180- 600 nm 


Orbital target visibility ' 


5-135ks 


5-135ks 


5-135ks 


5 - 145 ks 


Sensitivity ^ 


« 10-14 


« 10-14 


~ 8 X 10-5 


20.7 mag 


Field of view 


30' 


30' 


«5' 


17' 


Spectral resolution ^ 


~80eV 


«70eV 


0.4/0.025 A 


350 



Notes: i The orbital target visibility represents the total time available for scientific measurement 

per orbit. The value depends on what time the satellite is outside van Allen radiation belts which 

could damage the detectors. 

^ The sensitivity is a) for EPIC detectors: in units of erg s-i cm-^ measured in 0. 15 - 15 keV range 

after 10 ks exposure, b) for RGS detectors: in units of counts s~i cm"^ measured in O VII 0.57 keV 

line flux with a background of 10-4 counts s-i cm-^, c) for OM: 5 cr detection of an AO star in 1 ks 

exposure. 

^ The spectral resolution of EPIC detectors is at 1 keV (at the energy of Fe Ktc line E = 6.4 keV 

the energy resolution of both detectors is « 150 eV. The resolution of OM is the resolving power 

(A/AA). 

main properties of the instruments are summarised in Table 14. ll 

Observation modes of EPIC detectors 

The EPIC detectors are suitable for our research of broad iron lines with their energy 
range (see Table 14. ll) and high effective area. There are several observation modes which 
differ from each other mainly in the covered field of view and the read-out time. The basic 
characteristics of the PN and MOS observation modes are summarised in Tables 143114. 31 
The "Full Frame" or "Extended Full Frame" modes are useful for faint sources where 
long read-out time is necessary to reach a reasonable signal to noise ratio. Oppositely, 
"Timing" mode or "Burst" mode are dedicated for observation of very bright sources to 
avoid the problem with pile-up (see Section 14.1. 21) . The duty cycle in the "Burst" mode is 
only 3% which means that the rest 97% of counts are lost. 

SAS - tool for reduction of the XMM-Newton data 



The Science Analysis System (SAS) is a collection of tasks, scripts and libraries, specif- 
ically designed to r educe and analyse data collected by the XMM-Newton observatory 
(|Gabriel et a/.ll2004|) . It is publicly available at http://xmm.esac.esa.int/sas, where also the 
instructions for install and usage are described. In this Section, we briefly summarise 
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Table 4.2: EPIC PN observation modes. 



Modes 


Time resolution [ms] 


Live time [%] 


Max. count rate* [cts s" 


-'] 


Extended Full Frame 


200 


100 


2 




Full Frame 


73.4 


99.9 


6 




Large Window 


48 


94.9 


10 




Small Window 


6 


71 


100 




Timing 


0.03 


99.5 


800 




Burst 


0.007 


3 


60000 





Table 4.3: EPIC MOS observation modes. 



Modes 


Time resolution [ms] 


Live time [%] 


Max. count rate* [ctss ^] 


Full Frame 


2600 


100 


0.7 


Large Window 


900 


99.5 


1.8 


Small Window 


300 


97.5 


5 


Timing 


1.75 


100 


100 



*The maximal count rate values are shown for a point source. Only observations of sources with 
lower flux avoid problems with an irremovable pile-up. 

what is needed in preparation of the X-ray data before the spectral analysis and which is 
all accessible with the SAS software: 



1. The raw XMM-Newton data, called Observation Data Files (ODE), are provided 
in EITS (Elexible Image Transport System) format. The SAS software allows to 
create a calibrated and concatenated event list where events from all detector chips 
are included in one single file. The information about the instrumental calibration 
is necessary to properly interpret the data. The researchers are encouraged to use 
the most updated calibration data set. 

2. The SAS software helps its users to create scientific products according to their 
own taste and data screening criteria. High X-ray background flare events are stan- 
dardly removed by creating of good-time intervals of the observation where the 
background level is under a certain criterion. 
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3. After choosing good-time intervals, an X-ray image might be created. The source 
spectrum is extracted from the region on the chip where the point spread function 
(PSF) is the brightest (it is usually a circle around the centre of the PSF). The 
background spectrum is extracted from the outer regions of the same chip where 
the PSF is weaker. 

4. Each X-ray instrument has a typical response on the incoming radiation and has 
a different sensitivity in different wavelengths. The instrumental spectral transfer 
functions (response matrices) are needed for the quantitative spectral analysis. The 
SAS software allows the users to generate it for the individual observation. 

5. The SAS software contains also a tool for estimation of the pile-up fraction in the 
data (see Section |4. 1.21) . and many other useful tasks (for more information see 



Gabriel et al.\ (|2004|) or look at the SAS web pages). 



4.1.2 Pile-up 

Pile-up is one of the biggest problems of the detection technique using a CCD camera 
when the observed source is too bright. Any CCD camera is composed from individual 
pixels which are collecting photons. The energy of an incoming photon invokes a charge 
in the photo-active region (in typical CCDs, epitaxial layers of silicon are used). The 
charge is then read-out during the duty cycle. Problem with pile-up occurs when two 
or even more photons deposit charge packets in a single pixel ("photon pile-up"), or in 
neighbouring pixels ("pattern pile-up") during one read-out cycle. 

In general, the pile-up makes spectra harder because two or more soft photons are de- 
tected as only one hard photon. In addition, when the summed energy is above a rejection 
threshold, it leads to a more or less pronounced depression of counts in the central part 
of the point spread function implying undesirable flux loss. The Tables |42ll4.3l contain 
the values of the maximal count rates of the observed source fluxes for different obser- 
vation modes to avoid the problem with pile-up. In general, the MOS detector is more 
susceptible to pile-up than the PN camera. 

If the pile-up is only moderate, it might be removed by: 

1. using only single events, i.e. events where no signal is detected in neighbouring 
pixels. The single events are less sensitive to the pile-up than any multiple events. 

2. excising of the PSF core. However, this procedure may cause a drastic loss of 
counts. 

A SAS task epatplot is a suitable tool for constraining the pile-up level. It measures 
the "pattern pile-up" which is proportional to the "photon pile-up". The pattern pile-up 
occurs when two single events are interpreted as one double event with twice the energy. 
The probability of single and double events fractions can be modelled for a given flux. If 
there is a pile-up the double events are more frequent than predicted. The rate of diversion 
of the data from the theoretical model curves corresponds to the pile-up level. 
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We encountered the pile-up problem in the observation of a stellar-mass black hole 
binary, GX 339-4, performed by the XMM-Newton satellite in 2004, see Section |4Xn for 
details. 



4.1.3 Re-binning of the data 

There is not a unique way of the treatment with the data after their reduction described 
in Section 14. 1.11 Some astronomers use the data for the spectral analysis without any 
re-binning. Some re-bin the data but with different conditions. The common condition 
for grouping is to have at least a certain number of counts per bin, e.g. 20. The reason is 
that the commonly used statistics for fitting the data is the Gaussian one with;^'^ values 
describing the fit goodness, and the derived formulae of the Gaussian statistics are based 
on the assumption of a sufficient number of counts (see Section|4.1.4| ). However, when the 



number of counts per bin is not adequate C-statistics (|Casblll976|) may be used, instead. 
Further, we discuss re-binning of the data because of the energy resolution rather than an 
insufficient number of counts. 

The counts are intrinsically distributed in many channels of the instrument according 
to the detected energy of photons. However, the energy resolution of currently used X-ray 
instruments is not so brilliant and it is largely exceeded by the total number of channels. 
For instance, the PN detector has the energy resolution of order of « lOOeV (see Ta- 



ble [471]) and the energy range « lOkeV which is about « lOOx larger than the energy 
resolution. Thus, the corresponding number of energy bins of the width equal to the order 
of the energy resolution is ^ 100. However, the PN detector has 4095 intrinsic channels 
in total, i.e. more than one order of magnitude more. 

Hence, we find necessary to re-bin (group) the data with respect to the intrinsic res- 
olution of the instrument. We find sufficient re-binning which over-samples the energy 
resolution by a factor of three. We use the pharbn script provided us by M. Guainazzi 
for this purpose. The energy resolutions of the instruments are typically approximated as 
power laws of the energy (with an exponent e = -0.5 for the PN detector and e = -0.46 
for the MOS detector), and with a known energy resolution at some reference energy 
(2.26% at 6keV for PN, 6.31% at 1 keV for MOS). 

We emphasise that the interpretation of spectral results can be influenced by the 
adopted procedure of the data re-binning, and so different authors may reach different 
conclusions. We will illustrate this fact on an X-ray spectrum of GX 339-4 (see the de- 
tails in Section |?.2.2D . 

4.1.4 Goodness of the fit 

In this section, we br i efly d escribe how to constrain the goodness of the fit. Readers are 



referred to e.g. iPressI (|2002|) for more information about different statistical techniques to 
analyse the data. 

X-ray spectra are measured with spectrometers which do not obtain the actual spec- 
trum but rather photon counts C within specific instrument channels /. This observed 
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spectrum is related to the actual spectrum of the source f{E) by: 

Jr»cx> 
f{E)R{I,E)dE, (4.1) 



where R{I, E) is the instrumental response which is proportional to the probability that an 
incoming photon of energy E will be detected in channel /. 

The instrume ntal spectral transfer fun ctions (response matrices) are not generally in- 



vertible (see e.g. iBlissett & Cruiselll979r) . and it is therefore insecure to reconstruct the 
original spectrum (so called unfolded spectrum). Instead of deriving f{E) for a given 
set of C(/), model spectra fmiE) are calculated how they would look like when passing 
through the detector. The predicted counts Cp{I) are: 

fM(E)R(I,E)dE, (4.2) 



and then, they are compared with the data counts C(/). 

The quality of the fit, i.e. how the model is able to interpret the data, is usually 
characterised by the;^^ value given by: 

where C(I) is the number of counts in the energy bin / (DATA), Cp(I) is the number of 
predicted counts in the energy bin / (MODEL), and cr(I) is the (generally unknown) error 
for the channel /. The error cr(/) is usually estimated as « y/C(I). 

The x^ statistics provides a goodness-of-fit criterion for a given number of degrees 
of freedom v, which is calculated as the number of channels minus the number of model 
parameters, and for a given confidence level. If ;^^^ value exceeds a critical value one can 
conclude that the fit /m(£) is not an adequate model for C(/). As a general rule, one 
wants to have a reduced ;^f^gj = x^/y to be approximately equal to one (x^ ~ v). A reduced 
X^ value that is much greater than one, x^^d ^^ ^' iiidicates a poor fit, while ;^f^g^ << 1 
indicates that the errors on the data have been over-estimated. 

^^^red ~ 1' so the best-fit model /«(£') pass the "goodness-of-fit" test, one still cannot 
say that /m(£^) is the only acceptable model. For example, if the data used in the fit are not 
particularly good, one may be able to find many different models for which adequate fits 
can be found. In such a case, the choice of the correct model to fit is a matter of scientific 
judgement. 



There are several software packages to analyse X -ray s pectra - XSPEC (lArnaiid 



1996i) . SPEX (iKaastra. Mewe & Nieuwenhuiizenlll996h . ISIS (iHouck & Denicolall2000D . 



For all the spectral analysis presented in this Thesis, we used the XSPEC fitting package. 
The best-fit values of the model param eters were found using the modified Levenberg- 



Marquardt algorithm (iBevingtonI 119691) . which is the default one in the XSPEC. This 
algorithm is local rather than a global one, so we were aware of the possibility to get 
stuck in a local minimum instead of the global minimum by starting of the fitting process 
from reasonable initial values of model parameters. Therefore, we standardly repeated 
fitting procedures for several different initial parameter values. 
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4.2 GX 339-4 



GX 339-4 belongs to the Galactic black hole candidates (JRemillard & McClintoc5l2006|) . 
Dynamical properties of the system are still not completely known. The companion star 
stays unresolved by the present measurements even during the quiescent state of the ob- 
ject. However, the optical sp ectroscopy reveals periodic behaviour of He I I and Bowen 



blend N III lines in spectra ( Hynes et al.\ l2003l : iBuxton & VennesI l2003h . iHynes et al. 
(|2003[) determined the orbital period of the system as P « 1 .Id. From the given period 
the radial velocity curve K2 of the secondary can be fitted. The mass function /(M) is 
given by eq. (fT35]) . For GX 339-4: f{M) = 5.8 + 0.5. The mass of the black hole can be 
determined from the knowledge of the inclination of the system and the ratio of the mass 
of the companion to the mass of the black hole (see eq. 11.55^ 



However, these quantities are not yet well constrained. Hynes et al.\ (l2003h estimated 
the mass ratio q < 0.08 from a small modulation of the wings of the He II line. The incli- 
nation angle is unknown. However, a jet is resolved from radi o and infrared obse rvations. 
The derived inclination from these measurements is i < 26° (|Gallo et a/.ll2004|) . but this 
need not be necessarily related with the inclination of the orbital plane. 

We chose this object for our analysis because it exhibite d a relativistic b roadened iron 
line p rofile in the X-ray spectra in both high/soft state (.Miller et Q/.ll2004h and low/hard 
state (|Miller et al.luOOw . see Figure l24l This fact rather contrasts with the interpretation 
of the low/hard state via a truncated disc (Section [T3]2]). The extremely skewed line pro- 
files suggest that the radiation comes from the innermost parts of the accretion disc within 
a few gravitational radii. Assuming accreting material radiates only above the marginal 



stable orbit it implies a high value of the spin of the black hole (Section O. iMiller et al 
(|2008|) determined the value of the spin a = 0.93 + 0.01 (statistical) ±0.03 (systematic) by 



joint fitting of two XM M-Newton and on e Suzaku observations. The XMM-Newton spec- 
tra were re-analysed by lReis et a/.l(|2008|) who measured the spin value as a = 0.935+0.01 
(statistical) ±0.02 (systematic). 



4.2.1 Low/hard state observation 

The low/hard state observations of GX 339-4 were performed by the XMM-Newton satel- 
lite in March 2004 (Obs. ID #0204730201 and #0204730301). The MOS detectors were 
operating in the "Full Frame" mode. Although the object was in the low/hard state, it was 
still bright with the count rate > 20 cts s~\ i.e. more than one order of magnitude greater 
than the maximal recommended value to avoid the pile-up problem (0.7cts s ^ see Ta- 
ble |4]3l). Indeed, we found a substantial level of the pile-up in the data using the epatplot 
tool (see Section |4T2l) . 

Figure 14.11 shows the central chip of the MOS 1 detector of the first observation (it 
is similar in the second observation). The black core itself indicates the over-exposition 
of the detector. Figures 14.21 and 14.31 show pattern distribution with the energy for two 
distinctly extracted regions from the central MOS chip of the two observations. The first 
extracted region is a circle around the core of point spread function (PSF) which is the 
standard way of the region extraction to obtain the source spectra if there is no pile-up 
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Figure 4.1: Central chip of the MOS 1 detector (zoomed on PSF) of the observation 
ID #0204730201. Radii of the green circles are 18, 80 and 120arcsec. 



problem. The radius was chosen as 120arcsec. The second region is an outer annulus 
between 80 and 120arcsec. In this case, however, the most of the flux is lost, and the 
signal to noise ratio becomes very low. 

The more deviated from the expected curves the pattern fractions are (bottom panels 
of the figures) the more significant the effect of pile-up is. It is clear that the pile-up 
is highly presented when the extracted region is the whole circle including the core of 
PSF. However, the pile-up is not completely removed nor for a faraway outer annulus 
though the level of the pile-up is much lower. Hence, the spectral analysis might be 
possible in combination with the use of only single events. However, the drastic excision 
of such a large region causes that the signal to noise is so low that the quantitative spectral 
analysis gives very poor constraints about the spectrum and is not suitable for a detailed 
investigation of any model parameters (not speaking of a very sensitive spin parameter). 

Besides the two extreme extraction regions, we checked several annuli with the inner 
radius between a nd 80 arcsec. We found that the pile-up is significant up to » 80 arcsec. 
Miller et al\ (120061) were satisfied with an annulus 18 - 120 arcsec. However, Figure BT^ 
shows very strong indications of the pile-up in both observations for this extraction region. 
Because of the pile-up problem, we left these data sets as unsuitable for our purposes, 
regarding any spectral results obtained with t hese data as untrustworthy. 

Meanwhile, the data were re-analysed by Done & Diaz Trigol (|2009() who also found 
the pile-up level significantly high. Moreover, they performed a more comprehensive 
analysis including FN data, as well. They showed that the FN data are less affected by 
the pile-up than MOS data, and so they can be used for the spectral analysis if only single 
events are considered and the brightest core of FSF is excised (in the FN "Timing mode", 
it means to excise several central columns of the source image). They concluded that the 
FN data do not support the presence of a broa d ir on line as sugges ted from the piled- 
up MOS data (presented by JMiller g? g/.l l2006l and iReis grQ/]|2008b . They argued that 
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Figure 4.2: Quantitative pile-up diagnostics in the MOS 1 camera of the observation 
ID #0204730201 using the SAS task epatplot. The data points correspond to single events 
s (red), double events d (dark blue), triple events t (green), and quadruple events q (light 
blue). Up: Spectra related to the individual patterns. Bottom: Pattern fractions. The solid 
lines represent expected fractions if there is no pile-up effect. Left: Data extracted from a 
circle with radius of 120arcsec. Right: Data extracted from an outer annulus between 80 
and 120arcsec. 



missing relativistic smearing supports the truncated disc int erpretation of t he low /hard 
state. Similar issues with the pile-up were also reported by lYamada et all (|2009h with 
Suzaku data of GX 339-4. 



Evidently, the pile-up issues and the possible ways of removing its degradation of 
the spectra are subject of discussions. Different researchers seem to agree on the general 
conclusion that the pile-up can be a serious problem. However, in individual observations 
the impact of the effect is often estimated differently. It appears that the problem can be 
definitively avoided only by future high throughput detectors. 
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Figure 4.3: The same quantitative pile-up diagnostics as in Figure |4~2l but for the sec- 
ond observation (ID #0204730301). Left: Data extracted from a circle with radius of 
120 arcsec. Right: Data extracted from an outer annulus between 80 and 120 arcsec. 



4.2.2 Very high state observation 

We also re-analysed the observation of GX 339-4 in the very high state (VHS). The PN 
detector was operating in the "Burst mode" which avoids the pile-up problem (see Ta- 
ble |4]2l). Using the SAS task epatplot we found the pile-up to be insignificant with these 
data. Therefore, we used these data for the further spectral analysis and also for the 
comparison of the relativistic line models in Section 12.3. 1[ However, we obt ained a n 



unacceptable fit with A'reduced ~ 6.9 when we applied the same model as in Miller et al. 



(12004) to the data in the energy range 0.7 - 9 keV. The difference of the results is likely 
due to the different re-binning of the instrumental energy channels (see Section |43]). 



Fig ure l43]shows ratios o f differently binned data to the same model as used in lMiller et al 



(|2004|) and lReis et al.\ (|2008|) . The upper panel shows the data re-binned with the only con- 
dition to have at least 20 counts per bin, which was used in the mentioned works. This con- 
dition is, however, very weak with respect to the total number of counts Accounts ~ 1 .0 x 10' 
and the total number of energy channels A^chan = 1 .5 x 10^ in the 2 - 10 keV energy range. 
As a result, the data stayed practically un-binned and the energy resolution was exces- 
sively over- sampled. 

In the lower panel of Figure |431 the data are re-binned with respect to the instrumental 
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Figure 4.4: The same quantitative pile-up diagnostics as in Figure |43] and Figure l43] but 
for the same extracted region of the source as in Miller et aL (,2006,) . i.e. an 1 8 - 1 20 arcsec 
annulus. Left: Observation #0204730201. Right: Observation #0204730301. 



energy resolution (the energy resolution is over-sampled by a factor of three). An absorp- 
tion feature at high energies which was concealed in large error bars of the ungrouped 
data is now clearly revealed. The discrepancy is even more visible in the quantitative 
spectral analysis, i.e. after constraining the fit-goodness in the 0.7 - 9 keV energy range 
by the reduced ;t^ values (see Section |4T4l). While ;t'^^j = x^/v - 2963/1659 = 1.8 
for the ungrouped data x^^^ - 1427/208 = 6.86 for the grouped data. After adding 
a Gaussian line to model a sp ectral feature at 2.3 keV, possibly being some calibration 
issue only (iMiller gra/.ll2004 . Afred - 2453/1657 = 1.48 for the ungrouped data, and 



^reduced ~ 921/206 = 4.47 for the grouped data in the energy range 0.7 - 9keV. So, the 
model would still have not been acceptable if applied to the appropriately grouped data. 

Further, we analysed the spectrum of the r e-binne d data. First, we removed an arti- 
ficial SMEDGE model which was used by .Miller et al\ (i2004.) for fitting a smeared edge. 
As a next step, we allowed the model parameters of the continuum and also the iron line 
to float. As a result, we found a better fit with;tVv = 277/202 = 1.7 in 0.7-9 keV. The 



model p arameters are shown in Table 14.41 and the y are compared there with the values 
found bv lMillerg?aZl(l2004h and iReis grail (l2008b . 



4.2. GX 339-4 



83 



o 

"a 
o 

E 

■(5 

T3 




Energy [keV] 



Fig ure 4.5: R atios of the PN data of GX 339-4 to the model used in iMiller g? g/.l (|2004() 
and lReis et al\ (|2008|) with ungrouped data {top) and grouped data {bottom) as described 
in the main text. 
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Figure 4.6: Left: The X-ray spectrum of GX 339-4. The residuals from the model are 
plotted in the bottom panel of the graph (see the main text for further details). Right: 
More detailed view of the iron line band. 



Our fit is characterised by a steeper power law with F = 3.08 + 0.05, a slightly hotter 
disc with fcTin = (0.87 + 0.01) keV, and a less prominent broad iron line. The column den- 



sity hh = (0.63 + 0.01) X 10^^ cm"2 is slightly higher than the value m = .374 x 10^^ cm"^ 



from the Leiden/Argentine/Bonn H I measurements (|Kalberla et q/.II2005|) . which may be 
explained by an additional local absorption. We tried also to add a reflection component 
into the model, as pexrav (pexriv) or refsch model, but without any improvements 
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Table 4.4: Results for PN spectrum of GX 339-4 (VHS) in 0.7-9 keV. 



model 


parameter 


Miller et al. (2004) 


Reis et al. (2008) 


our fit 


PHABS 


nullO^^cm-^] 


0.51+0.01 


0.53 


0.63 + 0.01 


POWERLAW 


r 

normalisation 


2.60-0.05 

9 9+0.3 
^■^-0.1 


2.553 + 0.003 

9 C9+0.01 
^•-^^-0.02 


3.08 + 0.05 

5 64+°- *2 
-'•"^-o.is 


DISKBB 


kT [keV] 
normalisation 


0.76 + 0.01 

2300!^°° 


726+°*'°* 

"•'^"-0.001 

2545 + 20 


0.87 + 0.01 

1450!^! 


LAOR 


E [keV] 

^[deg] 

1 
normalisation 


6.97_o.2 

2.1!- 

11-1 
5 5+°-^ 

■^■■^-0.1 

077+0005 

"•"' ' -0.015 


f. OO+0.09 
"•^^-0.008 

1 Q 1+0.02 
^•^^-0.005 

18.8 + 0.5 

f. OC+0.09 
"•^^-0.008 

0.113+0.002 


6.97-0.2 

-'•"-o.e 
17 + 4 

3.3 + 0.1 

006+°°°'* 

"•"""-0.002 


SMEDGE 


E [keV] 

T 

W [keV] 


7 9+0.1 

"•"-0.1 

1.0 + 0.3 


7.1+0.1 

2 5+°-2 

^•-^-0.1 

4 5+°-5 

^•-■-2.0 


- 




x'/v 


3456.5/1894 


2344.6/1652 


277/198 



Note: iMiller et al\ (12004) required the photon index of the power law to lie within AF < 0.1 of the 
value from the simultaneous RXTE measurements. The value F = 2.6 is their upper limit. 



of the fit. i Reis et al\ (|2008|) used the refhidden model for fitting of the reflection com- 
ponent from the hot disc. However, we did not have the opportunity to use it because it is 
not a publicly available model. 



The spectrum of GX 339-4 is shown in Figure 14.61 with a detailed view of the iron 
line band in the right panel. A broadened iron line feature is still present. However, due 
to the differ ent slope of the und erlying power law the line is muc h weaker than the o ne 
presented in iMiller et aM (|2004l) . The equivalent width found by ' MiUer et al\ (12004) is 



200^40 eV while our value is 154^^^ eV. The main difference is in the determination of the 



inner disc radius inferred from the iron line modelling. While Miller et al.\ (|2004r) and 
Reis et al\ (|2008h concluded that the black hole in GX 339-4 is very rapidly rotating with 
a ~ 0.9 we found an intermediate value for the spin from our fitting, a ~ 0.7 (see also 
Section [23J]). 
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4.3 MCG -6-30-15 



MCG -6-30-15 is a nearby Seyfert galaxy (z = 0.00775) of the type I (see Section [jAll)- 
This object is one of the most studied objects in X-ray doma in since it exhibited a very 
broad iron line in the spectrum of the ASCA observation (JTanaka et all Il995h . The 
broadness of the line has been attributed to its origin from the relativistic accretion disc 
around a rapidly spinning black hole. Further investigations of this object with fol- 
lowing X-ray mis sions confirmed the p resence of the relativistically broade ned iron Kg 
line: BeppoSAX (lOuainazzi et a/.lll999l), XM M-Newton (.Wilms gf al. 2001; Fabian e t al. 



2002 



Martocchia . Matt & Karas 2002; Vaug han & Fabianll2004l:lB renneman & Reynolds 



2006h . Chandra + RXTE (JLee g?Q/.ll2002h . and Suzaku (JMiniutti g?Q/.ll2007b . All cited 

authors concluded that the extension of the red wing can be reasonably interpreted only 

with mo dels employing a very high value of the black hole spin. However. i Miller. Turner & Reeves 
(|2008bh showed in their analysis that the substantial X-ray data set for MCG -6-30-15, 
comprising all the mentioned missions, could be equally well fitted by an absorption- 
dominated model with no relativistically smeared emission. The main idea of this work 
is that the source is obscured by several clumpy absorption zones ("partial covering" ab- 
sorption). 

A possible way how to distinguish between the two different models is to study the 
time varia bility. The lack of variabi lity in the energy range of the broad iron line was 
noticed by Fabian & VaughanI (l2003h in the XMM -Newton spectrum and explained by 
the light bending model (|Miniutti & FabianI 12004) . They constrained the inner edge 
of the disc to be wi thin two gravitational radii (the spin value a > 0.94). However, 
Niedzwiecki & Zycki, (.2008.) pointed out that the line profile has always a pronounced 
blue peak in the light be nding model, wherea s the blu e peak is missing in the deep min- 
imum state observation dlwasawa et al.\\l99w . Hence, iNiedzwiecki & Zyckil (|2008|) con- 
cluded that the innermost region of the disc is not closer than 2-3 gravitational radii. 
It has been also shown that a possible warm absorber may have an important effect on 
the spectral variability ( "Miller. Turner & Reevesll2008bl : iMiyakawa et al.,, 200% . Within 
these models, the variability is explained by changing of the ionisation state of the warm 
absorber. 

The uncertainty which model is more realistic pertains to all similar sour ces. A strong 



argument against partial-covering warm absorber model was pronounced by Fabian et al. 



([2009). In the X-ray spectrum of an extragalactic source, 1H0707-495, the authors de- 
tected an iron L line as well (thanks to extraordinarily high abundances of iron), and 
measured a time reverberation lag between the direct X-ray continuum and its reflection 
from matter falling into black hole. The time lag of 30 s is comparable with the light 
travel time within one gravitational radius, implying that the radiation must come from 
the closest neighbourhood around a maximally rotating black hole. 



Long XMM-Newton observation of MCG -6-30-15 



A long XMM-Newton observation took place in summer 2001 (31st July - 5th August) 
with the acquired exposure time was about 350ks (Obs. #0029740101, #0029740701, 
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#0029740801). The EPN and both MOS cameras were operating in the "Sm all Window" 



mode (see Section 14.1.11). The spectra were analysed by several authors jFahian et al. 



2OO2I : iBallantvne. V aughan & Fabian'20031: IVaughan & Fabianll2004l : iBrenneman & Reynolds 
20061 : iMiller. Turner & Reeves. ,2008b.. etc.). We re-analysed the data and then, we used 
them for analyses presented in the previous sections of the Thesis (Section 123) and [3!Z2]) . 
In this section, we describe the details of our re-analysis. 

We reduced the EPN data from three sequential revolutions (301, 302, 303) using 
the SAS software version 7.1.2. Intervals of high particle background were removed by 
applying count rate thresholds on the field-of-view (EPIC, single events) light curves of 
0.35cts/s for the PN, and 0.5cts/s for the MOS. The patterns 0-12 were used for both 
MOS cameras, and patterns 0-4 (i.e. single and double events) for the PN camera. The 
source spectra were extracted from a circle of 40 arcsec in radius defined around the 
centroid position with the background taken from an offset position close to the source. 
The two MOS spectra and the related response files were joined into a single spectrum 
and response matrix. 

Using the FTOOL MATHPHA (part of the HEASOFT) we merged together the three 
spectra to improve the statistical significance. Further, we used the XSPEC version 12.2 
for the spectral analysis. We re-binned all the data channels to over-sample the instru- 
mental energy resolution maximally by a factor of 3 and to have at least 20 counts per 
bin. The first condition is much stronger with respect to the total flux of the source - 
4 X 10"'' ergcm"^ s~' (1.1 x lO^cts) in the 2-lOkeV energy interval. 



X-ray continuum 



The resulting X-ray spectrum of MCG -6-30-15 is shown in Figure 14771 Above E x 
2.5 keV, the X-ray continuum can be described by a power law component with the photon 
index F = 1 .9 absorbed by Galactic gas matter along the line of sigh t with colurnn den sity 
nu = 0.41 X 10^' cm"^ revealing a broad iron line at £ « 3 - 7 keV (JFabian et a/.ll2002h . It 
is clear from the figure that other components need to be added into the model to describe 
the spectrum at lower energies as well. 

The shape of the data residuals in lower energies suggests a combination of two ef- 
fects: an additional absorption and a soft excess. The local additional absorption is usually 
attri buted to a warm ab s orber, i.e. absorption by totally or par t ially i onised matter: see, 
e.g., iNetzer et al.\ (|2003|) : iBlustin et al.\ (|2005|) : iKrongold et al.\ (1200 7') for more informa- 
tion about warm absorbers in Seyfert galaxies . The soft excess can be caused by reflection 
on the ionised surface of the accretion disc dCrummy et a/.ll2006h . Alternative explana- 
tions for the soft excess are black-body radiation or complex partially ionised absorption 
(se e Section 1 1.4.31 for brief discuss ion of the origin of the soft excess). 



Brenneman & Reynolds] (|2006|) included warm absorber and ionised reflection in their 



model and they showed that the extreme relativist i c smea ring is still required. However, 
in a spectral analysis by Miller. Turner & Reevesl (|2008ar) . the spectrum is characterised 
by a complex partial covering absorption in four different zones which affects the higher 
energy band more significantly implying that even no relativistic smearing was needed in 
the model. Our aim is not to discuss the appropriateness of the particular models. Instead, 
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Figure 4.7: The X-ray spectrum of MCG-6-30-15 observed by the XMM-Newton satel- 
lite. The continuum is characterised by a power law with the photon index F = 1.9 
absorbed by Galactic gas matter along the line of sight with column density n^ = 
0.41 X 10^^ cm^. The residuals from the model are plotted in the bottom panel clearly re- 
vealing features of a local absorption and soft excess in the soft X-ray band, and a feature 
at around 6 keV which can be explained by the presence of a broad iron line. 



we relied on the model of a broad line which we used in Section 12.31 for a comparison 
of two iron line models, kyrline and laor, and in Section I3.2.2[ for a comparison of 
different angular emissivity laws. 



Iron line complex 

The iron line complex (see Figure 1481) consists of three components: 



1. a narrow line at E = 6.4 keV possibly originating by reflection on a distant torus 
and/or outer parts of an extended accretion disc 

2. a broad relativistic line 

3. an additional narrow line which can be either an emission line E « 6.97 keV (of 
hydrogen-like ionised iron atoms), or an absorption line at E x 6.7 keV 

It is impossible to distinguish from the spectral analysis the statistical preference be- 
tween the models with the additional emission line of fully ionised iron atoms and with 
the absorption line of mildly ionised iron atoms. The values of the broad iron line pa- 
rameters obviously depend on the parameters of the underlying model. However, it was 
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Figure 4.8: A detailed view of the iron line band: the data residuals from the power 
law model are shown in terms of sigmas with error bars of size 1 cr. Energies of narrow 
components considered in the model are indicated by red lines: E = 6.4 keV (observed 
at £■ = 6.35 keV due to the cosmological redshift) of the neutral Kor emission line, and 
E = 6.77 keV (observed at E = 6.72 keV) of a mildly ionised absorption line. See the 
main text for further details. 



discussed by iFabian et al.\ (|2002h that the energy of the broad relativistic line is consistent 
with neutral iron Kor line originating in the innermost regions of cold disc independently 
of the i nclusion of the narrow ernission or absorption line. This was confirmed in the anal- 
ysis by Brenneman & Reynoldsl (|2006h who preferred the emission line at E x 6.97 keV 
becau se they found it to be more consistent with a Chandra observation (lYoung et al\ 
20051). However, they considered also the model with the absorption line. They found the 
equivalent w i dth of the absorption line EW = -21.3 eV much less than the one derived by 
Fabian et al.\ (|2002|) . where EW = -138 eV, and so, they concluded that it is not robustly 
required by the fit. 

We considered both cases and realised that the parameters of the broad iron line are 
more sensitive to the underlying model than it was suggested from the previous works, 
including the broad line energy, i.e. the ionisation state of the disc. 

The energy of the narrow absorption line was constrained to be £ = 6.77;^° IJ3 keV, 
see the left panel of Figure 14. 9[ where the dependence of the goodness of the fit is 
plotted against the energy of the absorption line. The photon index of the power law, 
spin, inclination, energy and radial emissivity of the broad iron line, and all normali- 
sations were allowed to vary during the fitting procedure. We measured the equivalent 
width of the absorption line EW = -(21 + 5) eV, consistently with the result obtained by 
Brenneman & Reynoldsl ( 2006 ). 



4.3. MCG -6-30-15 



89 




210 



200 




6.72 6.75 6.78 6.81 

Energy [keV] 



6.84 



6.84 6.87 6.9 6.93 6.96 6.99 
Energy [keV] 



Figure 4.9: Confidence plot of the energy value of a considered narrow absorption (left) 
or emission (right) spectral line. The dashed line represents 90% confidence level. 
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Figure 4.10: Contour plots of the energy of the broad iron line and inclination angle 
for two model alternatives. Left: narrow absorption line at E = 6.77 keV is included. 
Right: narrow emission line at E = 6.97 keV is included, instead. The individual curves 
correspond to 1 cr, 2 cr, and 3 cr, respectively. The best fit values are indicated with a 
small cross (in the right panel, the energy has the lowest possible value, E = 6.4 keV, so 
the cross is situated in the graph boundary). 



The alternative model employs the emission line, instead. The energy of the nar- 
row emission line was constrained to be £■ = 6.88 ± 0.02 keV, a rather lower value than 



6.97 k eV considered by lFabian et a/.l(|2002|) . but consistently with iBrenneman & Reynolds 
(|2006|) . The total goodness of the fit against the energy of the emission line is shown in 



the right panel of Figur e|4.9[ The equivalent width is E W = (28 + 4) eV, i.e. about twice 
larger than the value by lBrenneman & Reynoldsl(|2006h . 

Further, we investigated how the energy of the broad iron line depends on the choice 
of one of the model alternatives. We constructed contour plots for the energy of the broad 
line and inclination angle for the two model alternatives - see Figure l4?T0l The parameters 
of the narrow lines modelled as Gaussian profiles were fixed. We used E = 6.77 keV for 
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the energy of the absorption line, and E = 6.97 keV for the energy of the emission line. 
The photon index of the power law, spin, inclination, energy and radial emissivity of the 
broad iron line, and all normalisations were allowed to vary during the fitting procedure. 
The values for the broad line energy are relatively consistent, similarly depending on 
the inclination angle value. Nevertheless, in the alternative with the absorption line (left 
panel of Figure |4~T0l ), the ionised case is preferred rather than the case of a cold disc. 



4.4 IRAS 05078+1626 



IRAS 05078+1626 is a nearby Seyfert 1.5 galaxy. Before its identification as an infrared 
source it was also known as CSV 6150 (Catalogue of Suspected Variables). Its position on 
the sky is / = 186.1 and b = - 13.5 in the Galactic coordinates. The cosmological redshift 
of this galaxy is z « 0.018 ( "Takata et al.\\l994 ). It had never been spectroscopically 
examined in X-ray prior to the observation discussed in this Thesis. However, it was 
detected in a number of X-ray surveys, such a s the all- sky monitoring of the INTEGRAL 



IBIS/ ISGRI instrument (Sazonov et a/.ll2007h . the SWIFT BA T instrument (lAiello et al. 



20081 : iTueller et a/.ll2008h . and the RXTE Slew Survey (XSS) (iRevnivtsev et q/.II2004 . 

The X-ray spectroscopic properties of intermediate Seyferts are rather elusive: both 
obscu red Type 1 and unobscured Type 2 active galac t ic nuclei (AGN) have been reported 
(e.g., ICappi et all l2006l : iBrightman & Nandral l2008l : iBianchi et all l2009bh . It has been 
suggested that intermediate Seyfert galaxies are seen at intermediate inclination angles 
between pure "face-on" Seyfert Is and pure "edge-on" Seyf ert 2s, which follows di 



rectly frorn the orientation-based AGN unification scenarios ( Antonucci & Milleiill985 



Antonuccilll993l : IUrry & Padovanilll995|) . For this reason. X-ray spectroscopy of type 1.5 
Seyferts may provide clues to the nature and geometrical distributions of optically thick 
gas surrounding the active nucleus, the latter being the fundamental ingredient behind the 
unification scenarios. 

IRAS 05078-1-1626 is i ncluded in the FERO project ("Finding Extreme Relativistic 
Objects"; i Longinotti et al.WlOOS) with the aim of establishing the fraction of a relativis- 
tically broadened Ka iron lines in the spectrum of a complete flux-limited sample (2- 
10keVflux> ImCrab). 



4.4.1 Results from XMM-Newton observation in 2007 

Observations and data reduction 

The XMM-Newton observation of IRAS 05078-1-1626 was performed between 2007 Au- 
gust 21 UT 22:24:49 and 22 UT 15:35:43 (Obs. #0502090501). The EPN and both MOS 
cameras were operating in the "Small Window" mo de (see Section|4.1.1|). T he RGS cam- 
eras were operating in the "Spectroscopic" mode Men Herder et a/.ll200lh . The spectra 
were reduced with the SAS software version 9.0.0 (|Gabriel et q/.II2004|) . Intervals of high 
particle background were removed by applying count rate thresholds on the field-of-view 
(EPIC, single events) and CCD#9 (RGS) light curves of 0.35 cts/s for the PN, 0.5 cts/s for 
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Figure 4.11: EPIC-PN light curves in the 0.2-2 keV band (upper panel) and 2-10 keV 
band (middle). The hardness ratio HR is defined as the ratio of the counts at 2-10 keV to 
the counts at 0.2-2 keV and presented as a function of time in the lower panel. The bin 
time is as 2048 s. 



the MOS and 0.15 cts/s for the RGS. The exposure time after data screening is ~ 56 ks for 
MOS, ~ 40 ks for PN and « 58 ks for RGS, respectively. The patterns 0-12 were used for 
both MOS cameras, and patterns 0-4 (i.e. single and double events) for the PN camera. 
The source spectra were extracted from a circle of 40 arcsec in radius defined around the 
centroid position with the background taken from an off"set position close to the source. 
The two MOS spectra and the related response files were joined into a single spectrum and 
response matrix. Finally, the PN and MOS spectra were re-binned in order to have at least 
25 counts per bin and to over- sample the energy resolution of the instrument maximally 
by a factor of three, while the RGS spectra were left un-binned. Consequently, different 
statistics were used in fitting the s pectra - the traditional x^ statistics to fit the PN and 
MOS spectra and the C-statistics (|Cashlll976|) for all fits including RGS data. For the 
spectral analysis, we used XSPEC version 12.5, which is part of the HEASOFT software 
package version 6.6. 



Timing properties 



The PN light curve of the source is shown in Figure I4.11[ We have divided the energy 
range into two bands and checked the light curve behaviour in each of them, as well as 
a hardness ratio, which we defined as the ratio of the counts at 2-10 keV to the counts 
at 0.2-2 keV. The energy ranges were chosen for sampling different spectral components, 
as indicated by the energy where the continuum starts deviating from a power law model 
that describes the hard X-ray spectrum (see Sect. 14.4.11) . The hardness ratio stays almost 
constant during the observation, suggesting that no significant spectral variations occur. 
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Figure 4.12: PN spectrum extracted from the first half of the observation with the lower 
source flux (blue) and from the second half of the observation with the higher source flux 
(black). The ratio of the two spectra is presented in the lower panel. 

although the source flux increased by around 20%. 

To confirm this conclusion also for narrow spectral features, such as the iron emission 
line, we compared the PN spectra extracted during the first and the second halves of 
the observation (see Figure |4~T21) . The spectra correspond to the lower/higher source flux 
because the flux is increasing nearly monotonically during the observation. We calculated 
the ratio values of the two data sets and fit them with a simple function f(E) = aE + b 
using the least square method. The fitting results are a = -0.004 + 0.003 and b = 1.12 + 
0.01 with the sum of the residuals x^ = 111 for 190 degrees of freedom. When we set 
a = 0, the fitting results are comparably good with b = l.ll + 0.01 and;\f^ = 113. The 
ratio of the spectra is plotted in the lower panel of Figure |4.12[ Because no significant 
spectral diff'erences are evidentl, we analyse the time-averaged spectra hereafter. 



Mean spectral properties 



The signal-to-noise ratio is very good up to high energies (Figure B~131 ), so we fit the EPIC 
spectra spectra in the full energy range where they are well calibrated (0.35-12 keV). The 
X-ray continuum is described by a power law model at energies above 2 keV, although 
the iron line at E = 6.4 keV is present (Figure 14.141 ). The photon index of the power law 
is F ^ 1.49(1). In this and all subsequent models we included absorption by Galactic gas 
matter along the line of sight with column density n^ = 0.188 x lO^^cm"^. This value 
is from the Leiden / Argentine/Bonn HI measurem ents (IKalberla et Q/.ll2005h . We used 
the TBABS model (| Wilms. Allen & McCrayll2000|) to fit the absorption produced by the 



^Meanwhile, Ide Marco et all (120091) reported study of variability of FERO-sample galaxies, including 
IRAS 05078+1626. They also found that the spectral variability of this source is diminutive. 
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Figure 4.13: Total PN spectrum with the background level showing that the signal-to- 
noise ratio is very good up to high energies. 



Galactic interstellar matter. 

We applied the simple TBABS*POWERLAW model to both PN and MOS spectra. The 
X^ value is 3557 with 528 degrees of freedom ix^ Iv = 6.7) in the 0.35 - 12.0 keV energy 
range. The spectra differ from the power law model not only around E = 6.4 keV but 
also at lower energy band 0.35 - 2.0 keV (adding a Gaussian line model to fit the iron line 
improves the fit only to;^^/v = 3384/524 = 6.5). Residuals against this model are shown 
in Figure HTTSl They may be due to a warm absorber and/or soft excess. 

The spectral residuals reveal certain discrepancies between the PN and MOS spectra 
(see Figure |4?T5] with the data/model ratios of both spectra, PN and MOS, with the iden- 
tical model parameters of the spectra). The level of discrepancy is, however, comparable 
to th e level of systernatic u ncertainties in the cross-calibration between the EPIC cam- 
eras (IStuhlinger et a/.l 1200 6*). Nonetheless, we conservatively analyse the EPIC spectra 
separately. We use the same models for both spectra but allow the values of the model 
parameters to be different. The values of the photon index using the simple power law 
model differ from each other when fitting the spectra independently, resulting in a harder 
PN spectrum with F = 1.60(1) compared to the MOS spectrum with F = 1.54(1), ignor- 
ing the energies below 2keV and also between 5.5-7.5 keV. Although the absolute value 
of these spectral index measurements does not have a direct physical meaning, given the 
simplicity of the model applied on a small energy band, the comparison between them is 
illustrative of the quality of the cross-calibration between the EPIC cameras. Differences 
of the order o f AF ^ 0.06 in the har d X-ray band are consistent with current systematic 
uncertainties (IStuhlinger et a/.ll2006|) . 
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Figure 4.14: XMM-Newton PN (black) and joint MOS (red) spectrum of 
IRAS05078+1626 described by a simple power law model absorbed by Galactic neu- 
tral hydrogen in the line of sight with hh = 0.188 x 10^^ cm"^. The photon index of the 
power law is F = 1.49. 



RGS spectrum 

We jointly fit the un-binned first-order spectra of the two RGS cameras with the same 
model's parameter values except the overall normalisations. The continuum is well-fitted 
by the simple power law model with the photon index F = 1.57. We searched further 
for narrow emission and absorption lines in the spectrum using several ZGAUSS models 
with the intrinsic width cr set to zero. We calculated the 90 % confidence interval for 
a blind search, as P = Po/Mriah where A/^triai = A/^bins/3 = 3400/3 and Pq = 0.1. For 
the RGS data P = 8.8 x 10"^, to which AC = 22.4 corresponds assuming the Student 
probability distribution. The only line fulfilling this criterion by improving the fit about 
AC = 31.7 is an emission line at the energy E = 0.561 + 0.001 keV (22.10 + 0.04 A) and 
the equivalent width 7^3 eV. We identified it with the forbidden line of the O VII triplet 
(£LAB=0.561keV). 



EPIC spectrum 



The forbidden line of the OVII triplet is clear signature of a photoionised plasma. No 
significant features were detected that may be expected alongside the O VII (f) line, if it 
were produced in a collisionally ionised plasma, such as the resonance line in the OVII 
triplet or the OVIII Lyor. This led us to try and explain the residuals against a power law 
model in the soft X-rays as eff"ect of intervening ionised absorption gas. We used the 



4.4. IRAS 05078+1626 



95 




1.6 



0.9 1.2 

Energy [keV] 



1.4 ■ 



1.2 ■ 



0.8 



i 



-f 



tt 



It 



'k 



-im 1)% 



5.6 



6 6.4 

Energy [keV] 



6.8 



Figure 4.15: Ratios of the simple power law model (the same as in Figure l4T4l) to the 
data in different energy bands: left: at lower energies, right: in the iron line band, where 
the narrow Kor line at the rest energy E = 6.4 keV is prominent (observed at £ = 6.29 keV 
due to the cosmological redshift). Black crosses correspond to the PN data points, while 
the red crosses correspond to the MOS data points. 



XSTAR model version 2.11n7c (|Kallman & Bautistall200l|i ^l to calculate a grid of tabular 
models with the input parameters constrained from the preliminary data analysis with sim- 
ple models whenever possible (photon index Y ^ 1.7, density p < lO^'* cm"^, luminosity 
L < lO'^'^ergs"^ column density 10'^ cm"^ < ^h < 10^^ cm"^, and ionisation parameter 
-5<log^<5). 

A single-zone warm absorber component modifying the power law continuum dra- 
matically improved the fit from;tVv = 1850/270 = 6.9 to ;tVv = 402/270 = 1.5 for 
the PN spectrum. The ionisation parameter converged to a very low value, and we found 
that this almost neutral absorption can be successfully reproduced with zphabs, which 
is a simpler model than XSTAR, so we preferred this possibility. The addition of another 
warm absorber zone improves the fit to;t'^/v = 320/266 =1.2 for the PN spectrum, and 
it requires the ionisation parameter log^ = 3.9. We checked that adding another warm 
absorber zones does not improve the fit significantly. The residuals from the model (see 
Figure 14361 upper panel) reveal an extra emission that remains at low energies, as well as 
around the iron Kor line band. 

These features can come from reflection of the primary radiation on the surface of 
the accretion disc, so we added the reflion model (Ross & Fabianll2005 o. which calcu- 
lates the ionised reflection for an optically thick atmosphere with constant density. We 
examined the significance of the addition of the reflection component into the complex 
model of the PN spectrum by the statistical F-test. The low value of the F-test probabil- 
ity (5 X 10"^^) strongly favours this additional model component. The best fit was now 
X^/v = 246/264 = 0.95 for the PN spectrum. We hereafter call this model the 'baseline' 
model; in theXSPEC notation: tbabs x zphabSn x xstar x (powerlaw -i- reflion 

-I- ZGAUSS). 

The parameter values of the 'baseline' model are presented in the Table |4.5[ The 



http://heasarc.gsfc.nasa.gov/docs/software/xstar/xstar.html 



Table 4.5: Parameters of the 'baseline' and 'final' models. 



OS 



Model 


Model 


'baseline' model 


'final' 


('double reflection') model 


component 


parameter 


PN 


MOS 


PN 


MOS 


PN+MOS+RGS 


ZPHABS 


«h[1022 cm-2] 


n 1^4+0.005 
"•^"^-0.007 


n 1 9Q+0.007 
"•^^^-0.007 


102+"""'' 

"•^"^-0.005 


120+"""*^ 


i06+"""4 


XSTAR 


«H[1022cm^2j 


130+2" 


170+20 


120+^" 


150+™ 


130+2" 

IJU 20 




log^ 


o O+O.l 

-^•-'-0.1 


2 4+''-i 


2 2+^-4 


2 5+1-" 


2 5+1" 




He/Hesolar - Ca/Casolar 


1(f) 


1(f) 


1(f) 


1(f) 


1(f) 




Fe/Fesolar - Ni/Nisolar 


o-2!!ii 


l+^i 


1 2+"-^ 

^•^-0.3 


9+"-2 

"•^-0.2 


1 l+"-2 
^•^-0.2 


POWERLAW 


r 


1 81+0.03 


1 on+005 

^•^0-0.05 


1 75+""' 

^•'-'-0.03 


1 74+0.07 
^•'^-0.03 


1 76+""4 

^•'°-0.02 




normalisation 


(7 ± 1) X 10-4 


(7 ± 1) X 10-4 


(6 ± 1) X 10-4 


(7 + 2) X 10-4 




REFLION 


r 


1.81 (b) 


1.80(b) 


1.75 (b) 


1.74(b) 


1.76 (b) 




log^ 


3 O^''^ 

-'•"-0.2 


3 2+"-2 

•'•^-0.2 


3 0+"-2 

J.U_Q2 


3-1 


3 0+"-i 




Fe/Fesolar - Ni/Nisolar 


0.2 (b) 


0.1 (b) 


1.2 (b) 


0.9 (b) 


1.1 (b) 




normalisation 


(3 ± 2) X 10-'' 


(3 ± 2) X 10-'' 


(2 ± 1) X 10-^ 


(1 ± 1) X 10-'' 




REFLION 2 


normalisation 


- 


- 


3 + 1 X 10"' 


4 ± 1 X lO-'^ 




ZGAUSS 


E [keV] 


6 40+"*" 


6 44+""2 

"•^^-0.02 


- 


- 


- 




(T [keV] 


06+""^ 


02+''"4 

"•" -0.02 


- 


- 


- 




z 


0.018 (f) 


0.018 (f) 


- 


- 


- 




normalisation 


(3 ± 1) X 10-^ 


(3 ± 1) X 10-^ 


- 


- 


- 


x^ly 




246/264 


405/243 


256/266 


404/244 


C/v = 1551/1347 



Note: The sign (f) after a value means that it was fixed during the fitting procedure. The sign (b) means that the parameter was bound to the value 
of the corresponding parameter of the previous model component. The sign "-" means that the model component was not included while dots only 
mean that there are more values related to the individual spectra which are not necessary to be all shown in the table. 
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Figure 4.16: Residuals of the PN data from the 'baseline' model including both ionised 
reflection and absorption {lower), only reflection {middle), and only absorption {upper). 



quoted errors of the parameters represent a 90 % confidence level for a single interest- 
ing parameter. The measurement is obviously affected by a much larger systematic error, 
which, however, could be properly quantified only if we knew the "right" model. The 
value of the power law photon index increased to F a: 1.8 compared to the simple model 
applied to the data in Sect. 14.4.11 because we included the additional local absorption in 
the model. The data residuals from the model are shown in the lower panel of Figure l4~T6l 
In the same figure, we also show residuals from the best fit performed with the 'baseline' 
model, excluding the ionised absorption (middle panel) and the ionised reflection compo- 
nent (upper panel). 

The narrow iron Ka line with the rest energy E = 6.40 + 0.01 keV, the width cr = 
0.06 + 0.03 keV, and the equivalent width EW = 82 + 15eV evidently represents cold 
reflection. This suggests an origin of this spectral component in the outer part of the 
disc, or from the torus. The cold reflection is also supposed to contribute to the soft 
part of the spectrum with the individual emission lines. For this reason, we replaced 
the Gaussian profile in the 'baseline' model with another reflion component (called as 
REFLION 2 in the Table 14.51) with the same values for the photon index and abundances 
as the REFLION 1 model component. The ionisation parameter was kept free during the 
fitting procedure, but it very quickly converged to its lowest value ^ = 30 (log^ = 1.477). 
The advantage of the REFLION model compared to the other available reflection models 
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Figure 4.17: The 'final' model. The total model is shown in black (solid line), the primary 
radiation is red (dashed), the reflion components are blue (dotted) for cold reflection 
and magenta (dot-dashed) for ionised reflection. 



is that it also includes the soft X-ray lines, with the disadvantage in this case that the 
ionisation parameter cannot be set to zero. 

This 'double reflection' model, in the XSPEC notation tbabs x zphabSn x xstar 
X (POWERLAW -I- REFLION -I- REFLION), docs not significantly improve the fit goodness 
over the 'baseline' model (with x^/v = 256/265 = 0.96 for the PN spectrum), but it 
represents a more self-consistent astrophysical picture. Therefore, we call the 'double 
reflection' model as 'final' model. In contrast to the 'baseline' model, it does not require 
sub-solar iron abundances, see the Table |431 where the parameter values for this model 
are presented. The 'final' model with each component separately drawn is shown in 
Figure I4.17[ All the plotted components are absorbed by a warm absorber surrounding 
the central accretion disc and two kinds of cold absorber - one from Galactic interstellar 
matter and one from local absorber in the host galaxy. 

Some model parameters were not allowed to vary during the fitting procedure. The 
redshift of the ionised absorber was fixed to the source cosmological value, because leav- 
ing it free yields a negligible improvement in the quality of the fit. Second, we used the 
same iron abundances across all the components in the model. 

In the 'final' model, the warm absorber ionisation parameter is consistent with the 
ionised reflection component. This result is also presented in Figure |4.18[ where the 
contour lines related to the Icr, 2cr, and 3cr levels of ;^^^ between the ionisation parameters 
of the two model components are presented. 

Table 14.6! summarises flux values of the individual components of the 'final' model 
for both PN and MOS spectra for two energy bands, 0.5 - 2keV and 2 - lOkeV, and 
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Figure 4.18: The contour plot of the ionisation parameters of the REFLION model, rep- 
resenting the ionised accretion disc, and of the XSTAR model, representing the warm 
absorber in the 'final' model. The individual curves correspond to the 1 cr, 2 cr, and 3 a 
levels. The position of the minimal value of ;^^ found by the fitting procedure is marked 
by a cross. The corresponding^^ values are given in the plot. 



also shows fractions of the reflection radiation to the total emission (sum of the primary 
and reprocessed radiation). The flux ratio is almost equally shared between the cold and 
ionised reflection components, and its value is in total 7? < 0.2. The absorption-corrected 
luminosity values of the source in the same energy bands are L()5_2kev = (1-21 + 0.02) x 
10"^^ erg s"^ and L2-iokev = (1-87 + 0.02) x 10^^ erg s~\ respectively. 

We also used the 'final' model for a joint fitting of all the XMM-Newton instruments 
(PN, MOS and both RGS spectra) together. The parameter values were bound among all 
the spectra, only normalisation factors were allowed to vary. The goodness of the joint 
fit is given in C-statistics because the RGS data are un-binned (and each individual bin 
contains only a few counts). The result is C = 1551 for a number of degrees of freedom 
V = 1347. All the spectra, together with the residuals, are shown in Figure BTTyi and the 
corresponding parameters in the last column of Table 14.51 



4.4.2 Discussion of the results 

Constraints on the location of the absorbers 

In this section, we discuss a possible location of the absorber's system in the ' final' model 
Photoelectric absorption is almost invariably observed in Type 2 Seyferts (|Awaki et al 



19911 : [Turner g?a/.l 1 19971 : iRisaliti. Elvis & Nicastrd 120021) and generally attributed to an 



optically thick m atter responsible for orientation-dependent classification in AGN unifi- 
cation scenarios ( Antonucci & Milleiil 19851: lAntonuccill 1993 ). 
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Table 4.6: Flux values of the 'final' model and its individual components. 



Model 


Flux at 0.5- 2 keV 


Flux at 2 


- lOkeV 


component 


[lO-i^erg 


cm ^ s '] 


[10-'2erg 


cm~^s~'] 




PN 


MOS 


PN 


MOS 


total model 


7.05!-^ 


7 00+°°^ 

' -""-0.03 


25.0!°:^ 


25.4!- 


unabsorbed model 


I6.6:li 


16 7+°-^ 

^"•'-0.2 


25 7+°-2 


26.5_+°:^ 


POWERLAW 


13 6+°-' 


14 3+°' 


22 6+°-^ 


23.8_+°:} 


REFLIONion 


1 9+0.2 


1 1+°' 


1 6+°' 

^•"-0.2 


o.9!°:i 


REFLIONcoid 


1 1+°1 


1 3+°' 

^•-'-0.1 


1 5+°' 


1.8!°:? 




12+°°' 

"•^^-0.01 


0.07:°;°! 


06+°°' 

"■""-0.01 


03+°°' 

"•"-'-0.01 


■^cold * 


07+°°' 

"•"'-0.01 


0.08t°I 


06+°°' 

"•""-0.01 


07+°°' 

"•"'-0.01 



* the ratios of the reflection component flux values to the flux value of the total unabsorbed 
model (sum of the primary and reflected radiation). 



Because the IRAS 05078+1626 galaxy is probably viewed under an intermediate in- 
clination between unobscured Seyfert Is and obscured Seyfert 2s, the torus rim may also 
intercept the line of sight to the AGN and absorb the radiation coming from the centre. 
The col d absorption can, however, als o be associated with the interstellar matter of the 
galaxy (ILamastra. Perola & Mattll2006h . 

Both reflection components are inside the ionised absorber in the 'final' model. The 
geometrical interpretation is that the cold reflection occurs on the outer parts of the disc 
or the inner wall of the torus. Reflection on the nearer part of the torus is heavily absorbed 
by the torus itself, so only radiation reflected on the farther peripheral part of the torus 
can reach the observer after passing through the warm absorber. However, an alternative 
scenario, in which the cold reflection is unafl"ected by the warm absorber, i.e., tbabs x 

ZPHABS X [REFLIONcoid + XSTAR X (POWERLAW + REFLIONdisc )], is also acceptable 

with;tVv = 265/265. 



The lack of constraints on the variability in the warm absorbed features (JKrongold et al. 



2007h . caused by the moderate dynamical range of the primary continuum, as well as sta- 
tistical limitations in our spectra, prevents us from precisely constraining the location of 
the warm absorber. 



Constraints on the location of the ionised reflector 

The ionised reflection might occur either at the inner wall of a warm absorber cone or on 
the accretion disc. Even in the latter case, the reflection cannot occur arbitrarily close to 
the black hole. In this section we investigate the constraints of the accretion disc location 
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Figure 4.19: The joint fit of all spectra of the XMM-Newton instruments - PN (black), 
MOS (red), RGS 1 (magenta), and RGS 2 (blue), together with the model residuals. 



and structure, which can be drawn from the lack of the significant relativistic blurring of 
the disc reflection component. 

We convolved the ionised reflection component with the fully relativistic kyconv 
model (|Dovciak. Karas & Yaqoobll2004l) . Two assumptions about the disc emissivity were 
considered. First, the radial part of the intensity decreases with the power of the disc ra- 
dius q{I oc r"*), where the value of ^ was allowed to vary between 2 and 3.5. Second, the 
angular dependence was assumed to be isotropic which seems to be appropriate approxi- 
mation for our situation of an X-ray irradiated accretion disc (Svoboda et al. 2009) . 

We examined the expected confidence levels of the best-fit values of the disc's inner 
radius by stepping this parameter in the whole range of its possible values - from the 
horizon to the outer disc radius, which we set to 400 gravitational radii {Rg = GM/c^). 
The results are shown in Figure l4~20l At the 90 % confidence level, the accretion disc is 
not allowed to extend closer to the black hole than 60Ra. 

The "relativistic blurring method" would be less appropriate in looking for the im- 
prints of the innermost parts of the accretion disc if the disc w ere too highly ionise d 
(log ^ « 4) and the narrow reflection features were not present (|Ross & Fabianll2005|) . 
However, the ionisation parameter value of the reflection component is not so high in 
the 'final' model and the dominant feature is the intermediately ionised iron line (E « 
6.7 keV). If we assumed a stratified disc with the ionisation state decreasing with the ra- 
dius from the centre, the hydrogen-like iron line would be also expected to appear in the 
spectrum (as an intermediate stage between the over-ionised and mildly ionised contribu- 
tion). Because it is not detected in the data, the accretion disc truncation provides a more 
reasonable explanation of missing signatures of the relativistic blurring. 
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Figure 4.20: The best-fit values of ;t'^ statistics for the inner disc radius parameter, which 
we obtained by gradually stepping it from the horizon radius to the outer radius of the 
disc (400i?g). The dashed line is the 90% confidence level for one interesting parameter. 



Mass accretion rate 

Disc truncation is ex pected in low-luminosity AGN where the inn er accreti on flow is 



advection-dominated (INarayan & Yi 1994 : Esin. McClintock & Narayan 1997l : lNarayan & McClintock 
2008L and references therein). The transition fro m the outer standard accretion disc may 
occur, e.g., via the disc evaporation mechanism (JMeyer. Liu & Meyer-HofmeistenllOOOl : 
Liu & TaamI |2009|) . The observationa l evidence of a tru n cated accretion disc in low- 
luminosity AGN was reported e.g. by Lasota et al.\ (119961) : lOuataert et al.\ (119991). How- 
ever, i ts presence is also suggested in some observations of Seyfert galaxies (ILu & Wang 
2000l : bone. Madejski & Zvcki|l2000l : IChiang & Blaesll2003l : iMarkowitz & Reevesll2009[) 



and even a quasar (JMatt et al. 



20051) where the luminosity value is estimated as a half 



of the Eddington value. General ly, it is expected that t he lower the luminosity, L/L^m, 
the larger transition radius (see lYuan & NarayanI 12004 . and references therein). Fur- 
thermore, we investigate whether the disc truncation hypothesis is consistent with the 
IRAS 05078-1-1626 luminosity. To have these quantities in Eddington units, we first esti- 
mated the mass of the black hole. 

IRAS 05078-1-1626 belongs to the s ample of the infrared-se lected Seyfert 1.5 galaxies 
observed by a 2. 16 m optical telescope (i Wang. Wei & Hell20 06^ where, among others, the 
velocity dispersion in the O III emission line was measured. The correlati on between the 
O II I line width and the mass of active galactic nucleus was discussed in iNelsonI (|2000[ ) 
and iBorosonI (|2003|) . The value from the optical measurements, crnm ~ 130km s"\ 
corresponds to the mass M « 4 x IO'Mq using a correlation plot in lBorosonI (|2003b . The 



scatter of the correlation is somewhat large with the reported limit of a factor of 5 for an 
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uncertainty in the black hole mass determination, so the value only provides an order of 
magnitude estimation. 

The value of the Eddington luminosity is Lgdd - 1-3 x lO^^M/MoCrg s"^ « 5 x 
10"^^ ergs'' for the giv en value of the mass. We used luminosity-dependent corrections by 



Marconi et al\ (|2004|) to estimate the bolometric luminosity of IRAS 05078+1626 from 
the X-ray luminosity. Its value is L a: 5 x 10'*'*ergs"' ~ lO'^Lndd- Correspondingly, the 
mass accretion rate, M = L/c^, is sub-Eddington with M » O-IMphh- Th is value is typical 
of less luminous Seyfert galaxies (see for example iBianchi et a/.ll2009ah . and is consistent 
with the disc truncation hypothesis. 



Chapter 5 
Conclusions 



We employed the relativistic models of iron line to probe the innermost regions of black 
hole accretion discs with the current X-ray data by XMM-Newton, and with the simulated 
data for next generation X-ray missions. We found the modelling of broad iron lines to be 
a suitable method to measure the angular momentum of black holes at all scales - from 
stellar-mass microquasars to giant black holes of billions solar masses of distant quasars. 
In this section, we summarise the main conclusions of our investigation presented in the 
previous chapters. 

5.1 Relativistic line models 

In Section [21 we described how the observed profile of an intrinsically narrow emission 
line is distorted due to the effects of rapid orbital motion and strong gravity. We inves- 
tigated the iron line band of the X-ray spectrum for two representative sources - MCG - 
6-30-15 (active galaxy) and GX 339-4 (X-ray binary). The iron line is statistically better 
constrained in the active galaxy MCG -6-30-15 due to a significantly longer exposure 
time of the available observations. The spectra of both sources are well described by a 
continuum model plus a broad iron line model. We determined the spin values using the 
KYRLINE model as: 

• MCG-6-30-15:a = 0.9-1.0, 

• GX 339-4: a = 0.56 - 0.85. 



T he value for MCG -6-30- 15 is co nsistent with previously obtained results (|Fabian et al. 



20021 : iBrenneman & Reynoldsl 120061). However , in the case of GX 3 39-4, our best-fit 



value is lower than the one by iMiller et al.\ (12004) and lReis et al.\ (|2008b . The model with 



their derived values does not provide an acceptable fit when applied to the data re-binned 
with respect to the instrumental energy resolution. We found that the iron line is not as 
extremely broad as previous analyses suggested, and the spin value has rather an interme- 
diate value a ~ 0.1. 

We compared two relativistic models of the broad iron line, LAOR and KYRLINE. 
In contrast to LAOR, the KYRLINE model has the spin value as a variable parameter. 
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However, the laor model can still be used for evaluation of the spin if one identifies the 
inner edge of the disc with the marginally stable orbit. We realised that the discrepancies 
in the results between the KYRLINE and laor models are within general uncertainties of 
the spin determination using the skewed line profile when applied to the current data (for 
MCG -6-30-15: aiaor = 0.94-0.9982, and for GX 339-4: aiaor = 0.63-0.87). This means 
that the spin is currently determined entirely from the position of the marginally stable 
orbit, as it is done with the LAOR model. 

However, the results are apparently distinguishable for higher quality data, as those 
simulated for the next generation X-ray missions, which will be sufficiently sensitive to 
resolve the slight variations in the overall line shape due to the spin. We found that the 
LAOR model tends to over-estimate the spin value and furthermore, it has insufficient 
energy resolution which affects the correct determination of the high-energy edge of the 
broad line. The discrepancies in the overall shape of the line by the LAOR model are more 
visible especially for lower values of the spin. 

Another advantage of the kyrline model over the LAOR model is that the KYR- 
LINE model gives better pronounced minima of ;^f^ for the best-fit values. The confidence 
contour plots for a/M versus other model parameters are more regularly shaped. This 
indicates that the kyrline model has a smoother adjustment between the different points 
in the parameter space allowing for more reliable constraints on spin. 

5.2 Emission directionality 

In Section [31 we investigated how the black hole spin measurements using X-ray spec- 
troscopy of relativistically broadened lines depend on the employed definition of the an- 
gular distribution of the disc emission. We considered three different approximations of 
the angular profiles, representing limb-brightening, isotropic and limb-darkening emis- 
sion profiles. We studied the role of the emission directionality in the spin determination 
with the current data - we used the XMM-Newton observation of MCG -6-30-15, and 
with the artificial data simulated with energy resolution and sensitivity of future X-ray 
missions, such as IXO (International X-ray Observatory). 

We realised that the use of an improper directionality profile could affect the other 
parameters inferred from the relativistic broad line model. Especially for the frequently 
used case of limb darkening, the radial steepness can interfere with the line parameters of 
the best-fit model by enhancing the red wing of the line. The limb darkening law favours 
higher values of spin and/or steeper radial dependence of the line emissivity; vice-versa 
for the limb brightening profile. 

Steeper radial emissivity in the innermost region of an accretion disc has b een de- 



tected in several other sources, both in AGN and black hole binaries (lMillenl2OO70 . The 
improperly used limb darkening in the reflection model represents one of the possible ex- 
planations. For this reason, constraining the "correct" directionality in more sophisticated 
models of reflection spectra of accretion discs is desirable. 

We applied the NOAR radiation transfer code to achieve a self-consistent simulation 
of the outgoing spectrum reflected from a cold disc without imposing any ad hoc formula 
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for the emission angular distribution. We found that the isotropic directionality reproduces 
the simulated data to the best precision. On the other hand, the model with the limb- 
darkening profile was not able to reproduce the simulated data especially for the higher 
values of black hole spin. 

We consider this to be an important result since much of the recent work on the iron 
lines has revealed a significant relativistic broadening near rapidly rotating central black 
holes. In some of these works, the limb darkening law was employed without testing 
different options. The modelled broad lines are typically characterised by a steep power 
law in the radial part of the intensity across the inner region of the accretion disc, as in the 
mentioned MCG -6-30-15 observation. We conclude that the significant steepness of the 
radial part of the intensity also persists in our analysis, however, the exact values depend 
partly on the assumed angular distribution of the emissivity of the reflected radiation. 



5.3 Data reduction 

We re-analysed XMM-Newton observations of MCG -6-30- 15 and GX 339-4. According 
to the previous X-ray spectroscopy analyses, these sources exhibited extremely broad 
iron lines. The results of our spectral analysis are discussed in the context of the aspects 
of iron line modelling in Sections [IE We stressed that the correct re-binning of the 
data, which reflects the energy resolution of the used instrument, is necessary to obtain 
statistically relevant results. The un-grouped EPIC/PN or MOS data are characterised by 
large error bars which implies that any applied model produces lower x^ values when 
fitting such data (see Section 14.2.21 ). As a consequence, a larger set of difl^erent models 
may be accepted, including the incorrect models which would be excluded if the data 
were properly grouped. 

Photon pile-up is another problem which may aff"ect the spectral analysis significantly. 
It occurs during exposures of too bright sources, such as some Galactic black hole bina- 
ries. We encountered a serious problem with the pile-up in the data re-analysis of the 
low/hard state observations of GX 339-4 (see Section 14.2.11 ). We found that the elimi- 
nation of the pile-up by excising the core of the point spread function of the source is 
possible only at the cost of a drastic loss of counts which makes any interpretation of the 
spectra rather disputable. 



5.4 X-ray spectrum of IRAS 05078+1626 

In Section I4.4[ we presented a spectral analysis of the XMM-Newton observation of a 
Seyfert galaxy IRAS 05078-1-1626 being the first X-ray spectroscopy study of this source. 
The X-ray continuum spectrum of IRAS 05078-1-1626 is dominated by a power law 
with a standard value of the photon index (F = 1.75 in the 'final model'). The residuals 
from the power-law continuum can be interpreted in terms of a warm absorber surround- 
ing an accretion disc, and a reflection of the primary radiation from an ionised matter and 
on a cold torus. The outgoing radiation is absorbed by cold matter (hh « 1 x 10^' cm"^). 
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which can be either located in the inner side of the torus or caused by gas in the host 
galaxy. The type of the galaxy determined from the previous infrared and optical re- 
search is Seyfert 1.5, suggesting that the active nucleus could be seen at large inclination, 
consistent with either interpretation or even allowing a combination of both. 

The ionised warm absorber occurs in the central part of the AGN. Its column density 
was found to be rh > 1 x 10^"* cm"^, which i s a rather high value compared to the warm 



absorbers detected in other Seyfert galaxies (jBlustin et a/.l 120051) . This may be because 
we are looking through a longer optical path of a conical non-relativistic outflow due to 
the high inc lination of the system. Such a conical outflow is suggested in the model by 



Elvisl (I2000|) (see Figure [LTl) . The ionisation parameter of the warm absorber is log ^wa = 
2.5 + 1.0, which is comparable to the value related to the ionised reflection log ^reflection = 
3.0 + 0.2, suggesting a link between them. 

If the ionised reflection is associated to the warm absorber (e.g. the inner walls of 
a conical outflow), the lack of spectral features associated with the accretion disc is a 
natural consequence thereof. If, instead, the ionised reflection occurs at the accretion disc, 
it cannot extend up to the marginally stable orbit. The lack of the significant relativistic 
blurring of this model component requires the disc to be truncated (inner disc radius 
^in > 60 Rg). This idea is also supported by the low ratio of the reflection radiation to the 
primary one, R < 0.2, and also by the relatively low mass-accretion rate M « O.lMedd 
determined from the source luminosity. 



Chapter 6 

Future perspectives 



X-ray spectroscopy is a feasible method to probe the innermost regions of black hole 
accretion discs in active galactic nuclei and stellar-mass black hole binaries. Especially, 
it provides a unique way to measure black hole spin which plays an important role in 
black hole energetic balance and evolution. The existence of powerful relativistic jets 
is frequently attributed to rapidly rotating black holes where the energy extracted from 
the black hole rotation generates and maintains the far-reaching coUimated ejections of 
plasma. But is it really so? Is there any correlation between the spin and radio-loudness of 
galaxies and black hole binaries? What is the statistical distribution of the black hole spin? 
Is the black hole spin natal or are black holes spun up via accretion? Such questions could 
be answered if we reliably knew the values of the spin for a sufficiently large number of 
black holes. Hence, precise measurements of the black hole spin is among the important 
future tasks of X-ray astronomy. 



^lack hole angular momentum has been measured since late 1980s (|Fabian et al. 



19891). At the beginning with a large uncertainty, but following the fast development 



of X-ray detectors, these measurements have become increasingly precise. The origi- 
nal method dealt with a relativistically broadened iron line. This method is still one of 
the most suitable methods of the spin measurements. Its advantage is its applicability to 
black holes at all mass-scales. However, it depends on the accretion regime whether a 
broad relativistic line may occur in the spectrum or not. The necessary conditions are that 
the accretion flow is in the "thin-disc" regime (the disc is geometrically thin and optically 
thick) and not over-ionised (^ < 5000). 

The advantage of stellar-mass black hole binaries is that their spectral states evolve 
on relatively short time-scales so that it is possible to detect them in different states. A 
relativistically broadened iron line should be common fo r stellar-mass blac k hole bi- 



naries when observed in the appropriate accretion state (iMiller et al.\ |2009|) . On the 
other hand. X-ray spectra of the brightest Galactic black hole binaries, such as Cyg X- 
1, GRS 1915-1-105, GX 339-4, etc., suffer significantly from the pile-up when observed 
with current X-ray satellites. Instruments with a more frequent and faster read-out cycle 
or with a different detection technology would be more suitable for observations of such 
bright sources. 

For stellar-mass black hole binaries, the spin measurements with a relativistic iron 
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line may be compared to the res ults obtained by modelling of the disc thermal radiation 



(|Remillard & McClintockl 120061) . An agreement of the results obtained by both methods 
will enhance the credibility of the spin measurements. Higher resolution X-ray detectors 
will allow other methods, such as X-ray reverberation or X-ray polarimetry, to indepen- 
dently measure the spin. These will also be applied for active galaxies. 

Many similarities between super-massive black holes in active galaxies and Galactic 
black holes were reported (see Section [1.4.31) . Do they also have a similar distribution 
in the spin? A statistical study of the presence of relativi stically broadened iron lines in 



spectr a of active galactic nuclei (AGN) was performed bylGuain azzi. Bianchi & Dovciak 



(|2006l) on the sample of a hundred AGNs observed by the XMM-Newton satellite. They 
estimated a fraction of 42+ 12 % of well exposed (> 10"* cts) AGNs that e xhibit a relativis 



tically broadened iron line. A continuation of this effort is represented by lLonginotti et al. 



(|2008j and de la Calle et al. (20 09, submitted to ASlA). Some other wo rks have be en done 



on dif ferent samples of galaxies dStreblyanska et a/.l2005l : lNandra et a/..2007, : iBrenneman & Reynolds 



2009) 



All of these researches have been done on well exposed observations of active galaxies 
in the local Universe. However, with more sensitive instruments, such as those planned 
for IXO, the relativistic iron lines will be detectable also in observations of more distant 
galaxies. Constraints on cosmic density of relativistic lines and its distribution on cosmo- 
logical length-scales w ill be available in near future. A pioneering work with the current 



data appeared recently (IB allanty ndl20 1 0|) . 

It is worth-mentioning that, like any other spectroscopic results, the spin measure- 
ments with a relativistic iron line are model-dependent. Models employing partial cover- 
ing absorption represent an alter native explanation of t he characteristic spectral curvature 



around the iron line energy (Miller. Turner & Reevesll2008a) . Observations with future 



X-ray satellites may be helpful to distinguish between the two models thanks to an in- 
creased sensitivity in combination with a broader energy range including the Compton 
hump as well. Self-consistent reflection models are then desirable to model all the reflec- 
tion components together. If the source is considerably variable then the timing studies 
may be determinative between the two scenarios. The first reverberation me asurement of" 



1^7 07-495 strongly favours the scenario with the relativistic line models (IFabian et al. 



2009) 



Beyond the current sensitivity of X-ray detectors, further details of accretion physics 
can be studied. Not only spin will be measured by relativistic line models but also the 
intrinsic assumptions inside the models will be tested. One of the uncertainties in the 
relativistic line models is due to the unknown emission directionality. We discussed its 
effect on the spin value determination for the case of an isotropically illuminated cold 
accretion disc (see Section O. This analysis can be enhanced in a future work by adopt- 
ing different assumptions, e.g. ionised disc, or an illuminating source localised above 
the black hole. More comprehensive simulations of the irradiated accretion discs taking 
the general relativistic effects on the radiation transfer into account would be desirable. 
Detailed knowledge of the model assumptions will be necessary to properly measure the 
spin. 

Currently, the uncertainties in the precise position of the inner edge represent an un- 



no 
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known error of the spin determination. However, it appears that the expected magnitude 
of these errors does not prevent us from setting interesting and realistic constraints on 
the spin parameter. Our present treatment of the problem is incomplete by neglecting 
the magnetohydrodynamical effects and their influence on the ISCO location. Future im- 
provements in our theoretical understanding of the inner edge location are desirable and 
will help to improve the confidence in the determination of the model parameters. 

Higher resolution X-ray detectors will possibly reveal more complex profiles of iron 
lines modified by a contribution of different non -axial patterns, for instance orbi ting hot 
spots (|Dovciak et Q/..,2004b) or spiralling waves te aras. Martocchia & SubnbOOli) . Some 
of the current obseryations of AGNs alre ady revealed periodic X-ray variability (see e.g. 



Iwasawa et a/.lll998l : [Turner et al.\ |2002[) or narrow emission lines in the iron line band 



which may be due to orbiting hot spots (see e.g. iDovciak et Q/.ll2004bl : IPiconcelli et al. 
20061) ■ T he study of non-axial patterns in the iron line profiles is possible using the KY 



package (|Dovciak. Karas & Yaqoobll2004|) . and will be promising with sensitive future 
X-ray missions. 

We would like to conclude by expressing our optimism regarding the future of X-ray 
research of black holes and accretion discs. New technologies that are currently being 
developed will enhance energy and time resolution of the observations and they will even 
open new channels of information (such as X-ray polarimetry). In addition, theoretical 
models and numerical codes are getting more advanced, thus allowing us to make inter- 
pretation of the observed data closer to their real meaning. 
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